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Preface 


This document stipulates protocols for maasming b.o-opfica, and 

of-view Sensor (SeaWiFS) algorithm development and ” [ validation of the Sea- 

with specifications of the measured variables which are to be included vn datoseto to^ ^ in 

WiFS system’s radiometric performance, a g°nt m eve op m e ’ . , r ,nl n>n Li te rs; optical instruments, and 

r=.a^ SSJ. - -alyl are described 

in Chapter 6. 

This volume. Revision 1, supersedes the earlier version < M " a “ £ 

SeaWiFS Technical Report Sene s. Instrument P e f a „ d thete wave length changes for 
previous version (e,g„ radiometric saturation levels have been to dlscrimi „ate divinyl and 

some channels). Pigment protocols have been strengthened (R. B d«am^r,nm y ^ ^ reorgimi „ d | C . 
monovinyl chlorophyll a. Spectrophotometnc a sorp ion pro , N protocols have been added for 

Trees, B. Mitchell, J. Cleveland, C. Roesler, C. Yentsch, and I D. ( 7 above . wa te, measurements 

situ measurements of spectral absorption and backscatter.ng l^ Z “' ve “l’ ^ c Davi8 ], and fo, 

of water-leaving radiance from ships and low-flying aircraft |R. Doe"". K L “ aer ’ “ 
samphng in cLl water [R Doerffer, K. Carder, C. Davis, and R. Arnonel- 

In general, the specifications and protocol, set forth here simply »■! -dap. 

and procedures that are common practice in the ocean optics com" challenging protocols should, 

call fo, significant improvements over today’s - adopting these 

at least for the present, be regarded more as g neressarv to meet the rigorous SeaWiFS uncertainty 

goals as protocols are that the lmprovemen s ca e closely approached with a significant but affordable 

s£ ~ 

^ ln M^sl!llity s^sTd IxpLnental verification are needed to 

2. laboratory research is needed to improve absolute 

absolute calibration procedures, to achreve or approach 1% internal consistency m P s 

calibrations of radiometers to be used in SeaWiFS v.lidatton exper.ments (Sectron The SeaW l S 
Intercalibration Round-Robin Experiment (SIRREX) program has made s.gmficant progress 
area (Mueller 1993 and Mueller et al. 1994). 

3 Radiometric linearity test procedures must be improved to extend linearity characterizations over the 

f 11 nnpratinff dvnamic ranges (Table 4) of the various irradiance and radiance sensors (Sect • • )• 

This is especially critical for downwelling irradiance measurements at the sea surface, where irra laI1 “ S 
of laboratory irradiance standards are only 2-15% (depending on wavelength) of saturation irradia - 

, instrument 

elrm can be corrected to within less than 5%. The maximum diameter which enables such a correc- 
tion varies with the absorption coefficient, and is therefore a function of wavelength, water mass, an 
the solar zenith angle. For oligotrophic to moderately turbid coastal water masses, wavelengths 


t Names in brackets are the principal contributors to the revised protocols. 
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than eOOnm, and solar zenith angles greater than 30°, upwelled irradiance and radiance data from 
many of the currently popular radiometers having diameters of 20-40 cm can be adequately corrected 
In more general conditions, however, new instruments must be developed to minimize self-shadine 

b^ naSlv ^ for near-infrared wavelengths and in Case-2 waters. Self-shading corrections have 
been partially verified in recent experiments and a provisional protocol is included in this revision 

5 ' the / ringent pr ° t0C0ls f0r avoidin S shi P sh ^ows and reflections will require 

Tethered f f profiIm S ra d>ometer configurations which are not widely used today (Section 5 1 1) 

More inh re ^ fa t T T t0 ^ the m ° St economical approach to meeting these requirements 

“ eX ?, enS1Ve a PP roaches include optical systems on either remotely operated 
vehicles (ROVs), or on small surface platforms with self-contained winches. 

6 ' Jserforte^WFrrS 12 ^ 1 ? 0 ^ P ,°\^ &tion s ™ si ^ty is critical for any airborne radiometer to be 
. SeaWlF ^ radiometric validation, or algorithm development and validation. Protocols and 
p edures for polarization sensitivity characterization must be developed in more specific detail than 
was accomplished here because of time constraints (Sections 3.3 and 4.3). 

7 ' tW, UnCertain ^ ies sP^ified here for cosine responses of irradiance collectors are significantly better 
than is typically realized in commercially available radiometers (Section 3.1.5). Moreover the spec 
ified uncertainties may challenge the precision of the laboratory procedures used to characterize an 
instrument’s cosine response (Section 4.1.5). Error in cosine response almost directly transkte to an 

Z V Ther e f rr0r m dow r mnS ff irradianCe f ° r the dear - 8kieS case so critical to the sLwTfI vahd a 
. herefore, a significant effort to carefully characterize the effect, and to work with instrument 

manufacturers to approach or achieve the specific uncertainties, is an important factor in our strategy 
for reducing the overall error budget for water-leaving radiance measurements to less than 5%. 

8 ' VtuA °\ 8 POrtab,e s ‘ aadard to trace a radiometer’s performance stability during the course of 
eld deployment is called for in the present protocols (Section 4.2.5). Several manufacturers offer 
reasonably portable radiometric sources, which may be suitable for this purpose but laboratory 
and field evaluations must be carried out to prove their suitability and to develop new and detailed 
procedures for their use in the field. The SeaWiFS Project OfTW fqpni P i , , 

National Institute of Standards and Technology (NIST) to develop a prototyp^^m. * ’ 

9 The present protocols for deploying and analyzing data from moored and frersdrifting optical systems 

u^succ^smy TZZT ‘"'° ,npl t te ' AI,hOUgh m ” red “ d *“‘“8 systems hav^n 

used successfully in several oceanographic experiments, there are no previous examples of their use 

the W „^,evM“ P * f “ "" king ,he * ““““‘s I* developed and distributed with 

STcriUci dl^der^ ^ ““ * tomThT.bove 

for ocean optical measuremeTT^llt appZSL TdeTell at 7 “fT" *° 

a,SOrithm ^liaatrr„ , ^S 


San Diego , California 
December 1994 


— J. L. Mueller 
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Abstract 

This report presents protocols for measuring optical properties, and other environmental variables, to validate the 
radiometric performance of the Sea-viewing Wide Field-of-view Sensor (SeaWiFS), and to develop and validate 
biooptical algorithms for use with SeaWiFS data. The protocols are intended to establish foundations for a 
measurement strategy to verify the challenging SeaWiFS uncertainty goals of 5% in water-leaving radiances and 
35% in chlorophyll a concentration. The protocols first specify the variables which must be measured, and briefly 
review the rationale for measuring each variable. Subsequent chapters cover detailed protocols for instrument 
performance specifications, characterizing and calibrating instruments, methods of making measurements in the 
field, and methods of data analysis. These protocols were developed at a workshop sponsored by the SeaWiFS 
Project Ofice (SPO) and held at the Naval Postgraduate School in Monterey, California (9-12 April 1991). This 
report began as the proceedings of the workshop, as interpreted and expanded by the authors and reviewed by 
workshop participants and other members of the bio-optical research community. The protocols are an evolving 
prescription to allow the research community to approach the unprecedented measurement uncertainties implied 
by the SeaWiFS goals; research and development an; needed to improve the state-of-the-art in specific areas. 
These protocols should be periodically revised to reflect technical advances during the SeaWiFS Project cycle. 
The present edition (Revision 1) incorporates new protocols in several areas, including expanded protocol 
descriptions for Case-2 waters and other improvements, as contributed by several members of the SeaWiFS 
Science Team. 


Prologue 

The NIMBUS-7 Coastal Zone Color Scanner (CZCS) 
introduced ocean color remote sensing as a powerful new 
tool for observing ocean bio-optical properties. Through 
the early 1980s, CZCS data were exploited by a growing 
number of scientists studying marine phytoplankton, ocean 
productivity, and ocean optical properties. Practical ap- 
plications to marine fisheries were also demonstrated. Un- 
fortunately, a successor ocean color imaging system was 
not developed before the CZCS ceased operating in mid- 
1986. At present, therefore, research in these areas is lim- 
ited to retrospective, albeit productive, investigations of 
the CZCS historical database. In mid-1995, the launch of 
the Sea-viewing Wide Field-of-view Sensor (SeaWiFS), tin; 
next generation ocean color sensor, will bring a welcomed 
and improved renewal of ocean color time-series observa- 
tions to the ocean research community. 

The National Aeronautics and Space Administration 
(NASA) SeaWiFS Prelaunch Science Working Group 
(SPSWG) has recommended baseline satellite ocean color 
products consisting of: 

1) normalized water-leaving radiance (Lwtv) at five 
wavelengths; 

2) aerosol radiance at three wavelengths; 

3) chlorophyll a concentration; 

4) a chlorophyll-like pigment (chlorophyll a plus 
phaeopigment a) concentration; 

5) the diffuse attenuation coefficient ( K ) at 490 nm, 

AT (490); and 

6) calibrated radiances observed at the satellite. 


The primary SPSWG goals for product uncertainty are 
derived water-leaving radiances to within 5% and chloro- 
phyll a concentration to within 35% in Case-1 waters, both 
globally and throughout the projected five-year mission. 
These goals have? been endorsed by the SeaWiFS Science 
Team and accepted by NASA. NASA has the responsibility 
to lead the product assurance, calibration, and validation 
program. NASA is also required to determine the degree to 
which the commercially procured ocean color data fulfills 
the contractually stated NASA requirements. 

The CZCS mission was unquestionably a scientific suc- 
cess, but it also taught the participants that the satisfac- 
tory performance of a remote sensing satellite system can- 
not be taken for granted. Immediately after launch, and 
periodically throughout its five-year mission, the SeaWiFS 
system performance- including algorithms must be in- 
dependently verified using in situ optical measurements of 
the ocean and atmosphere. It is imperative that these sup- 
porting optical measurements meet a uniform standard of 
quality and accuracy if the primary SeaWiFS goals of 5%. 
uncertainty in water-leaving radiance and 35% uncertainty 
in chlorophyll a concentration are to be met, or even closely 
approached. To that end, the NASA Goddard Space Flight 
Center (GSFC) SeaWiFS Project convened a workshop to 
draft protocols and define standards for optical measure- 
ments to be used in SeaWiFS radiometric validation, and 
algorithm development and validation. The original ver- 
sion of this document (Mueller and Austin 1992) reported 
the protocols agreed to by the participants, as expanded 
by the authors, made in consultation with both the work- 
shop participants and other members of the ocean bio- 
optics community. The present revision incorporates the 
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expanded protocols and revisions contributed by several 
members of the SeaWiFS Science Team and other col- 
leagues. 

This report specifies the type and quality of support- 
ing in situ optical measurements and analytical protocols 
that will be required to develop bio-optical algorithms and 


validate the SeaWiFS calibration. Economics dictate that 
these observations and data will accrue over several years 
from a variety of sources, using different instruments and 
approaches. These data must be internally consistent, of 
known and documented uncertainty, and in a form readily 
accessible for analysis by the SeaWiFS investigators. 
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Chapter 1 


Overview of SeaWiFS Bio-Optical Requirements 

Introduction 

The Ocean Optics Protocols for SeaWiFS Validation , Revision 1 are intended to provide standards, which if 
followed carefully and documented appropriately, will assure that any particular set of optical measurements 
will be acceptable for SeaWiFS validation and algorithm development. It is true that in the case of ship shadow 
avoidance, for example, there are some circumstances in which acceptable radiometric profiles may be acquired 
considerably closer to a ship than is specified here (Section 5.1.1). When the protocols are not followed in such 
cases, however, it is incumbent upon the investigator to explicitly demonstrate that the actual error levels are 
within tolerance. The most straightforward way for an investigator to establish a measurement that is both 
accurate enough to meet the SeaWiFS standards and uncontaminated by artifacts, such as ship shadow, will be 
to adhere closely to the protocols. 


1.1 OBJECTIVES 

Immediate concerns have focused this document on spe- 
cific preparations for the SeaWiFS mission. A longer-term 
intention is the development of bio-optical databases that 
are relevant to future needs, and therefore, this document 
also recognizes the capabilities of other planned or poten- 
tial ocean color sensors, including: 

a) the Japanese Ocean Color Temperature Sensor 
(OCTS); 

b) the Moderate Resolution Imaging Spectroradi- 
ometer (MODIS); 

c) the Medium Resolution Imaging Spectrometer 
(MERIS); and 

d) the German Reflecting Optics System Imaging 
Spectrometer (ROSIS). 

The key objective of the working group was to rec- 
ommend protocols and standards for supporting in situ 
optical measurements. These objectives addressed the fol- 
lowing subject areas: 

1. The required and useful optical parameters to be 
used for validation of SeaWiFS normalized water- 
leaving radiances and atmospheric correction algo- 
rithms, and for monitoring the satellite sensor’s cal- 
ibration and stability, will be defined. 

2. The instrumentation requirements, and standards 
for measuring the parameters in item 1, includ- 
ing definitions of measured quantities, wavelengths, 
field-of-view (FOV) and band specifications, sensi- 
tivity, uncertainty and stability, will be delineated. 

3. The optical instrument characterization, intercali- 
bration standards, and related protocols will be de- 


fined. This objective includes the following sub- 
jects: 

a) laboratory calibration and characterization 
measurements, uncertainties, and procedures 
to be applied to instruments used in Sea- 
WiFS validation and algorithm development 
activities; 

b) pre- and post-deployment measurements and 
procedures to be followed with moored in- 
strumentation; and 

c) methods for instrument calibration and char- 
acterization, and the requirements for record 
keeping and traceability, including intercali- 
brations of radiometric and optical standards 
between participating laboratories. 

4. The at-sea optical sampling strategy and protocols 
will be standardized. This objective includes such 
considerations as: 

a) the rationale and justifications for moored, 
underway, drifting, shipboard, and airborne 
measurements; 

b) ship shadow avoidance, depth resolution in 
optical profiles, and total sampling depths; 
and 

c) time of day, sky conditions, season, and geo- 
graphic considerations. 

5. The analysis approaches to be used shall be refined. 
This objective includes procedures and methodolo- 
gies recommended for generating variables from in 
situ observations, e.g., Lwn(z) from L u (z), K(z), 
remote sensing reflectance, etc., as well as error 
analysis. 
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6. Protocols for ancillary measurements, data archiv- 
ing, database population, and access to data will be 
standardized. 

7. The required atmospheric measurements will be de- 
fined, and the degree to which standard methodolo- 
gies are available will be evaluated. 

The development and validation of bio-optical algo- 
rithms for SeaWiFS will be addressed by a separate work- 
ing group, and thus, these topics are briefly examined in 
this report. Nonetheless, the SPSWG was charged with 
identifying data requirements arid sampling strategies for 
bio-optical support measurements in the context of the op- 
tical and radiometric measurements. This topic includes 
the following subjects: 

1. Discrete chlorophyll a and pigment concentra- 
tions will be measured using the US Joint Global 
Ocean Flux Study (JGOFS) program’s proto- 
cols and standards for high performance liquid 
chromatography (HPLC) pigment sampling and 
analysis, which are adopted by reference to the 
JGOFS Core Measurement Protocols (JGOFS 
1991), Chapter 9, “Pigments and Chlorophyll.' ” 

2. An assessment will be made of the roles of un- 
derway, moored, and discrete chlorophyll a fluo- 
rescence measurements, how such measurements 
are calibrated, and their usefulness for satellite 
data product validation. 

3. The need* for biogeochemical measurements of 
colored dissolved organic material (CDOM), coc- 
coliths, suspended sediment, detritus, etc., will 
be examined on the basis of baseline product 
requirements. This subject will also the ad- 
dress the extent to which standards and proto- 
cols have been defined for such biogeochemical 
measurements. 

This report is complementary to anticipated reports of US 
and International JGOFS working groups, which are con- 
currently evaluating bio-optical needs and sampling strate- 
gies for their respective science programs. 

1.2 SENSOR CALIBRATION 

The SPO must make every effort to track the sensor’s 
performance throughout the duration of the mission. Since 
the instrument will be designed for a five-year mission, it 
is certain that the sensor calibration at each wavelength 
will change in some unpredictable manner as a function 
of time. Experience with CZCS has shown it is very dif- 
ficult to determine a sensor’s calibration once it has been 
launched (Viollier 1982, Gordon et al. 1983, Hovis et al. 
1985, Mueller 1985, and Gordon 1987). Similar prob- 
lems have been encountered with other Earth observing 
systems, such as the National Oceanic and Atmospheric 
Administration (NO A A) Advanced Very High Resolution 


Radiometer (AVHRR) (Brown and Evans 1985 and Wein- 
reb et al. 1990). Because of the large atmospheric contri- 
bution to the total observed radiances (Gordon 1981) and 
the great sensitivity of the bio-optical algorithms to the 
estimated water-leaving radiances (Clark 1981), small er- 
rors in the calibration can induce sizable errors in derived 
geophysical products, rendering them useless for many ap- 
plications. 

By processing large quantities of so-called “clear w'a- 
ter” imagery, i.e., water with pigment concentrations less 
than 0.25 mg m -3 (Gordon and Clark 1981), Evans and 
Gordon (1994) were able to develop a vicarious calibration 
that was used in the global processing of the entire CZCS 
data set (Esaias et al. 1986 and Feldman et al. 1989). The 
approach, however, requires assumptions that may limit 
the scientific utility of ocean color imagery. Specifically, 
the normalized clear w r ater- leaving radiances, L\vn (443), 
L\vn (520). and Lwn (550), were assumed to be 1.40, 0.48, 
and 0.30 mW cm -2 /an - 1 sr -1 , respectively. The Ang- 
strom exponents were assumed to be zero and certain geo- 
graphical regions such as the Sargasso Sea were assumed to 
be clear water sites at all times. Under these assumptions, 
analyses of the derived (Lwn) values were used to calculate 
calibration adjustment coefficients to bring CZCS derived 
(Lwn) values into agreement for these regions. The vi- 
carious calibration of the 443 nm band is tenuous because 
of the great variability in Lwn (443) even in clear water. 
Additionally, certain command and engineering data from 
the NIMBUS- 7 platform were not archived, so that a de- 
tailed analysis of possible effects related to the spacecraft 
environment and the effects of spacecraft operation on the 
calibration could not be performed. 

Unlike CZCS, SeaWiFS will routinely produce geophys- 
ical fields in a near-real time, operational mode for distri- 
bution to the science community. This aspect of the mis- 
sion necessitates constant evaluation of the sensor perfor- 
mance and the derived products. Therefore, a multifaceted 
approach to address the problem of sensitivity degradation 
and sensor characterization is required during both the pre- 
and post- launch phases. The goal is to ensure that Sea- 
WiFS level -1 radiances are accurately known and meet the 
specifications of the SPSWG. 

The plan includes both onboard and vicarious calibra- 
tion approaches. SeaWiFS will have a solar measuring 
diffuser plate to reference the response to the sun (Gordon 
1981) and also will be capable of periodically imaging the 
moon by maneuvering the spacecraft. The vicarious cali- 
bration program will incorporate measurements of water- 
leaving radiances and other related quantities, from ships, 
drifting buoys, and fixed moorings, to develop time series 
and geographically diverse samples of oceanic and atmos- 
pheric data sets. Each approach has advantages and disad- 
vantages, but when combined, they should provide a com- 
plementary and comprehensive data set that will be suffi- 
cient to monitor short-term changes and long-term trends 
in the sensor’s performance. 
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1.3 BIO-OPTICAL ALGORITHMS 

The SPO will be responsible for producing a standard 
set of derived products and will produce both CZCS-type 
products and baseline products. The CZCS-type products 
will consist of pigment concentration, K (490), five normal- 
ized water-leaving radiances, and three aerosol radiances 
based on constant default wavelength dependence (epsilon) 
coefficients in aerosol corrections. The proposed baseline 
products will include five normalized water- leaving radi- 
ances, K( 490), chlorophyll a concentration, three aerosol 
radiances, and one or more error analysis products. 

The basic quantities to be computed from the sensor 
radiances are the water-leaving radiances, from which all 
other derived products except the aerosol radiances are 
computed. Every effort must be made to ensure these ra- 
diances meet the specifications of the SPSWG, i.e., ±5% 
in Case-1 waters. This requires the atmospheric correction 
algorithms to be considerably more sophisticated than the 
current CZCS algorithms. 

The baseline bio-optical products must meet the ac- 
curacy requirements established by the SPSWG over a 
variety of water masses. The current CZCS algorithms 
were based on a data set consisting of fewer than 50 data 
points (only 14 observations were available for the band- 
2-to-band-3 ratio algorithm) and performed poorly in re- 
gions of high chlorophyll a concentration, high suspended 
sediment concentration, high CDOM concentration, and 
coccolithophorid blooms (Groom and Holligan 1987). Ac- 
curate estimates of the baseline products are essential if 
SeaWiFS is to be useful in programs such as the Global 
Ocean Flux program [National Academy of Science (NAS) 
1984]. 

SeaWiFS will have the capability, due to improvements 
in the signal-to- noise ratio (SNR), digitization, dynamic 
range, and wavelength selection, to increase the accuracy 
of these products and to flag areas where anomalies or 
low confidence conditions exist. Clearly, a much larger 
database will be needed for developing a broader variety 
of bio-optical algorithms, some of which will be region spe- 
cific. The radiometric, optical, and chemical field observa- 
tions used in deriving bio-optical algorithms and for vicar- 
ious calibration of the sensor must, therefore, conform to 
stringent requirements w T ith respect to instrument calibra- 
tion and characterization, and must also conform to the 
observation protocols which have been specified to take 
advantage of SeaWiFS capabilities. 

The SPO will manage a program to compare the vari- 
ous atmospheric correction and bio-optical algorithms pro- 
posed by the science community. The purpose of this pro- 
gram is to independently evaluate suggested improvements 
or additions to the SeaWiFS products. This component of 
the calibration and algorithm development program will 
run in parallel with, but off-line from, operational process- 
ing and will provide an essential mechanism for incorpo- 
rating data and analyses from the community at large. 


1.4 COMMUNITY PARTICIPATION 

The SPO will rely on the oceanographic community to 
perform field research for atmospheric and bio-optical al- 
gorithm development, and for all of the in situ data collec- 
tion for the vicarious sensor calibration. A minimal subset 
of these observations will be sponsored by the SPO, but 
many projects sponsored by the NASA Research and Ap- 
plication Program and other government agencies are ex- 
pected to make major contributions to the global five-year 
effort. This requires close coordination of the various pro- 
grams involved and a clear definition of the observations, 
accuracies, and data collection protocols required for each 
type of activity. The purpose of this document is to clarify 
these requirements. 

1.5 VICARIOUS CALIBRATION 

For ocean observations, it is easy to show (Gordon 
1987 and Gordon 1988) that satellite sensor calibration 
requirements based on the quality of the existing CZCS 
pigment algorithms exceed currently available capabilities. 
Furthermore, the sensor calibration is unlikely to remain 
unchanged through launch and five years of operation in 
orbit. The only foreseeable way of approaching the ocean 
calibration needs is through vicarious calibration, i.e., fine 
tuning the calibration in orbit. 

The methodology used to achieve vicarious calibration 
for CZCS was described in detail by Gordon (1987). First, 
the calibration was initialized after launch by forcing agree- 
ment between the sensor-determined radiance and the ex- 
pected radiance based on radiometric measurements made 
at the surface under clear atmospheric conditions. Next, 
since the CZCS responsivity was observed to be time de- 
pendent, the algorithms were applied to other scenes char- 
acterized by bio-optical surface measurements and more 
typical atmospheres, and the calibration was adjusted un- 
til the measured water-leaving radiances were reproduced. 
Finally, the surface measurements of pigments were com- 
bined with satellite pigment estimates for a wide variety of 
atmospheric conditions, and the radiance calibration was 
fine tuned until the best agreement was obtained between 
the retrieved and true pigments. 

The CZCS vicarious calibration was not radiometric. 
It was a calibration of the entire system— sensor plus al- 
gorithms. To predict the radiance measured at the sat- 
ellite, L t , the water- leaving radiance, the aerosol optical 
thickness, and the aerosol phase function are all required. 
Also needed are ancillary data such as the surface pres- 
sure, wind speed, and ozone optical thickness. These data 
for vicarious calibration and validation will be obtained 
by measuring the upwelling radiance distribution just be- 
neath the surface, along with the aerosol optical thickness 
and the sky radiance, at the time of the satellite overpass. 
The sky radiance will be used to deduce the required infor- 
mation about the aerosol phase function (Voss and Zibordi 
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1989). The data set will be used to deduce L t , at the top of 
the atmosphere, coincident with a SeaWiFS overpass from 
which the calibration will be initialized. 

It must be stressed that this exercise is absolutely es- 


sential for calibrating the SeaWiFS system, i.e., sensor plus 
algorithms, and that it cannot be effected without a high 
quality surface data set obtained simultaneously with the 
satellite imagery. 
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Chapter 2 


Data Requirements 

Introduction 

The prime objective of in-water optical measurements for SeaWiFS is to derive accurate normalized water- 
leaving radiances that will be used both for direct validation comparisons with those derived from Sea\\iFS 
data, and to develop and validate in-water bio-optical algorithms. Therefore, a comprehensive field program 
measure optical and biogeochemical state variables will be required. 


2.1 DEFINITIONS 

The required and useful variables to be measured for 
SeaWiFS validation are listed in Table 1, classified into 
two discrete sets. The first set comprises variables that 
can be used for radiometric initialization and ongoing val- 
idation. The second set encompasses those variables that 
will be used for bio-optical algorithm development and val- 
idation. 

Surface incident spectral irradiance, E d { 0“,A); down- 
welled spectral irradiance, E d {z, A); and upwelled spectral 
radiance, L u (z, A), are the fundamental measurable quan- 
tities needed to derive normalized water-leaving radiances 
in most circumstances. Other ambient properties, like sky 
radiance, sea state, wind velocity, etc., are also useful ini- 
tialization and calibration measurements and are discussed 
below. 

Surface incident spectral irradiance , E d { 0“, A), is usu- 
ally derived from surface irradiance, E s ( A), measured on 
a ship well above the water, but the use of a radiometer 
floated just beneath the surface (z = O') may provide a 
better approach (Sections 3.2.1 and 5.1.4). E d ( 0 , A) varies 
due to fluctuations in cloud cover and aerosols, and with 
time of day, i.e., solar zenith angle. Profiles of E d (z, A) 
and L u (z, A) must be normalized to account for such vari- 
abilities during a cast. 

Downwelled spectral irradiance, E d (z , A), is required to 
compute the diffuse attenuation coefficient, K(z , A), which 
in turn, is needed for diffuse attenuation coefficient algo- 
rithm development (Austin and Petzold 1981), and for 
optically weighting the pigment concentrations to be es- 
timated from remotely sensed ocean color (Gordon and 
Clark 1981). E d (z, A) is also required to compute the spec- 
tral (remote sensing) reflectance, Rl{z, A), which is used 
to normalize L u (z , A) when developing and validating bio- 
optical algorithms. The need for this normalization arises 
because the spectrum of incident irradiance varies with 
changing solar zenith angle and atmospheric conditions. 
E d (0', A) can then be used, through Rl{z, A), to convert 


L u { 0~, A) measured under a given set of illumination con- 
ditions, e.g., overcast skies, to Lvvn(A), which would be 
measured under the restricted illumination and viewing 
conditions associated with SeaWiFS measurements. As 
with L u (0", A), E d { 0' , A) must be determined by extrapo- 
lation from a profile of E d (z, A) over the upper few optical 
depths and reconciled with the direct surface measurement 
of E s ( A). [Optical depth, r(z,A), in the context of this 
report is the integral of K(z, A), for either radiance or ir- 
radiance, depending on the context, from the surface to a 
given depth z.] 

Upwelled spectral radiance, L u { O', A), is the in-water 
variable which, when propagated upward through the sea 
surface, leads to the measured value of L w ( A)* L\v( A) 
is, in turn, adjusted using E d { A) to derive the normalized 
water-leaving radiance, Lwn{ A), for a clear-skv zenith sun 
at the mean Earth-sun distance. Unfortunately, it is not 
practical to measure L u (0~ , A) precisely at an infinitesimal 
depth below the surface. Therefore, the profile of L u (z, A) 
must be measured over the upper few optical depths with 
sufficient accuracy to determine K L {z , A) for L u (z, A), and 
to propagate L u (z, A) to the surface. At near-infrared (IR) 
wavelengths, the first optical depth is confined to the upper 
few tens of centimeters. Determination of L u (0” , A) in this 
situation is more challenging and will require special in- 
struments and experiment designs to accommodate the ef- 
fects of instrument self-shading, wave focusing, small-scale 
variability, possible fluorescence, Raman scattering, and 
extremely small working volumes. Careful measurements 
of inherent optical properties (IOPs), including a(z,A), 
c(z,A), and b b (z, A), and spectral fluorescence, may be 
useful, in addition to E d (z,\) and L u (z,A) measurements 
made with specially designed radiometers. 

Sky radiance is required to enable estimation of the 
aerosol phase function through inversion of the radiative 
transfer equation. It is also useful for estimating the mean 
cosine of the transmitted light field in the water. The sky 
radiance should be measured directly; for the latter appli- 
cation, however, it need only be estimated by occulting the 
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TaMe 1 Required observations for initialization and system calibration for satellite 
transfer (also ongoing calibration and atmospheric algorithm validation studies) and 
and validation. ' 


product verification 
bio-optical algorithm 


and radiative 
development 



Product 

Verification 

Radiative 

Transfer 

Bio-optical 

Algorithms 

Primary Optical Measurements 


incident spectral Irradiance, E d (Q ,\) 

Downwelled Spectral Irradiance, E d (z, A) 

Upwelled Spectral Radiance, L u (z, A) 

Spectral Solar Atmospheric Transmission, t s ( A) 
Submerged Upwelled Radiance Distribution, L(z,6,<f>) 
Spectral Sky Radiance Distribution 
Upwelled Spectral Irradiance, E u (z, A) 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Calculated or Derived Variables " “ 


water- leaving Radiance, L w { 0“,A) 

Attenuation Coefficient Downwelled Irradiance, K d (z, A) 
Attenuation Coefficient Upwelled Radiance, K L (z, A)' 
Spectral Reflectance, Rl(z, A) 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Ambient Properties 


Dea ana bky State Photographs 
Wind Velocity 

In Situ Fluorescence Profiles 
Aerosol Samples 

Temperature and Salinity Profiles 
Secchi Depth 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Primary Biogeochemical Measurements 


ruytopianKton figments (HPLC Technique) 
Phytoplankton Pigments (Fluorometric Technique) 
Total Suspended Material (TSM) Concentration 
Colored Dissolved Organic Material (CDOM) 


X 

X 

X 

X 

X 

X 

X 

X 

inherent Optical Properties 

O 1 T . - 


spectral oeain Attenuation Coefficient, c(z, A) 
Spectral Absorption Coefficient, a(z, A) 
Spectral Backscattering Coefficient, b b (z, A) 
Spectral Volume Scattering Function, 0(z,\,6) 
Red Beam Attenuation, c(z, 660 nm) 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Algorithm Specific Kesearch Measurements 

A • 1 _ 


.mi uorne riuorescence and Radiances 
Coccolith Concentration 
Detritus Absorption Coefficient 
Humic and Fulvic Acids 
Inorganic Suspended Material 
Organic Suspended Material 
Particle Absorption Coefficient 
Particle Fluorescence 
Particle Size Spectra 
Particulate Organic Carbon (POC) 
Particulate Organic Nitrogen (PON) 
Phycobilipigments Concentration 
Phytoplankton Species Counts 
Primary Productivity ( 14 C) 

Total Dissolved Organic Carbon (DOC) 
x - Needed for the indicated effort. 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 
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sun’s image on a deck cell measuring the incident spectral 
radiance from the sun and sky. The mean cosine at the 
surface can be used with profile measurements of £<*( A), 
E u (\), and c(A) to estimate h{\) (Gordon 1991). An abil- 
ity to exploit this and similar relationships will greatly 
enhance both development and verification of bio-optical 
algorithms. The spectral sky radiance distribution over 
zenith and azimuth angles is required to determine the 
aerosol scattering phase functions at radiometric compar- 
ison stations during the system initialization cruises and 
will be very useful if measured at all validation stations 
throughout the mission. 

Upwelled radiance distribution measurements just be- 
neath the sea surface will be required for quantifying the 
angular distribution of water-leaving radiance at stations 
used for system calibration initialization and long-term 
system characterization. These measurements will also be 
useful in relating radiance and irradiance reflectance, and 
K profiles, to IOPs and biogeochemical substances, e.g., 
chlorophyll a and CDOM, during bio-optical algorithm de- 
velopment and validation. 

Atmospheric transmittance spectra should be measured 
using a sun photometer in order to determine aerosol op- 
tical depths at each station. These data are particularly 
needed to verify the atmospheric corrections in direct com- 
parisons between SeaWiFS Lw(A) estimates and those de- 
termined from in-water L u ( 0",A). 

Sea state photographs are required to document surface 
wave conditions during radiometric measurements. This 
information is essential for identifying measurements made 
under questionable environmental conditions. 

Wind velocity is required to generate, through mod- 
els, estimates of the surface wave slope distribution, which 
will be used to calculate reflected skylight and sun glint 
in radiative transfer models (Cox and Munk 1954). Sur- 
face wave models driven by wind velocity may also be used 
to provide quantitative estimates of surface wave induced 
radiometric fluctuations. Qualitatively, wind velocity, and 
photographs or videotape recordings of sea state, will be 
useful for assessing station data quality. 

Upwelled spectral irradiance, E u (z, A), is a useful mea- 
surement, in addition to Ed and L u , because there ex- 
ist both empirical and theoretical relationships between 
IOPs, phytoplankton pigments, TSM, and irradiance re- 
flectance. L u (0“,A) and £ U (0“,A) are related by the fac- 
tor Q( A), which is not well determined at present, and has 
been shown to vary with solar zenith angle (Morel and 
Gentili 1993). Combined measurements of L u (0~,A) and 
£ u (0', A) will be extremely useful in determining Q{ A), 
which will, in turn, allow traceability of SeaWiFS mea- 
surements to previously derived irradiance reflectance re- 
lationships and algorithms. 

IOPs must be measured for development and valida- 
tion of the SeaWiFS semi-analytic Case-2 chlorophyll a 
algorithm. This algorithm is based on an explicit theo- 
retical function of the ratio of backscattering to absorp- 
tion, b b (X)/a{ A). In the first implementation, however, 


6 6 (A)/a(A) is modeled statistically as a polynomial func- 
tion of chlorophyll a concentration. This will introduce a 
strong, empirical component to the algorithm, suppress- 
ing the physical links between Ri( A) and b b (X)/a(X) and 
between b b ( A)/a(A) and phytoplankton pigment concentra- 
tion. A more robust algorithm would be based on direct 
measurements of absorption, scattering, pigments, and re- 
flectance. Due to recent advances in instrumentation, it 
is now practical to routinely measure profiles of a(z, A) 
(Section 5.2.4) and backscattering variables (Section 5.2.5) 
from which h b {z, A) may be approximated. Future develop- 
ment and validation experiments involving this algorithm 
must, therefore, include absorption, beam attenuation, and 
scattering measurements. 

Red beam attenuation coefficient, c(660), and in situ 
chlorophyll a fluorescence measurements are exceptionally 
useful in analyzing profiles of Ed(z), L u (z), and E u (z) to 
derive profiles of Kd(z), Kl(z), and K u (z), respectively. 
If these profiles are viewed in real time, they are also use- 
ful guides for taking water samples at depths that allow 
the vertical structure of chlorophyll a and suspended par- 
ticles to be accurately resolved in the top optical depth. 
Finally, the chlorophyll a fluorescence profile is used to in- 
terpolate HPLC and extracted fluorescence measurement 
of chlorophyll a concentrations from water samples at dis- 
crete depths. It is desirable to make these measurements 
simultaneously with irradiance and radiance profiles, if it 
can be done in a way to avoid self-shading of the instru- 
ment (Section 5.1.6). 

Secchi depth measurements are required for real-time 
assessment of water transparency during a station and as 
a quality check during analysis of radiometric profiles. 

Aerosol concentration samples using high volume tech- 
niques will be useful, in conjunction with aerosol optical 
depth spectra determined from sun photometer measure- 
ments, for chemical, size, and absorption characterization 
of aerosols, especially in studies of the effects of Saharan 
and Asian dust clouds on atmospheric corrections. 

2*2 BIOGEOCHEMICAL DATA 

Pigment concentrations will be determined using HPLC 
and fluorometric methods to develop and validate chlo- 
rophyll a algorithms, and to assess the effects of acces- 
sory pigment concentrations on water-leaving spectral ra- 
diances. These data will also be used to calibrate con- 
tinuous profiles of chlorophyll a fluorescence (Section 2.1). 
Phytoplankton pigment concentration will also be deter- 
mined using classical chlorophyll a and phaeopigment flu- 
orescence techniques that were used for CZCS pigment 
validation and algorithm development. The HPLC tech- 
nique provides more accurate and precise information for a 
greater number of pigments, but gives different values than 
does the fluorescence method for various species composi- 
tions and chlorophyll a-to-phaeopigment ratios. While the 
HPLC method is the primary pigment technique required 
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for SeaWiFS, the classical technique is still required to al- 
low the CZCS and SeaWiFS data sets and algorithms to 
be compared. 

Phycobilipigments , present in cyanobacteria and cryp- 
tophytes, are treated separately from the HPLC fat sol- 
uble pigments. Phycoerythrin and phycocyanin are the 
two major groups of phycobilipigments found in the ma- 
rine environment. The concentration of these water solu- 
ble pigments is important due to the contribution of solar 
stimulated phycoerythrin fluorescence to the underwater 
light field, and also to characterize the phytoplankton pop- 
ulation. At times, species which contain phycobilipigment 
can account for a large fraction of the primary productivity 
(especially in oligotrophic waters) and have been difficult 
to quantify due to their small size. These measurements 
are not required because SeaWiFS does not contain bands 
at their absorption or fluorescence peaks. The measure- 
ments are desirable, however, since several aircraft sensors, 
e.g., the Airborne Visible and Infrared Imaging Spectrom- 
eter (AVIRIS), Airborne Oceanographic Lidar (AOL), the 
Multispectral Airborne Radiometer System (MARS), and 
future satellite sensors, e g., MODIS and MERIS, do pos- 
sess bands at the absorption or fluorescence peaks of phy- 
cobilipigments. If such pigment information is required for 
a particular study, it is recommended that the glycerol- 
uncoupling technique described by Wyman (1992) be used 
to determine phycoerythrin concentrations. 

TSM measurements are required to assess the effect of 
suspended sediment on the derived products. TSM is of 
primary importance in coastal waters, where simple radi- 
ance ratio algorithms for TSM have uncertainties equiva- 
lent to, or better than, those for estimating chlorophyll- like 
pigment concentration. Organic suspended matter and in- 
organic suspended matter concentrations are subfractions 
of TSM; this partitioning of TSM is particularly useful in 
process studies. 

Particulates, both POC and PON, are required for pro- 
cess studies to help characterize the adaptive state of phy- 
toplankton and to inventory critical biogeochemical ele- 
ments. 

DOC has been shown to be a major pool of carbon 
in the oceans. Quantification of the transformations of 
this pool is crucial to understanding the marine carbon 
cycle. The colored fraction, CDOM, of the DOC is highly 
absorbent in the blue range, thus decreasing blue water- 
leaving radiances, and it must be taken into consideration 
for pigment concentration algorithms. DOC measurements 
are needed to develop robust relationships between CDOM 
and DOC, which are needed to evaluate the usefulness of 
ocean color observations for estimating DOC concentra- 
tions. 

CDOM concentrations are required to assess the ef- 
fect of Gelbstoff on blue water-leaving radiances and chlo- 
rophyll concentration. This is of primary importance in 
Case-2 waters, but is also relevant to phytoplankton degra- 
dation products. 


Humic and fulvic acids comprise the bulk of CDOM 
and have different specific spectral absorption coefficients. 
Their concentrations are useful for determining the correc- 
tion used for phytoplankton pigment concentration algo- 
rithms in Case-2 waters and for estimating CDOM from 
ocean color observations. 

Coccoiith concentration , which is the number density 
of small plates (coccoliths) composed of calcium carbon- 
ate (CaC0 3 ), is very important to light scattering. Coc- 
coliths are produced in copious amounts by marine phy- 
toplankton called coccolithophores. Scattering of light by 
coccoliths is highly apparent in visible wavelength satel- 
lite imagery, because they perturb the usual relationships 
between water-leaving radiances and chlorophyll a concen- 
tration and adversely impact atmospheric corrections. Ad- 
ditionally, coccoiith formation, sinking, and dissolution are 
significant factors in the ocean carbon flux budget. It 
is, therefore, necessary to measure coccoiith concentra- 
tion, both as number density and CaC0 3 concentration, 
to aid in 1) the correction of chlorophyll a concentration 
algorithms. 2) coccoiith algoiithm development, and 3) at- 
mospheric correction development and validation. 

The particle absorption coefficient , which is comprised 
of absorption by living, dead, and inorganic particles, is 
a useful variable for modeling the portion of solar energy 
that is absorbed by phytoplankton and bacteria. 

Detritus (or tripton) absorption coefficient , i.e., absorp- 
tion of light by detritus, represents a major loss of light 
which would otherwise be available to the phytoplankton 
component of the marine hydrosol. In many cases, absorp- 
tion by detritus is a significant term in the marine radiative 
transfer processes, and its determination is useful for phy- 
toplankton production models and for modeling the light 
field. 

Particle size spectra are very useful for in-water ra- 
diative transfer calculations, particularly if measurements 
include particles smaller than 1 pm. 

Particle fluorescence, measured using laser sources on 
single-cell flow systems, may be used to calculate particle 
scattering-to-fluorescence ratios for evaluating the popula- 
tion structure of the plankton (both phyto- and zooplank- 
ton). 

Phytoplankton species counts are important because 
species-to-species variability in optical and physiological 
properties represents a major source of variability in bio- 
optical algorithms and primary productivity models. This 
has been recognized, but it is generally ignored in remote 
sensing algorithms due to the tedious nature of species 
enumeration, the small sizes of many species, and the large 
number of species involved. This information, however, at 
various levels of rigor, is useful in evaluating the population 
and pigment composition. This is especially important for 
some groups, such as coccolithophores. 

Primary productivity, using the radioactive isotope 14 C 
estimation method, is not strictly required for validation of 
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water-leaving radiances or system initialization. Further- 
more, primary productivity is not a standard derived Sea- 
WiFS product, owing to the complexity of relating ocean 
color to production. It will, however, be extremely useful 
for process study applications of ocean color data if these 
measurements are made at the same time that the water 
column optical properties are determined. These data will 
aid in the development of models of primary production 
using satellite ocean color observations, a goal which is cen- 
tral to the overall SeaWiFS science mission. Of special im- 
portance are determinations of key photo-physiological pa- 
rameters derived from production measurements as func- 
tions of irradiance. If 14 C productivity measurements are 
made, they should conform to the JGOFS Core Measure- 
ments Protocols (JGOFS 1991). 

2.3 ABOVE- WATER TECHNIQUES 

Above-water ocean color radiance measurements from 
aircraft and ships can augment in-water measurements of 
L u (z, A) made to compare directly with SeaWiFS measure- 
ments for validation of its radiometric performance and 
algorithms. Above- water radiance measurements can, if 
they are accurately made, contribute an additional use- 
ful constraint in defining internally consistent sun-ocean- 
atmosphere-sensor models that will comprise the essence 
of SeaWiFS radiometric validation. For this application, 
above-water radiometers must meet the SeaWiFS specifi- 
cations for radiometric uncertainty, SNR, and spectral re- 
solution at a spatial resolution that will permit direct com- 
parisons with in-water measurements, and with SeaWiFS 
measurements by means of spatial averaging. Conversely, 
above-water radiance measurements made with less accu- 
racy than the SeaWiFS prelaunch specifications would in- 
troduce an unacceptable error source into the validation 
models and cannot be used for this purpose. 

Airborne ocean color data may also be used to deter- 
mine spatial variability in ocean optical properties during 
shipboard algorithm development and validation experi- 
ments. In this context, airborne ocean color measurements 
will be especially valuable in productive Case-1 and Case- 
2 waters, where variability in ocean optical properties can 
be large over mesoscale and smaller distances. Synoptic 
maps of ocean color distributions can be advantageously 
utilized to guide sampling by ships. It can also be used 
to place in-water data from an individual station in con- 
text with respect to nearby variability, and thus provide 
a basis for spatial interpolation and averaging when com- 
paring in-water bio-optical measurements with SeaWiFS 
image data. This application can be accomplished using 
aircraft radiometers meeting somewhat less stringent per- 
formance specifications than is demanded for direct vali- 
dation comparison between SeaWiFS and aircraft radiance 
measurements. 

Airborne measurements of fluorescence by chlorophyll, 
CDOM, and phycoerythrin, both by laser and solar excita- 
tion, are useful to evaluate spatial and temporal variability 


near ship and mooring stations and to provide independent 
assessments of bio-optical algorithms. 

2.4 ANCILLARY MEASUREMENTS 

Hydrographic data, water temperature (T), and salin- 
ity (S), derived from conductivity, temperature, and depth 
(CTD) profiles, are useful for characterizing the physical 
water mass regime in which an optical profile is measured. 
A T-S characterization is especially important near ocean 
fronts and eddies where interleaving water masses of very 
different biogeochemical composition, and therefore funda- 
mentally different bio-optical properties, can produce com- 
plex spatial and temporal patterns of near-surface optical 
properties. In these circumstances, T-S profiles can pro- 
vide an indication of whether a station location is suitable 
for reliable remote sensing validation and algorithm devel- 
opment comparisons. 

2.5 OPTICAL MOORINGS 

Optical moorings will be maintained in one or more re- 
gions of low optical variability to provide long-term time 
series comparisons between in situ and SeaWiFS measure- 
ments of normalized water-leaving radiance. Moored op- 
tical systems will also be extremely useful in a variety of 
oceanographic studies. For example, global satellite obser- 
vations of ocean pigment biomass and estimates of phyto- 
plankton production are essential to achieve the objectives 
of the JGOFS program (NAS 1984); SeaWiFS will play a 
key role in this effort. 

The oceans exhibit physical and biological variability 
over a wide range of space and time scales. This variabil- 
ity, and the need to synoptically measure distributions of 
physical and biological properties over large areas and long 
time periods, has motivated recent developments utilizing 
contemporaneous buoy, ship, aircraft, and satellite sam- 
pling strategies (Smith et al. 1987). In addition, lung-term 
mooring data are required to provide continuous observa- 
tions and permit an optimization of the accuracy of the 
derived satellite products (Booth and Smith 1988). 

There are two sources of systematic error in estimates 
of phytoplankton pigment biomass derived from satellite 
ocean color data. First, errors in satellite estimation of 
pigment biomass arise because physical forcing, biologi- 
cal properties, and ocean optical properties all vary sys- 
tematically with depth, and the upper layer optical signal 
observed by satellites may not adequately represent struc- 
ture deeper in the water column. In many circumstances, 
subsurface changes may go undetected unless contempo- 
raneous water column profile data from either shipboard 
or moored sensors is available. Second, visible wavelength 
sensor systems do not obtain data when the atmosphere 
is cloudy. Air-sea interactions giving rise to cloudiness are 
often closely linked to biological processes. For example, 
Michaelsen et al. (1988) showed that episodic wind events, 
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which give rise to coastal upwelling and subsequent phy- 
toplankton production along the California coast, cause 
cloudiness that will bias the statistics of pigment concen- 
trations derived from ocean color imagery. Similar biases 
of a factor of 3-4, due to wind mixing during cloudy pe- 
riods, were observed by Muller-Karger et al. (1990). In 
circumstances of this nature, visible and infrared satellite 
observations of the ocean are not random samples. Moored 
optical sensors can measure systematic temporal variabil- 
ity in the vertical distribution of pigment biomass and, 
at the same time, provide the continuous time series that 
may be used to remove the sampling bias associated with 
cloudiness. 

The detection and verification of intra- and interan- 
nual fluctuations in productivity and associated bio-optical 
variables are key goals of programs focused on studying 
global change. These goals place stringent requirements 
for long-term accuracy and precision on the measurement 
systems to be employed. The monitoring of bio-optical pa- 
rameters to resolve variability at global and decadal scales, 
as proposed by the SeaWiFS and MODIS missions, will 
require that moored in situ optical instruments be main- 
tained to supplement and support the satellite data sets. 
For example, the CZCS sensor degradation and the difficul- 
ties encountered in attempting to characterize that degra- 
dation, strongly point out how valuable in- water optical 
measurements would have been to that program. 

2.6 DRIFTING OPTICAL BUOYS 

Drifting optical buoys , which are expendable and anal- 
ogous to the ARGOS sea surface temperature (SST) drif- 
ters, represent a viable, cost-effective way to obtain sig- 
nificant numbers of daily optical observations to validate 
global data sets of water-leaving radiance. Due to high 
cost, there will probably never be more than a few per- 
manent optical moorings. Shipboard observations provide 
only a few measurements on any given day, due to both the 
cost of at-sea research and investigator availability. Op- 
tical data from drifting buoys, while less complete than 
measurements from ships and fixed optical moorings, can 
potentially surpass both in terms of global coverage and 
the number of near- real time comparisons with SeaWiFS 
observations. Judicious seeding of inexpensive drifters pro- 
vides one means of sampling conditions and regions critical 
to SeaWiFS product verification. For example, chronically 
high levels of aerosols caused by desert dust storms con- 
stitute a regional condition that will affect atmospheric 
corrections. These instruments may also allow study of 
possible variations in SeaWiFS performance as a function 
of latitude, due to orbital variations in sensor performance 
and sun angle dependency of algorithms. 

Optical drifter development activities occuring at Dal- 
housie University intend to develop instruments to measure 
seven and three upwelling spectral radiances, respectively, 
and a single downwelling irradiance. Along with temper- 
ature and barometric pressure, these data are transmitted 


over an ARGOS link. Storage procedures are designed 
to make maximum use of the limited ARGOS bandwidth. 
(An interrogating store and forward satellite system with 
greater bandwidth would be beneficial for this purpose.) 
The upwelling radiance data are obtained at a depth of 
approximately 0.5 m and must be propagated through the 
surface using Ki{ A) estimated from the relative spectral 
shape. These water-leaving radiance estimates will, there- 
fore, be inherently less accurate than surface reflectance 
observations made together with optical profiles and more 
complete ancillary observations. 

Both of these systems are in the test and evaluation 
phases. High risk areas which are being examined in- 
clude long-term stability, identification of bio-fouling ef- 
fects, operating lifetime, and validation of the techniques 
used for calculating the water-leaving radiance from the 
simple drifting sensors. The accuracy of the ARGOS sys- 
tem for drifter location is sufficient for global area coverage 
(GAC), but experience needs to be gained in analysis of 
the data to demonstrate the feasibility of using ARGOS 
positioning (150 m to 1km uncertainty) for system cali- 
bration and validation work in water mass regimes where 
mesoscale variability is significant. 

The cost of optical drifters will limit the number de- 
ployed. Proponents envision 50 such buoys adrift at any 
one time throughout the world — a number sufficient to pro- 
vide a large enough sample size to support viable global 
validation. In such a scenario, typically 60% of the drif- 
ters would be obscured by clouds during a SeaWiFS pass 
(those in some areas will have a much greater probability 
for clouds). Furthermore, current divergence areas will be 
systematically undersampled by drifters. 

The uncertainty of calculated Lw derived from drifter 
data has been estimated to be of the order of 15%, even 
though at the time they are deployed, the calibrated un- 
certainty of the instruments will be less than 5%. This 
estimate may be pessimistic, based on: 

1) the untest able possibility of drifts in radiomet- 
ric responses during long-term deployment of an 
expendable instrument, and 

2) errors associated with propagating L u (0.5 m, A) 
through the water column and interface to esti- 
mate LwW without benefit of measurements of 
K(z, A), surface roughness, and other ancillary 
measurements to be carried out at correspond- 
ing ship stations. 

If Lwn from drifters is only good to 15%, then they can- 
not be used to verify SeaWiFS radiometry within 5%, no 
matter how many drifter comparisons are made. Uncer- 
tainties must be less than 5% if this technology is to be 
used for SeaWiFS radiometric validation. Uncertainties 
of 15% may, however, be useful for validating SeaWiFS 
and derived products, and for interpolating SeaWiFS data 
through periods of extensive cloud cover. 
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Chapter 3 


Specifications 

Introduction 

This report describes measurements of optical properties, and other variables, necessary for validating data 
obtained with the SeaWiFS instrument, and for the development of in-water and atmospheric algorithms. The 
specifications herein are those required of instruments used on ships, or other platforms, to acquire that optical 
data. In some cases, the specifications have been selected to allow use of instruments that are affordable 
and that either currently exist, or that can be developed without major improvements in today’s state-of- 
the-art technology. In a few cases, new or improved instruments must be developed to realize the specified 
performance characteristics. The data uncertainty requirements for this program are more severe than those 
for a general ocean survey. Here, various investigators will use a variety of instruments that will be calibrated 
independently at a number of facilities, and contribute data to a common database which will be used to validate 
SeaWiFS measurements. The resulting radiometric and bio-optical database will provide an essential means 
of detecting and quantifying on-orbit changes in the SeaWiFS instrument relative to its prelaunch calibration 
and characterization. This chapter specifies instrument characteristics and data uncertainties thought by the 
SPSWG to be necessary, as well as sufficient, for this task. The validation analysis would be significantly 
degraded should calibration errors or differences of even a few percent, or wavelength errors or differences of a 
few nanometers, occur in (between) the instruments used to acquire the SeaWiFS bio-optical database. 


3.1 IN- WATER RADIOMETERS 

This section specifies radiometric characteristics for in- 
struments that are used to measure Ed(z, A), E u (z, A), and 
L u (z, A). The specifications are applicable to filter ra- 
diometers and to spectroradiometers based on monochro- 
mators. Minimum performance characteristics are spec- 
ified for spectral resolution, radiometric responsivity and 
resolution, SNRs, radiometric saturation and minimum de- 
tectable values, angular response, temporal sampling reso- 
lution, linearity, and stability. 

3.1.1 Spectral Characteristics 

In-water radiometers shall be capable, as a minimum, 
of making measurements at the wavelengths shown in Ta- 
ble 2. The table presumes the use of properly blocked in- 
terference filters to provide the required spectral bandpass 
and out-of-band rejection (10~ 6 or better). Care must also 
be taken to avoid possible out-of-band leakage due to flu- 
orescence by filter, or other optical component, materials. 
Filter radiometers should have channels with center wave- 
lengths, as measured in the assembled instrument, match- 
ing those given in Table 2 to within ±1 nm for 410 and 
443 nm, and within ±2 nm for all other spectral bands. 
Shifts of these magnitudes in center wavelengths will re- 
sult in changes in measured radiometric values of approx- 
imately ±1% or less (Booth pers. comm.) and this speci- 
fication should be met if possible. 


Table 2. Recommended spectral bands for discrete 
wavelength filter radiometers using 10 nm full- width 
at half-maximum (FWHM) bandwidths. In addi- 
tion, out-of-band blocking in the far tails of the in- 
strument response functions should be at least 10“ 6 . 


SeaWiFS 

Band 

Wavelengths 

[nm] 

F({ , E u , L u 
[nm] 

E, 

[nm] 

1 

402-422 

412m 

412 

2 

433-453 

443, 435 [2] 

443 

3 

480-500 

490 

490 

4 

500-520 

510 

510 

5 

545-565 

555 

555 

6 

660-680 

665, 683 

665 [3] 

7 

745-785 

a 

780 

8 

845-885 

a 

875(5] 


|T| A preferred option is to replace two separate 10 nm 
FWHM bands centered at 406 and 416 nm, with a 
single 412 nm channel. The two channels would al- 
low more accurate modeling of Lvw(412) matching 
SeaWiFS characteristics. 

\ 2 ] An optional extra band is used to improve modeling 
of Lwn(A) radiances to match the SeaWiFS 443 nm 
channel. 

f3l E s deck, only one channel in this band is necessary. 

[T| Due to the specialized nature of infrared in-water 
measurements, specialized sensors will be needed. 

O Optional for E s . 
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It is recognized, however, that enforcing ail nm hard- 
and-fast specification could be prohibitively expensive, and 
this tolerance should be regarded as a goal. With knowl- 
edge, to less than 0.2 nm, of the actual center wavelengths 
and complete spectral response functions, corrections prob- 
ably can be made to infer effective radiometric quantities 
for the SeaWiFS channels, when the spectral characteris- 
tics of SeaWiFS channels have also been measured, shortly 
before launch. Bandwidths must be 10 nm ±2 nm FWHM. 
They are made narrower than the SeaWiFS channels to re- 
duce the skewing of the parameters derived from underwa- 
ter irradiance or radiance profiles in spectral regions where 
absorption by natural sea water may exhibit rapid varia- 
tion with wavelength. 

To maintain the above tolerances, it is anticipated that 
filters will be ordered to a center wavelength with a toler- 
ance of Ao ± 1 nm and a FWHM bandwidth of 8.5 ± 1 nm. 
When the filter is installed in a radiometer with a 10° (half- 
angle) FOV, however, the spectral bandpass will broaden 
by 2-3 nm, and the center wavelength will shift. Further- 
more, as a filter ages in use, its transmission curve may 
undergo changes to further broaden the FWHM bandpass 
and shift the peak. The tolerances specified above include 
an allowance for some degradation before expensive filter 
and detector changes must be done. 


Table 3. High resolution spectroradiometric spec- 
ifications. 


Optical Sensors 

Spectral Range: 

380 to 750/900 nm 

Spectral Resolution: 

5 nm (or less FWHM) 

Wavelength Accuracy: 

10% FWHM of reso- 
lution (0.5 nm) 

Wavelength Stability: 

5% FWHM of reso- 
lution (0.25 nm) 

Signal- to- Noise Ratio: 

1,000:1 (at minimum) 

Stray Light Rejection: 

10" 6 

Radiometric Accuracy: 

3% 

Radiometric Stability: 

1% 

FOV Maximum: 

10° (for radiance) 

Temperature Stability: 

Specified for 0-35° C 

Linearity: 

Correctable to 0.1% 

Ancillary Sensors 

Temperature: 

0.2° C 

Pressure: 

0.1% (full scale) 

Horizontal Inclination 

1° over 40° range 


In a single instrument, all channels at a given nominal 
wavelength should match within 1 nm, if possible. It is 
desirable, therefore, to obtain all of the filters used by an 
investigator for measurements at any nominal wavelength 
(A n ) from a single manufacturing lot when possible. If this 
is done, E s ( A n ), Ed(X n ), E u ( A n ), L u { A n ), and any atmos- 
pheric radiometric quantities measured with that investi- 
gator’s systems, would all have a greater likelihood of be- 
ing measured over the same range of wavelengths, for each 


nominal wavelength (A n ). In any event, the actual spec- 
tral response function of each instrument channel must be 
measured and known to be accurate to less than 0.2 nm. 

High resolution monochromator- based spectroradiome- 
ters, with adequate sensitivity and stray light rejection 
characteristics, are also suitable instruments and are rec- 
ommended for many algorithm development studies. Suit- 
able specifications for such instruments are given in Ta- 
ble 3. (These instruments must also meet the specifications 
summarized in Tables 2 and 4.) 

3,1*2 Responsivity, SNR, and Resolution 

The expected operating limits for radiometric respon- 
sivities, SNR, and digital resolution are specified in Ta- 
ble 4, the limits for which were derived as follows: 

1. An Ed saturation value of 300/xWcm” 2 nm” 1 is 
assumed at all wavelengths. 

2. Implicit, but not stated, in Table 4 is that the min- 
imum required Ed{0) is 20/zWcm” 2 nm -1 ; it will 
not be appropriate to occupy validation stations 
when illumination is less than this minimum. 

3. The minimum Ed{ 0) implies a minimum detectable 
Ed(z) value of 1 p W cm” 2 nm” 1 at 3 optical depths 

( 3 /*). 

4. Digital resolution must be less than or equal to 
0.5% of the reading to maintain a 100:1 SNR. To 
permit a 1% uncertainty in absolute calibration, if 
that goal can be met in the calibration laboratory, 
the instrument must digitally resolve 0.1% of the 
irradiance (radiance) produced by the laboratory 
standards used; typical irradiance (radiance) values 
for calibration using 1,000 W FEL standard lamps 
traceable to the National Institute of Standards and 
Technology (NIST), and required digital resolutions 
at these signal levels, are given in Table 4 as “Cali- 
bration Irradiance” and “Digital Resolution (cal.),” 
respectively. A SNR of 100:1 requires a resolution 
in Ed{z) at three optical depths to 0.005 /iW cm” 2 
nm -1 per count, i.e., 2.5 digit resolution. At the 
surface, £^(0) should be resolved to 0.05 /xW cm” 2 
nm” 1 per count. 

5. The Case-1 saturation values of E u ( 0) represent 
the Instrument Specification Subgroup’s estimate of 
maximum reflectances to be expected in ordinary 
Case-1 waters: 12.5% at 410 nm, 7.5% at 488 nm 
and 0.5% at 670 nm. These saturation values will be 
too low for measurements in Case- 2 waters or cocco- 
lithophore blooms. In these situations, a maximum 
expected reflectance of 40% for A < 660 nm and 
20% for A > 660 nm is assumed. This implies that 
the expected maximum irradiance in £7 U (0) should 
be 80/xWcm“ 2 nm” 1 for A < 660 nm and 40 /xW 
cm” 2 nm” 1 for A > 660 nm. 
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■able 4. Required instrument sensitivities tor Se.WiFS validation and algorithm development as a function of 
radiometric measured variable and wavelength ’ 


Property 


Ed(z,X), 

Downwelled 

Irradiance 


E u (z,X), 

Upwelled 

Irradiance 


Lu 5 

Upwelled 

Radiance 


Variable 


Bd(0) m ax 


dE 

dN 

dE 

dN 


B u (0) ma x 

Eu ( t ^) 


dE 

dN 

dE 

dN 


Beal, Source 
Irradiance 

Leah Source 
Radiance 


L u (O)max 

Lu (/o) 


dL 

dN 

dL 

dN 


B C al 

dE 

dN 

Lq al 
dL 
dN 


410 nm 


488 nm 


300 

1 

5 x 10' 3 
5 x 1(T 2 

120 

37 

1 x 1CT 2 
5 x 1(T 4 
5 x 10" 5 

24 

7.5 

2 x 10' 3 
5 x 1(T 4 

5 x 10~ 5 

2 

2 x 10~ 3 

0.6 

6 x 10" 4 


665 nm 


300 

1 

5 x 1(T 3 
5 x 10" 2 

120 

22 

2 x 1(T 2 
5 x KT 4 
5 x 10 -5 

24 

4.5 

4 x 1CT 3 

5 x 1(T 4 
5 x 10~ 5 

5 

5 x 1(T 3 

1.5 

1 x 10“ 3 


300 

1 

5 x 1(T 3 
5 x 10 -2 

60 

1.5 

1.5 x 10 -3 
5 x 1CT 5 
5 x 1CT 6 

8 

0.3 

3 x 10~ 4 
5 x 10" 5 
5 x 10” 6 

15 

1 x 10 -2 
4.5 

4 x 10 -3 


Comment 


Saturation Irradiance 
Minimum Expected Irradiance 

Digital Resolution (profiles) 

Digital Resolution (surface unit) 

Saturation Irradiance (Case-2/coccoliths) 
Saturation Irradiance (Case-1) 

Minimum Expected Irradiance 
Digital Resolution (surface unit) 

Digital Resolution (profiles) 

Saturation Radiance (Case-2/coccoliths) 
Saturation Radiance (Case-1) 

Minimum Expected Radiance 
Digital Resolution (surface unit) 

Digital Resolution (profiles) 

Calibration Irradiance 

Digital Resolution {Ed, E s , E u cal.) 

Calibration Radiance 
Digital Resolution (L u cal.) 


~ -2 nm 1 sr 


Notes- 1 E u and E d are in units of /iW cnT 2 nm 1 and L„ is in units of /iWcm 

' 2. Responsivity resolution in radiometric units per digital count at the minimum requ.red s.gnal level. 

3. Specified ranges should maintain a 100:1 SNR. 


6. The minimum required irradiances at three opti- 
cal depths (as given in Table 4) assumes minimum 
reflectances of 1% at 410 nm, 2% at 488 nm, and 
0.15% at 670 nm. 

7 The saturation and minimum radiances, and radi- 
ance responsivity resolutions, for L u { 0) and L u (z) 
at three optical depths are 0.2 times the correspond- 
ing specification for E u { 0) or E u (z). This assumes 
Ey./L u = Q = 5 at all wavelengths and depths. 

The specifications in Table 4 are meant as guidance to 
interpret the following required performance requirements: 

a) The instrument must maintain a 100:1 SNR at every 
operating range encountered, during field measure- 
ments. 

b) The data for measurements obtained in the field 
must be recorded with a digital resolution less than 
or equal to 0.5% of reading. 


c) The dynamic range of the instrument s linear sen- 
sitivity must extend to include the signal levels en- 
countered during laboratory calibrations, and the 
calibration signals must be recorded with a digital 
resolution of 0.1% of reading to permit 1% uncer- 
tainty in calibration. 

In general, the above performance specifications do not 
pose exceptionally difficult engineering challenges, with the 
possible exception of the full dynamic range implied by 
Case-2 or coccolith saturation radiance L u (665) to mini- 
mum expected L u (665). In any event, this situation will 
require specially designed radiometers (Section 3.1.8). It 
is not necessary that every radiometer used for SeaWiFS 
validation operate over the full dynamic ranges given in 
Table 4. A radiometer is merely required to maintain the 
above performance specifications over the dynamic ranges 
of irradiance and radiance existing at locations and associ- 
ated illumination conditions where it is used for SeaWiFS 
validation or algorithm development. 
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3.1.3 Linearity and Stability 

Errors attributable to linearity or stability should be 
less than 0.5% of the instrumental readings over the dy- 
namic ranges specified in Table 4. This is a challenging 
goal, but one which must be met if the equally challenging 
goal of achieving 1% uncertainty in absolute calibration is 
to be meaningful. 


3.1.6 Operating Depth 

The instruments shall be capable of operating to depths 
of 200 m. Depths should be measured with an uncertainty 
of 0.5 m and a repeatability of 0.2 m for profiles in bands 
1 - 6 . 

3.1.7 Instrument Attitude 


3.1.4 Sampling Resolution 

Sampling frequency should be compatible with the pro- 
filing technique being used. For the preferred multispectral 
filter radiometers and spectroradiometric (dispersion) in- 
struments using array sensors, the minimum sampling fre- 
quencies are determined by the profiling rate and the depth 
resolution required. In general, five or more samples per 
meter should be obtained at all wavelengths. All chan- 
nels of E d (z, A), E u (z, A), and L u (z, A) at all wavelengths 
should be sampled within 10 -2 s at each given depth. 

The time response of the instrument to a full-scale (sat- 
uration to dark) step change in irradiance should be less 
than one second to arrive at a value within 0.1%, or one 
digitizing step, whichever is greater, of steady state. In 
addition, the electronic e-folding time constant of the in- 
strument must be be consistent with the rate at which the 
channels are sampled, i.e., if data are to be acquired at 
10 Hz, the e-folding time constant should be 0.2 s to avoid 
aliasing. Individual data scans may be averaged to improve 
signal-to-noise performance, provided adequate depth re- 
solution is maintained. 

3.1.5 Angular Response Characteristics 

The response of a cosine collector to a collimated light 
source incident at an angle 6 from the normal must be such 
that: 

1) for E u measurements, the integrated response 
to a radiance distribution of the form L(6) oc 
1 +4 sin 6 should vary as cos 6 accurate to within 
2%; and 

2) for E d measurement, the response to a colli- 
mated source should vary as cos# accurate to 
less than 2% for angles 0° < 6 < 65° and 10% 
for angles 65° < 0 < 85°. 

Departures from cos 8 will translate directly to approxi- 
mately equal errors in E d in the case of direct sunlight. 

The in-water FOV for upwelled radiance bands should 
be approximately 10° (half-angle). The resulting solid an- 
gle FOV (approximately 0.1 sr) is large enough to provide 
reasonable levels of flux, using silicon detectors, yet small 
enough to resolve the slowly varying (with #, for 0 < 30°) 
field of upwelled radiance. Smaller FOV sensors are ap- 
propriate, of course, if all of the other performance sped- 
fications are satisfied. 


The orientations of the instrument with respect to the 
vertical shall be within ±10°, and the attitude shall be 
measured with orthogonally oriented sensors from 0-30° 
with an uncertainty of ±1° in a static mode; it is not in- 
tended that this uncertainty be maintained while an instru- 
ment is subject to large accelerations induced by surface 
waves. These data shall be recorded with the radiometric 
data stream for use as a data quality flag. 

3.1.8 Red and Near-Infrared Wavelengths 

The fact that the SeaWiFS red and near-IR channels — 
bands 6, 7, and 8 at wavelengths of 665, 780, and 865 nm, 
respectively have such short attenuation lengths in water 
requires that special attention must be paid to these mea- 
surements. Problems due to instrument self-shading (Gor- 
don and Ding 1992) and very rapid attenuation of L u (z, A) 
must be considered at these wavelengths. Large instru- 
ments, such as the standard Marine Environmental Ra- 
diometer (MER) packages from Biospherical Instruments, 
Inc. (BSI), are not adaptable to these measurements. 

Suggested procedures for making the measurements are 
to use either fiber optic probes carrying light back to a re- 
mote instrument, or very small single-wavelength discrete 
instruments. Each of these concepts is adaptable to de- 
ployment from a small floating platform. Care must be 
taken to avoid direct shading by the supporting platform, 
but at these wavelengths, the large attenuation coefficients 
of water makes shadowing by objects more than a few me- 
ters away irrelevant. 

The minimum measurement scheme would be two dis- 
crete (10 nm FWHM) channels at 780 and 875 nm Ad- 
ditional channels at 750 and 850 nm, or more elaborately, 
high resolution spectroradiometry, would be useful in de^ 
termining the spectral distribution of the upwelling light 
field in these bands. 

These measurements should be performed as part of 
the standard validation data acquisition, because of their 
importance in the atmospheric correction algorithms. It is 
anticipated that in the majority of cases, and particularly 
in most Case-1 waters, these measurements will show neg- 
ligible upwelling light. In Case-2 waters, cases of extremely 
high productivity, or in coccolithophore blooms, L WN { A) 
at these wavelengths may be significant, and these mea- 
surements will become very important. 

When in-water measurements are performed at these 
wavelengths, the deck cell channels should be expanded to 
include bands at 750 and 875 nm (Table 2). 
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3.2 SURFACE IRRADI AN CE 

The spectral irradiance at the ocean surface shall be 
measured at wavelengths which correspond to the Sea- 
WiFS spectral bands (Table 2), but with 10 nm FWHM 
bandwidth. A total radiation pyranometer may provide 
helpful ancillary information, but this is not a required 
instrument. 

Instruments mounted aboard ships must be positioned 
to view the sky with minimum obstruction or reflections 
from the ship’s superstructure, antennas, etc. Particular 
care must be taken to prevent sun shadows from antennas 
falling on the irradiance collecting surface. Gimbal mount- 
ing of the deck sensor may be helpful to keep the surface 
of the sensor horizontal. Improperly designed gimbal sys- 
tems, however, can accentuate fluctuations caused by ship 
motion, and if there is obvious oscillation in the measured 
irradiance, the gimballing should be improved to eliminate 
the problem. 

An intuitively attractive technique is to measure irra- 
diance with a sensor floated a fraction of a meter below 
the sea surface, far enough away from the ship to avoid 
ship shadows. The flotation assembly should be designed 
to avoid shadowing the radiometric FOV and to damp 
wave-induced motions. This type of arrangement has an 
additional potential for supporting a small sensor to also 
measure upwelling radiance, L u ( A), just below the surface. 
Unfortunately, the oceanographic community has only had 
very limited experience with this approach for measuring 
E s ( A) (Waters et al. 1990) and its attendant difficulties 
with wave- induced fluctuations in near-surface Ed- Addi- 
tional research should be performed to evaluate the use of 
a floating surface unit as the potentially preferred method 
for measuring E s ( A) in future revisions to these protocols. 

3-2.1 Surface Radiometer Characteristics 

The specified number of channels and spectral charac- 
teristics of deck cells are the same as those for subsurface 
irradiance measurements as shown in Table 2. Saturation 
irradiances are the same as for Ed(X) (Table 4). The dy- 
namic operating range for these sensors needs to only be 
25 db, with a SNR of 100:1 but must include the nominal 
calibration irradiance (Table 4). Linearity must be within 
±0.5%. Sampling frequency should match the frequency 
of the underwater radiometer, which should be 1 Hz or 
faster, and all wavelengths should be sampled within an 
interval less than or equal to 10~ 2 s. Cosine response char- 
acteristics should give relative responsivity to a collimated 
source (in air) which matches cos#, accurate within 2% 
for 0° < 6 < 65°, and within 10% for 65° < 6 < 90°. If 
a floating surface radiometer is used, its cosine response 
and immersion characteristics must meet the same specifi- 
cations as those for profiling irradiance meters. 

For some oceanographic process studies, it may be ac- 
ceptable to use a radiometer system measuring E s { X) at 


only a single wavelength. If only a single channel deck 
radiometer is available, its spectral characteristic should 
closely match one of channels 2-5 with a 10 nm FWHM 
bandwidth. A broad-band, or photosynthetically available 
radiation (PAR), radiometer should never be used for this 
purpose. 

3.3 ABOVE- WATER RADIOMETRY 

The performance characteristics to be specified for an 
above-water ocean color radiometer will vary, depending 
on how a particular instrument is to be employed in Sea- 
WiFS validation experiments. For radiometric compar- 
isons with SeaWiFS and in-water measurements, the fun- 
damental criterion to be met is that estimates of spectral 
normalized water-leaving radiance derived from shipboard 
or airborne measurements must have the same uncertainty 
specified for those derived from in-water measurements of 
L u (z, A) (Table 4). A less accurate radiometer may be used 
to semi-quantitatively characterize spatial variability near 
ship stations. 

In general, the spectral characteristics of above-water 
radiometers should match those specified for L U (A) in Ta- 
ble 2. In some cases, however, it may be acceptable for 
a radiometer to match the SeaWiFS specifications, which 
specify center wavelength within 2 nm and 20 nm FWHM 
bandwidth. Recalling the sensitivity of solar radiometry to 
the exact center wavelength and detailed spectral response 
function (Sections 3.1.1 and 4.1.2), any use of airborne 
radiometers must quantitatively account for the different 
spectral responsivity functions between measurements of 
radiance by SeaWiFS, in-water radiometers, and above- 
water radiometers at each channel’s nominal center wave- 
length. 

A high-altitude imaging radiometer must have a radio- 
metric uncertainty and SNR in all channels equal to those 
of the SeaWiFS instrument if its imagery is to be used 
for direct radiometric verification of SeaWiFS radiometric 
performance. In some cases, the requisite SNR may be 
realized through pixel averaging to a 1 km spatial resolu- 
tion commensurate with that of SeaWiFS. Direct radio- 
metric comparisons between aircraft and SeaWiFS radi- 
ances, however, also require that the different atmospheric 
path effects be carefully modeled, and that the uncertainty 
in those modeled adjustments be independently estimated. 
This can be done most effectively when the aircraft, mea- 
surements are combined with the full suite of shipboard 
in- water, atmospheric, and ancillary measurements (Ta- 
ble 1). In this case, direct comparisons between aircraft 
and ship radiometry may require that both the SNR and 
the uncertainties realized in combined analyses of the two 
data sets will represent a smaller spatial resolution than 
the nominal 1 km instantaneous field-of-view (IFOV) for 
SeaWiFS. 

Performance characteristic specifications are similar for 
ocean color radiometers used to measure water-leaving ra- 
diance from either the deck of a ship or an aircraft flown 
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at low altitude, i.e., 200 m altitude or lower. Radiomet- 
ric characteristics should match the criterion set forth for 
in- water L u { A) radiometers in Sections 3. 1.1-3. 1.4 and Ta- 
bles 2-4. The instrument FOV should be between 5° and 
10° (full angle), and all wavelengths must be coregistered 
within 10% of the IFOV. All channels must be scanned 
simultaneously, or within less than 10 _2 s (depending on 
the digitizing design), to avoid aliasing due to varying wave 
reflectance in shipboard measurements, and to avoid time- 
space aliasing in airborne measurements. This constraint 
precludes use of filter wheel radiometers and others which 
scan channels sequentially over a time interval greater than 
10“ 2 s. Sampling over longer periods of time may be done 
by either electronic integration of all channels simultane- 
ously, or by averaging multiple scans. 

A radiometer’s sensitivity to the polarization of aper- 
ture radiance is critical for ocean color remote sensing ap- 
plications. Polarization sensitivity is likely to be present 
in any radiometer having mirrors, prisms or gratings in 
its optical path. To measure accurate water- leaving radi- 
ances using instruments of these types, it is necessary to 
depolarize aperture radiance using either fiber-optics or a 
pseudo-depolarizer . Shipboard and airborne ocean color 
radiometers must have a polarization sensitivity of less 
than 2% in all channels. The sole exception to this rule 
will occur in the case of instruments designed to actually 
measure the polarization components of aperture radiance, 
e.g., the polarization channels of the French Polarization 
and Directionality of the Earth’s Reflectances (POLDER) 
instrument. 

Each application of a particular above-water radiome- 
ter system, if it is proposed for SeaWiFS validation, must 
be evaluated on its own merits. The instrument’s respon- 
sivity, uncertainty, stability, FOV, and spectral character- 
istics must be evaluated in the context of the models to be 
used to compare its radiance measurements to in- water, or 
SeaWiFS, radiance measurements. The suitability of spa- 
tial averaging to improve SNRs must be evaluated in terms 
of the spatial variability prevailing in the experiment site, 
particularly when in-water and aircraft radiances are to 
be directly compared. Finer resolution aircraft imagery, 
or low-altitude trackline data, will often be essential for 
determining the validity of attempts to directly compare 
in-water and SeaWiFS radiances measured at a particular 
site. 

In summary, airborne and shipboard above- water ra- 
diometry can obviously contribute extremely valuable data 
for validating the radiometric performance of the SeaWiFS 
instrument and the algorithms employed with SeaWiFS 
data. There is, however, a wide possible range of radiome- 
ter characteristics that can be applied to this program, and 
detailed specification of required characteristics can only 
be done in the context of each particular experiment’s de- 
sign. Only the guiding principals and desired end-to-end 
performance are specified here. 


3.4 IOP INSTRUMENTS 

The IOPs are: 

1) the beam attenuation coefficient, c(z, A), in units 
of m" 1 ; 

2) the absorption coefficient, a(z, A), in units of 
m -1 ; and 

3) the volume scattering function, b(z,\o,0), in 
units of m _1 sr" 1 . 

The integral of the volume scattering function over 47r 
steradians is the total scattering coefficient, 6(2, A), with 
units of m -1 . The integral of the volume scattering func- 
tion over the back hemisphere is the backscattering coeffi- 
cient, 65(2, A), with units of m” 1 . 

It will be possible to measure the spectral attenuation 
and absorption coefficients in situ at the time of SeaWiFS 
deployment. The instruments for the measurement of the 
spectral absorption and attenuation coefficients should, at 
a minimum, have the characteristics given in Table 5. 

Spectral resolution at more than SeaWiFS wavelengths 
would be desirable to deduce pigment concentrations. In 
the case of beam attenuation coefficients, the requirements 
for uncertainty and precision correspond to changes in c(A) 
resulting from changes in concentration of approximately 
5 and 2 fig 1 _1 of suspended mass, respectively. Stability 
should be tested with instruments connected to the data 
acquisition system. Stability with time should be better 
than 0.005 m -1 between calibrations. 


Table 5. Minimum instrument characteristics for 
the measurement of the spectral absorption and at- 
tenuation coefficients. 


Instrument Characteristics 

Spectral Resolution: 

410, 443, 490, 510, 
555, and 670 nm 

Bandwidth: 

10 nm 

Uncertainty: 

0.005 m- 1 

Precision for A < 650 nm: 

0.002 m- 1 

Precision for A > 650 nm: 

0.005 m- 1 

Stability with 

0.005 m -1 over 

Temperature: 

0-25° C 

Sampling Interval: 

> 4 samples m -1 

Source Collimation Angle: 

< 5 mrad 

Detector Acceptance Angle: 

< 20 mrad 

Depth Capability: 

200 m 


The spectral total scattering coefficient cannot be mear 
sured directly. It can be obtained from 6(A) = c(A — a(A), 
provided c(A) and a(A) are determined with the appropri- 
ate uncertainty. The spectral backscattering coefficient, 
65(A), has the same requirements for spectral resolution, 
bandwidth, and linearity as a(A) and c(A). Since 65(A) is 
not a transmission-like measurement, however, the uncer- 
tainty of its determination will be approximately 10%. 

The shape of the volume scattering function can, at 
present, be determined in situ only crudely with devices 
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like the ALPHA and Scattering Meter (ALSCAT) and the 
General Angle Scattering Meter (GASM), which were built 
more than a decade ago at the Visibility Laboratory of 
the Scripps Institution of Oceanography. These are single 
angle measurement devices, which must be scanned as a 
function of angle and wavelength. Because measuring scat- 
tering with these instruments is a slow process, they do not 
lend themselves readily to incorporation into other instru- 
ment platforms. Since it will be possible to independently 
determine b{ A) and b b ( A) when SeaWiFS is deployed, the 
determination of the shape of the volume scattering coeffi- 
cient could possibly be determined with acceptable uncer- 
tainty by means of measurement of a few moments of the 
scattering function. A new instrument development effort 
would have to be initiated to pursue this approach. 

3.5 ATMOSPHERIC AEROSOLS 

Sun photometers should be used to measure atmos- 
pheric aerosol optical thickness. These sun photometers 
should have specifications in agreement with the World 
Meteorological Organization (WMO) sun photometer spec- 
ifications (Frohlich 1979). Specifically, the instruments 
should have a 2° FOV, temperature stabilization, and a 
precision of ±0.01%. The specific wavelengths of channels 
should correspond to the recommended WMO wavelengths 
of 380, 500, 675, 778, and 862 nm. For SeaWiFS valida- 
tion, additional channels at 410, 440, 490, 510, and 555 nm 
should be added to the WMO set. 

3.6 SPECTRAL SKY RADIANCE 

Measurements of spectral sky radiance distribution 
should be made using a photoelectric all-sky camera. Spec- 


tral characteristics of the sky radiance camera channels 
are those specified for E s ( A) (Table 2). Data should be 
in a format such that absolute radiance values can be ob- 
tained with an uncertainty of 5% and sky irradiance can 
be determined from integrals of the data to within 10%. If 
the dynamic range of the camera is insufficient to capture 
both the sun and sky distribution, neutral density filters 
(or some other method) should be used so that radiance 
from both the sun and sky can be measured. 

3.7 PHYTOPLANKTON PIGMENTS 

HPLC equipment and associated standards must con- 
form to protocols specified in Chapter 9 of the JGOFS 
Core Measurement Protocols (JGOFS 1991). In situ chlo- 
rophyll ffuorometers should have a resolution of at least 
0.001 mg of chlorophyll a per m 3 . 

3.8 HYDROGRAPHIC PROFILES 

A calibrated CTD system should be used to make pro- 
files to maximum depths between 200 and 500 m. The 
instrument should meet the minimum specifications given 
in Table 6. 


Table 6. The minimum instrument characteristics 
for the measurement of hydrographic profiles are 
listed. 


Parameter 

Range 

Uncer- 

tainty 

Resolu- 

tion 

Pressure [dbars] 

0-500 

0.3% 

0.005% 

Temperature [°C] 

-2-35 

0.015° C 

0.001° C 

Salinity [PSU] i 

1-45 

0.03 PSU 

0.001 PSU 
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Chapter 4 


Sensor Characterization 


Introduction 

This chapter details procedures to characterize the sensor performance specifications prescribed in Chapter 3. 
Procedures for characterizing environmental radiometers are presented in Section 4.1 including special char- 
acteristics of underwater radiometers. Section 4.2 describes radiometric standards used as a basis for absolute 
responsivity characterization of radiometers. Special consideration for radiometers measuring Lvv(A) above 
water either from a ship or aircraft, are discussed in Section 4.3. The remaining Sections ,4.4^4.8, address cah- 
TatTons of IOP instruments, sun photometers, radiance distribution cameras, HPLC and fluorometric systems, 

and CTD profiles. 


4.1 RADIOMETRY 

The characterization of radiometric instruments used 
for the acquisition of field data for SeaWiFS validation and 
algorithm development shall include the determination of 
those instrument characteristics that affect its calibration 
as used in the field environment. In addition to the ob- 
vious radiometric calibration, it is therefore necessary to 
determine: 

a) the spectral sensitivities of the various measure- 
ment channels; 

b) the angular sensitivities of an irradiance or ra- 
diance sensor in the medium, i.e., air or water, 
in which it is to be used; 

c) the temporal response of the system; 

d) the effects on responsivity caused by water im- 
mersion; and 

e) the effects of temperature and pressure on the 
above characteristics. 

The elements of radiometer characterization and calibra- 
tion are outlined schematically in Fig. 1. 

For any instrument to provide suitable data for SPO 
use, the investigator must be certain that the instrument 
characterization has not changed beyond accepted limits 
and that the time history of the calibration is traceable. 
Certain attributes, e.g., angular response characteristics, 
are sufficiently constant that they only need to be deter- 
mined once, unless the instrument is modified. The ex- 
act nature of instrument modifications during maintenance 
will determine which characterization procedures must be 
repeated. On the other hand, radiometric calibrations and 
the assessment of system spectral characteristics of filter 
radiometers, must be repeated before and after each major 
field deployment. 


4.1.1 Absolute Radiometric Calibration 

Determination of the absolute radiometric responses of 
the irradiance and radiance sensors requires the availability 
of a properly manned and equipped radiometric calibration 
facility. Such a facility must be equipped with suitable 
stable sources and sensors, e.g., lamp standards of spec- 
tral irradiance and flat response radiometers, respectively. 
Either the sources or the sensors must have defined spec- 
tral radiometric characteristics that are traceable to NIST 
(Section 4.2). The calibration facility must also have a va- 
riety of specialized radiometric and electronic equipment, 
including: reflectance plaques, spectral filters, integrating 
spheres, and highly regulated power supplies for the opera- 
tion of the lamps. Precision electronic measurement capa- 
bilities are also required, both for setting and monitoring 
lamp current and voltage and for measuring the output of 
the radiometer. 

It is not expected that every investigator will be able 
to perform his own radiometric calibrations; because of 
this, a few centrally located facilities will be equipped and 
staffed to perform these calibrations as a routine service for 
the community. The facilities will perform frequent inter- 
comparisons to assure the maintenance of the radiometric 
traceability to the NIST standard of spectral irradiance. 
The goal shall be to provide reproducible calibrations from 
400-850 nm to within better than ±1%; the minimum re- 
quirement for radiometric data to be used in SeaWiFS val- 
idation is for repeatable calibrations within less than 5%. 

Radiometric calibrations of irradiance sensors will be 
performed after it has been ascertained that: the confor- 
mity of the sensor angular response to the required cosine 
function is satisfactory, the sensor linearity is satisfactory, 
and the spectral sensitivity, including out-of-band block- 
ing, is known and satisfactory. 

Radiometers shall be calibrated using a 1,000 W FEL 
standard of spectral irradiance, with calibration traceable 
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Fig. 1. Elements of radiometer characterization and calibration. 


to NIST and lamp operation in accordance with Walker 
et al. (1987). The irradiance collector is placed normal to, 
and at the prescribed distance from, a working standard 
lamp of spectral irradiance. The lamp should be of appro- 
priate size to provide an irradiance at the sensor that will 
be at least 30%, and preferably above 50%, of full-scale for 
the sensor channel being calibrated, although this is not 
always achievable in practice (Table 4). The lamp-sensor 
space shall be appropriately baffled and draped so that oc- 
culting the direct path between lamp and sensor will result 
in a response of less then 0.1% of the response to the lamp 
flux. 

For multispectral instruments, all channels may be cal- 
ibrated simultaneously if sufficient flux is available at all 
wavelengths. The instrument response is recorded for all 
channels together with associated dark responses. Am- 
bient and photosensor temperatures are recorded, where 
available. For characterization, the radiometric calibration 
should be performed at temperature extremes of -2° C and 
40° C for in- water sensors, and at -10° C and 45° C for 
irradiance sensors used above the surface. If responses dif- 
fer significantly at temperature extremes, responses should 
also be determined at intermediate temperatures. 

The portable irradiance and radiance reference stan- 
dard to be used to trace instrument stability during field 
deployments (Section 4.2.5) should be placed in position 
on the sensor immediately following the calibration to es- 
tablish the instrument response to this reference unit. 

Radiance calibration activities require a uniform source 
of known radiance that will fill the angular field of view of 
the radiance sensor. The two procedures that may be used 
are given below. 

Calibration Methods 

1. A working lamp standard of spectral irradiance is 
placed at the prescribed distance from a plaque of 


known Lambertian reflectance. The plaque is nor- 
mal to, and centered on, the lamp calibration axis. 
The radiance sensor is positioned to view the plaque 
at an angle of 45° from the plaque normal (any other 
angle at which the diffuse reflectance of the plaque 
is known is acceptable also). It must be established 
that the plaque fills the sensor’s FOV and that the 
presence of the sensor case has not perturbed the 
irradiance on the plaque. The instrument response 
and dark signal is recorded. It must be verified 
that the plaque fills the FOV with uniform radi- 
ance for each channel of a multichannel radiance 
sensor. Separate calibration setups may be required 
for different channels and the lamps may have to 
be moved as much as 3 m away from the plaque to 
assure uniform illumination. This procedure is dif- 
ficult to apply to sensors with a large FOV. 

2. An integrating sphere with an exit port of sufficient 
size to fill the FOV of the radiance sensor may be 
used if the radiance of the exit port, at the chan- 
nel wavelengths, can be determined with sufficient 
uncertainty. 

These methods are discussed more fully in Section 4.2.2. 

4.1.2 Spectral Bandpass Characterization 

These instruments should be characterized to define the 
nominal wavelengths and bandwidths, defined as the full 
width of the passband as measured to the FWHM intensity 
points. The nominal, or center wavelength, will usually be 
defined as the wavelength halfway between wavelengths at 
which the normalized response is 0.5, and the channel is 
characterized by this wavelength and the FWHM band- 
width. The determination of the spectral response func- 
tion, i.e., the passband, will be made for each channel with 
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a scanning monochromatic source, with a bandwidth less 
than 0.2 nm; the source output must be normalized to a 
detector of known spectral sensitivity. The response func- 
tion thus measured is then normalized to the maximum 
(peak). 

Although the results of this characterization are usually 
represented by only the nominal wavelength and FWHM 
bandpass, the full normalized response function should 
be recorded for use in detailed wavelength adjustments 
and comparisons with the SeaWiFS channel response func- 
tions, which will not be known until shortly before launch. 

It is further recommended that the internal instrument 
temperature be monitored during these tests, and that 
the test be repeated at two temperatures at least 15° C 
apart, e.g., 10° and 25° C. If a significant shift, greater than 
1.0 nm, with temperature of either the center wavelength 
or bandwidth is detected, then additional temperature cal- 
ibration points are recommended. Dark offsets must be 
recorded during each test. 

For spectral characterizations of irradiance diffusers, 
the entire surface of the diffuser should be illuminated by 
the monochromator’s output. In the case of radiance de- 
tectors, a diffuser should be used to diffuse the monochro- 
mator slit image and uniformly fill the instrument’s FOV. 

The wavelength response of a monochromator- based 
radiometer is calibrated by scanning over line sources, with 
sharp peaks at well known wavelengths. Suitable spectial 
calibration sources, such as, mercury, cadmium, and neon 
lamps, are provided by several vendors, together with tab- 
ulations of the wavelengths of the emission lines generated 
by each source. 

The width of the slit function of a monochromator may 
be estimated by scanning over a laser line, e.g., helium- 
neon, at a very small wavelength interval. I he instrument 
FOV must be filled during the test. 

It is anticipated that the monochromator- based spec- 
tral characterization will not be able to adequately measure 
leakage of broadly distributed out-of-band radiation; there- 
fore, blocking of blue light in channels longer than 540 nm 
must be routinely tested. Where continuous wave (CW) 
argon lasers are available, out-of-band response should be 
measured at 488 nm. One recommended test that can be 
performed during the absolute calibrations at A < 640 nm 
is the sequenced measurement of three Schott BG-18 fil- 
ters, each 1 mm thick, using an FEL-type light source. The 
procedure is to measure the channel signal using each filter 
separately, then in combination, and comparing the com- 
puted and measured transmissions. If a significantly higher 
combined transmission of the three filters, when they are 
used in combination, is measured relative to the calcu- 
lated transmittance, then spectral leakage is present. At 
wavelengths greater than 640 nm, other filters that atten- 
uate the wavelength of interest, with a transmission value 
of less than or equal to 0.1 and which pass shorter wave- 
length light with significantly greater transmission, should 
be substituted for the BG-18. 


Consideration must also be given to unblocked fluores 
cence by the filters, or other optical elements, as a possi- 
ble source of light leaks. Methods to test for fluorescence 
contamination specifically, are not well established at this 
time. 

While leakage of blue light into red channels is the most 
significant oceanographic optical problem, the leakage of 
red and IR light into blue channels can cause significant 
errors when the instrument is calibrated using a red-rich 
source. A convenient way to measure this leakage is to 
place a long wavelength- pass, sharp-cut, absorbing glass 
filter that does not exhibit fluorescence between a broad 
band (e.g., incandescent) source and the sensor. A non- 
zero response indicates unwanted out-of-band red response 
and the need for improved red blocking. 

4.1.3 Temporal Response 

The temporal response of a spectrometer may be ex- 
amined by introducing a step function of near full-scale 
flux to the system using an electrically operated shutter 
and measuring the system’s transient response at 0.1 s, or 
shorter, intervals. The response should be stable within 
one digitizing step, or 0.1%, whichever is greater, of the 
steady state value in one second or less. 

4.1.4 Radiance Field-of-View 

It is required that the radiance FOV of the instrument 
be known. The FOV should not normally enter into the 
absolute calibration, however, if the FOV is fully filled by 
a calibration source of uniform radiance. 

In this test, the instrument is placed on a rotational 
stage with the entrance aperture of the radiometer over the 
rotation axis. A stable light source with a small filament 
is placed several meters in front of the instrument, which 
is then scanned from —30° to +30° in 2° increments. The 
angle positioning should be within ±0.1°. The on axis, 0°, 
mechanical alignment is made using the window surface 
as reference, by adjusting to get the reflection of the lamp 
filament to return on axis. The error in this alignment 
is approximately ±0.1°. The in-air measurement angles, 
0 O , are converted to corresponding angles in seawater, 0 W1 
using the relation 0 W = 0 a jn w , where n w is the index of 
refraction of seawater at the particular wavelength of each 
channel. 

4.1.5 Collector Cosine Response 

The directional response of cosine collectors must be 
characterized. The directional response of the deck cell is 
determined in air, and the in-water instruments are mea- 
sured immersed in water. Full spectral determinations are 
required. For instruments measuring upwelling irradiance 
E u (z, A), it is recommended that the cosine response of 
each instrument be measured individually. For downwell- 
ing irradiance Ed(z , A) instruments, checking a production 
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run may be satisfactory if the vendor’s material and design 
are demonstrated to be uniform throughout the duratio 

° f t Absolute responsivity calibration of an irradi ^ nce 
ter is done in air, using light incident normal to the co lec- 
tor To properly measure irradiance incident on the pla 
at all angles /(relative to the normal), the instrument s 
response should follow a cosine function. In other wor s, 
for an instrument response V(0) to a given collimated ir- 
radiance incident at 0 = 0, if the ins trument is routed 
to the angle 0 away from the original normal axis, the re 
sponse should be V(8) = V(0)cos8. If this ^™ ent J 
met then the on-axis calibration is sufficient and the de- 
vice’ will correctly measure irradiance arriving at the plan 
of the collector, regardless of the directional distribution 

at which the light arrives. . . 

The preferred irradiance collector design 
proved cosine response over that of a simple flat plate d - 
fuse collector (Boyd 1955 and Tyler and Snnth lOT^ J us 
improvement is mostly for near-grazing angles (6 approach 
ing 90° to the normal) and is particularly impor 
measurements of the upwelling underwater irradiance are 
made, i.e., with the collector facing downward. In that 
case, most of the light is from the sides, in the region 

these near-grazing angles. , 

Since EJz, A) measurements are to be made u 
ter the tetttg to determine the fidelity of the instrumen 
to ihe cosine function must be made with the instrument 
submerged. A description of the suitable experimental pr 
cedure follows Petzold and Austin (1988). 

The instrument is suspended in a tank of water while 
supported by a fixture designed to allow rotation about 
an axis through the surface and center of the collector A 
tungsten-halogen lamp with a small filament is enclosed 
a housing with a small exit aperture and 
mately 1 m from a large window in the tank. The collec 
tor is placed approximately 25 cm behind this window, an 
equivalent lamp distance of 1.25m or more is required. A 
circular baffle should be placed immediately m front of the 
window to reduce stray light. The water should be hig y 
filtered to the extent that the effects of scattered light are 

mdiscerniquivaient air path lamp dista nce should be ap- 
proximately 1.25 m or greater. At this distance, the a- 
off at the outer edge of a 6 cm diameter diffuse collector 
would be 0.9994, or -0.06%, when the diffuser is at 8 - 0 
with the normal. The net effect e ? ire area of 

the diffuser would be 0.9997 or -0.03%. When 8 - 90 , 
with the diffuser edge-on to the lamp the distance to the 
lamp varies for different points on the surface. T 
error over the entire surface for this condition is 0.99997 
or -0.003%. All other angles fall between these limiti g 

C3S The signals from the instrument are recorded for 8 = 0° 
and at 5° intervals to 8 = ±75° and 2 5° intervals over 
75° < 8 < 90°. The readings at 8 = 0° are recorded at 


the beginning, the middle, and the end of each run and 
examined a/a measure of lamp and -«„« stab,!,., 
over the time involved. At least two runs should be made 
about different axes through the surface of the diffuse^ 
All readings are normalized to 1.000 at 6 an 

compared with the value of the cosine of each angle. If 
V(9) is the normalized measured value, relative local error 

at angle 0 is given as (V(0)/cos0) -l- 

Assuming the average response to the four measur r - 
ments made at each 8 (four separate azimuth ang es <p) 
adequately represent the overall mean cosine respond , o 

the collector, then the error, £, in measuring irradiance 

over the interval 8 n < 8 < 8 N for a uniform radiance dis- 
tribution is approximately 


N 


£ V(8i) sin 0 t A0 

_ — 1, 

e ~ N 

£ cos 8 t sin 0, A0 


( 1 ) 


using a simple trapezoidal quadrature. Similarly, for a 
radiance distribution of the form 1 + 4sm«, to Simulate 
upwelled irradiance 

£ V(0i) (1 + 4 sin 0.) sin 8 l A 0 
_ t=o 1, (2) 

£ cos0 i (l + 4sin0 l ) sin0, A0 
i=0 

where 0o = 0, 0 n = 2 and — 2 N' Pn 11 i valent 

The asymmetry of the cosine response, 6, is equivalent 

to an effective tilt of an ideal cosine collector with respec 

to the instrument’s mechanical axis, which can be quanti- 

fied as ^ 

f cos(0 4- 0t) si n 0 

, (3) 


6 = 


#2 

f cos (0 
0\ 


■ 9 t ) sin 0 d6 


where 0t is the tilt angle. . r 

The measured asymmetry is computed as the ratio of 

sums of measurements at opposite <t> (0 > 0) and tt ( 

0) in the same plane, that is, 


0 N = 7t/2 

V(0i,O) sin 0i A0 

c _ i=0 

0 0 N =_7r/2_ 

^2 V(0 l )sin0 l A0 

t=0 


1, 


(4) 


for A0 = ±277* 


Variations in asymmetry from channel to channel may 
be due to the placement of the individual detectors behin 
the diffuser. Any offset of the average asymmetry with th 
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mechanical axis could be due to any one of a variety of 
causes: 

1) the alignment on the rotating test fixture not 
being correct, 

2) tilt of the diffuser, 

3) the detector array not being centered, 

4) nonuniformity of the reflectance of the internal 
surfaces of the instrument between the diffuser 
and the sensor array, or 

5) nonuniformity of the diffuser. 

4.1.6 Immersion Factors 

When a plastic, opal-glass, or Teflon diffuser is im- 
mersed m water, its light transmissivity is less than it was 
in air. Since an instrument’s irradiance responsivity is cali- 
brated in air, a correction for this change in collector trans- 
missivity must be applied to obtain irradiance responsivity 
coefficients for underwater measurements. 

The change in a collector’s immersed transmissivity is 
t e net effect of two separate processes: a change in the 
reflection of light at the upper surface of the collector, and 
internal scattering and reflections from the collector’s lower 
surface. A small part of the light flux falling on the col- 
lector is reflected at the air-plastic, or water-plastic, inter- 
ace, and the majority of the flux passes into the collector 
body. The relative size of this reflectance, called Fresnel 
reflectance, depends on the relative difference in refractive 
indices between the diffuser material and the surrounding 
medium. & 

The refractive index of the collector material is always 
larger than that of either water or air, and because the 
refractive index of water is larger than that of air, Fresnel 
reflectance is smaller at a diffuser-water interface than at 
a diffuser-air interface. Thus, the initial transmission of 
ig . rough the upper surface of an irradiance collector 
is larger m water than in air. The immersed upper sur- 
face is, on the other hand, also less effective at reflecting 
e upward flux of light backscattered within the diffuser 
body and light reflected at the lower diffuser-air interface 
m the instrument’s interior, processes that are not affected 
y immersion. Therefore, a larger fraction of the internally 
scattered and upwardly reflected light passes back into the 
water column than would be lost into air. Because the in- 
creased upward loss of internally reflected flux exceeds the 
gain in downward flux through the diffuser-water interface 
the net effect of these competing processes is a decrease in 
the collector s immersed transmissivity. 

To measure this effect, a suggested and acceptable pro- 
cedure is as follows: The instrument is placed in a tank of 
water with the irradiance collector level and facing upward 
A tungsten-halogen lamp with a small filament, powered 
by a stable power supply, is placed at some distance above 
e water surface. The depth of the water is lowered in 
steps and readings are recorded for all wavelengths from 
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each carefully measured depth. A final reading is taken 
with the water level below the collector, i.e., with the col- 
lector m the air. The amount of energy arriving at the 
collector vanes with the water depth and is a function of 
several factors: 

a) the attenuation at the air-water interface, which 
varies with wavelength; 

b) the attenuation over the water pathlength, which 
is a function of depth and wavelength; and 

c) the change in solid angle of the light leaving the 
source and arriving at the collector, caused by 
the light rays changing direction at the air- water 
interface, which varies with wavelength and wa- 
ter depth. 

Usi [ lg ( Fresnel reflectance equations, the transmittance 
through the surface is 


Ts( X) = 


4nu,(A) 


(l + n u ,(A)) 2 


( 5 ) 


where n w (\) is the index of refraction of the water at wave- 
length A. 

. The transmittance through the water path, 7 UI (A), will 
T W (X) = (6) 

where K(\) is an attenuation coefficient of the water and 
2 is the path length in corresponding units. 

The change in solid angle with water depth x is given 
by the factor 


G(z, A) = 



-2 


( 7 ) 


where d is the distance of the lamp source from the collec- 
tor surface. 

The immersion correction factor F t ( A) for irradiance is 
then calculated for each depth z as 

Fi ^ = A) Ta ^ Tw A ) G ( 2 - A ) - (8) 


where E a ( A) and E w (z, A) are the irradiance in air and the 
irradiance underwater at depth z, respectively. 

There are two unknowns in (5)-(8): the attenuation co- 
efficient of the water K( A) and the immersion factor FAX) 
A minimum of three measurements must be made to solve 

H t0 g6t E *W’ and two different 

water depths for E w (z, A). The recommended method is 

to take readings of E (z,X) at many depths. Then, using 

for the ^7“ d im a IeaSt ~ squares egression is solved 
for the F t ( A) and K( A) terms giving the best fit. The com- 

plete jderivatioH of (5)-(8) is given in Petzold and Austin 
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The absolute calibration for the spectral radiance chan- 
nels is found by viewing a surface of known radiance in air 
in the laboratory. When the instrument is submerged in 
water, a change in responsivity occurs and a correction 
must be applied. This change in responsivity is caused by 
the change in the indices of refraction of the different me- 
dia in which the instrument is immersed — in this case air 
and water. Two optical changes occur, both of which are 
caused by the change in refractive index. The two effects 
to be corrected are: 

1) the change in transmission through the inter- 
face between the air and the window during cal- 
ibration, and the same effect through the water- 
window interface during data measurement, and 

2) the change in the solid angle included in the un- 
derwater FOV relative to that in air. 

Since the refractive index of seawater, n u; (A), is a function 
of wavelength A, the correction factor F l (X) will also be a 
function of wavelength. 

If the refractive index of air is assumed to be 1.000 at 
all wavelengths, and if %(A) is the index of refraction for 
the (glass) window and n^A) is the index of refraction for 
water, then, as shown in Austin (1976), the correction for 
the change in transmission through the window, T g { A), is 

r s « _ i ^> +-<*»;„ (9) 

(l 4- n 5 (A)) 

and the correction for the change in the FOV, F v , is 

Fy( A) = (n™(A)) 2 (10) 

The index of refraction of a Plexiglas™ window, n g ( A), 
may be computed using an empirical fit to the Hartmann 
formula, that is, 


n s(A) 


1.47384 + 


7.5 

A - 174.71’ 


( 11 ) 


where A is the wavelength in nanometers (Austin 1976). 
The index of refraction for seawater n^(A) may be similarly 
computed using an empirical fit of the data from Austin 
and Halikas (1976), 


n^(A) = 1.325147 + 


6.6096 

A - 137.1924 


( 12 ) 


4.1.7 Linearity and Electronic Uncertainty 

The linearity of the radiometric channels must be de- 
termined over their expected range of use. The above- 
surface (deck cell) and underwater irradiance sensors in- 
tended for the measurement of downwelling irradiance have 
full-scale (saturation) values that are not readily obtained 
with the usual incandescent blackbody sources, such as 
1,000 W 3,200 K tungsten-halogen projection lamps. The 
linearity at the high end of the calibrated range may be 
determined by using 900-2,000 W high pressure xenon arc 
lamps, which provide a small, stable source of high inten- 
sity (approximately 6,000 K) radiation. With such lamps, 
irradiance levels approximating full sunlight can be at- 
tained. Using such sources for the high end, and the more 
easily managed tungsten-halogen lamps over the range be- 
low 20-30% of full scale, the linearity of the response char- 
acteristic of the radiometric channels can be assessed. The 
flux should be changed in 5 db (0.5 log), or less, steps using 
a proven and accepted procedure for controlling irradiance 
such as inverse square law, or calibrated apertures. These 
suggested procedures for testing linearity at the higher lev- 
els are not well established in practice, and research is 
needed to determine the precision which can be attained. 

If departures from linearity are found, they must be in- 
corporated into the calibration function for the instrument 
and be properly applied to the raw level -1 data to obtain 
calibrated level -2 irradiance and radiance data. Level- 1 
and level-2 data are defined in Section 6.1. 

It is recommended that all instruments utilizing inputs 
from ancillary sensors, e.g., transmissometers, be charac- 
terized for the linearity and uncertainty of the voltage mea- 
surement covering the full output range of the ancillary 
sensor. For instruments with range dependent gain chang- 
ing, either manual or automatic, the scale offset and linear- 
ity for each range should, at a minimum, be tested annu- 
ally. Uncertainties exceeding 0.1% of any reading within 
the normal working range must be investigated and cor- 
rected. 

Other characteristics of electronic sensor systems may 
adversely affect measurement uncertainty. During the de- 
sign and engineering prototype development of a radiome- 
ter, the design and implementation must be analyzed to 
characterize, and correct as needed, possible effects of hys- 
teresis, overload, recovery times, cross talk between either 
optical transducers or electronic channels, and sensitivity 
to orientation in the Earth’s magnetic field, which is par- 
ticularly likely with photomultiplier tubes. 


The immersion factor F t ( A) is then obtained as 


4.1.8 Temperature Characterization 


Fi( A) - T g (X)F v (X) 

__ n ™(^) { n w (A) 4- n 9 (A)) 
(1 + n g ( A)) 2 


Two major types of temperature-induced variation may 
, . be seen in an optical radiometric instrument: 1) offset or 

dark changes, and 2) scale responsivity changes. Each 
underwater instrument must be individually characterized 
over the range of -2-40° C. In the case of deck cells, the 
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temperature range for testing should be extended to -10- 
45° C. Sensors exhibiting temperature coefficients greater 
than 0.01% per °C over this temperature range, should be 
fully characterized over their respective ranges to estab- 
lish the means and precision with which post-acquisition 
processing can be used to correct for temperature depen- 
dency. Although knowledge of the zero, or dark current, 
drift is essential for working at the lowest radiances or ir- 
radiances, it should be emphasized that more significant 
near-surface errors may be induced by temperature varia- 
tions in responsivity. These possible responsivity changes 
must be individually determined across the spectrum. 

In the above discussion, the temperatures cited are en- 
vironmental temperatures, but it should be emphasized 
that any correction must use the temperature of the af- 
fected element, which is normally in the interior of the 
instrument. This is best accomplished by routinely using 
temperature sensors placed at critical locations within the 
instrument. For highest precision, dynamic temperature 
testing involving temporal transients, as well as possible 
temperature gradients within an instrument, may be ap- 
propriate. 

4.1.9 Pressure Effects 

Pressure can cause radiometric measurement errors by 
deforming irradiance collectors. Pressure coefficients asso- 
ciated with polytetrafluoroethylene (PTFE) based irradi- 
ance diffusers are known to exist, but they are not uniform 
and there may be hysteresis effects. It is recommended 
that each type of irradiance detector be examined for vari- 
ations in responsivity with pressure. If a significant effect 
is observed, then pressure-dependent responsivity coeffi- 
cients should be determined separately for each instrument 
and collector. The pressure characterization should also 
test for, and quantify, hysteresis and temporal transients 
in responsivity under a time varying pressure load. The 
characterization of pressure effects has not previously been 
common practice, and the requisite procedures are there- 
fore poorly defined; new protocols must be developed. 

4.1.10 Pressure Transducer Calibration 

The radiometer’s pressure transducer, which is used to 
measure instrument depth during profiles, should be tested 
and calibrated before and after each major cruise. 

4.2 RADIOMETRIC STANDARDS 

This section describes sources and methods by which 
the NIST scale of spectral irradiance is transferred to cali- 
brate irradiance and radiance sensors. The principal work- 
ing standards used for spectral irradiance responsivity cal- 
ibration are FEL lamp working standards (Section 4.2.1). 
The spectral irradiance scales of the FEL lamps are in 
turn transferred to spectral radiance scales using plaques 
of known bidirectional reflectance, or integrating spheres, 


or both (Section 4.2.2). An ongoing series of SeaWiFS In- 
tercalibration Round-Robin Experiments (SIRREXs) has 
been initiated by the SPO to assure internal consistency 
between the laboratories which calibrate radiometers for 
SeaWiFS validation (Mueller 1993 and Mueller et al. 1994). 

4.2.1 Lamp Irradiance Standards 

The options available for radiometric calibration stan- 
dards are limited to standard sources or standard detec- 
tors. Lamp standards of spectral radiance and irradiance 
are provided by NIST and various commercial standard- 
izing laboratories and manufacturers who furnish NIST 
traceable secondary standards. The uncertainty cited for 
these standards by NIST is, at best, 1% in the visible and 
3% is a more realistic estimate of absolute uncertainty at- 
tainable using lamp standards alone. Over the calibration 
range from 250-2,500 nm, the uncertainty is approximately 
6% at the endpoints. 

The lamp standard of spectral irradiance is tradition- 
ally used for radiometric calibration, mainly because of its 
ease of use compared to the spectral radiance lamp. NIST 
publishes guidelines for the setup, alignment, and use of 
these standards. The vendors that manufacture and cali- 
brate these lamps also issue guidelines for their use. 

4.2.2 Radiance 

Spectral radiance may be obtained by using an irradi- 
ance standard lamp and a Lambertian reflecting plaque. 
The standard lamp is positioned on axis and normal to 
the center of the plaque at the calibrated distance. The 
instrument or detector package to be calibrated is nomi- 
nally positioned to view the plaque at 45° measured from 
the axis. The radiance, then, is given by 

L( A) = -p(\)E(\), (14) 

7T 

where p( A) is the bidirectional reflectance of the plaque for 
0° incidence and 45° viewing, E( A) is the known spectral 
irradiance from the lamp during calibration and the total 
FOV of the instrument being calibrated is filled by the 
illuminated plaque. 

The known radiance of the plaque provides an uncer- 
tainty comparable with that of the irradiance standard 
lamp, i.e., less than or equal to 3%, for calibrating a ra- 
diance detector with a very narrow FOV (~1°). Large 
plaques, e.g., 40 cm 2 , have been successfully used to cali- 
brate radiance sensors having up to 25° full- angle FOVs. 
Intercomparisons of calibrations on underwater radiance 
sensors (possessing in-air full- angle FOVs ranging from 20- 
24°), made using this technique at different laboratories, 
have generally agreed within approximately 5%. 

A better approach to calibrating multispectral radiance 
sensors is to view an integrating sphere that is uniformly il- 
luminated by stable, appropriately baffled lamps, and that 
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i? 5 FOV n POr h‘ lar8e : n °“8 h O-Mew, fill the ten- 
sor s FOV. The sphere and exit port must be large enough 

to place the radiance sensor far enough away to prevent sfg 

nificant secondary illumination of the sphere walls due to 

S ot r off : he T or ’ s entrance ^ if the -- 

distort the retro -'-eflected light will both increase and 
distort the uniformity of the radiance distribution within 
the sphere. Traditionally, the calibration of an integrating 

atelvT r ce "T* has bmi »««mplished by appropri 

" OU,P,,, ^ a 8t “ d “ d lam P 


Sphere Calibration Methods 
1. The approach used at NASA/GSFC is to view 
the irradiance output of the lamp, initially, and 
then the sphere, with a spectroradiometer equip- 
ped with integrating input optics (McLean and 
Guenther 1989 and Walker et al. 1991). The 
spectral irradiance responsivity of the radiome- 
ter is calibrated using the lamp data, and the 
(assumed) Lambertian radiance of the sphere is 
determined by dividing the measured spectral 
irradiance output of the sphere by tt. 

2. An alternative method is to calibrate a stable 
narrow FOV radiometer by viewing the stan- 
dard lamp output reflected from a plaque, as de- 
scribed above. The output from the sphere’s exit 
port is then viewed within this radiometer The 
radiometer should also be used, at this point, to 
map the angular distribution of radiance in the 
sphere as viewed through the exit port. This im- 
portant verification of a uniform radiance distri- 
bution is not possible when Method 1 is used to 

nf M r ft. e n P o ,ere radiance A Promising variant 
of Me hod 2 is to calibrate the sphere using a 
sell-calibrating radiometer (Palmer 1988). 

4.2.3 Radiance Standardization 

Detectors of the type embodied in the United Detec- 
Effic eC t n n 0gy ^ ED ' 200 ( w here QED stands for Quantum 

h T ter are 99 " % f l uantum efficient. 
Palmer (1988) shows how such a detector may be com- 

ned with precision apertures and well-characterized fil- 

rs to measure self-calibrated spectral radiance with an 

absolute uncertainty less than 1%. A calibration approach 

based °n such radiometer standards is essential to achieve 

1% internal consistency in the radiometric uncertainty of 

measurements made for SeaWiFS radiometric validation. 

It is worth emphasizing here that the essential objec- 

tkll dat r iu lntemal consistenc y ^ the SeaWiFS op- 

t^hnkui r d g Unif ° rm application of calibration 
chmques based on a common radiometric standard with 

tantf ? appr ° aching 1% > or less if possible. An impor- 
, u not essential, goal is to establish NIST traceable 
absolute uncertainty of less than 1% with this standard 


calihraT H ,br f radiometer ma y be used directly to 
librate and map the radiance distribution of integrating 

sphere sources (Method 2 in Section 4.2.2 above). The self 

ferred a to n a S Standard of radia nce may be trans- 

ferred to a stable lamp source of irradiance through the re- 
fer I °kk reflectance plaque technique, described above 

diance* source* ' * ,1SOr8 * S,a " d “ d '“»P i™- 

tica^T "',77 hava "°* yet b “" into a prac- 

tical and widely accepted set of procedures for the calibra- 
tion of oceanographic or airborne radiometers using self- 
calibrating radiometric standards. A significant level of 
laboratory work must be done to establish the repeatability 
results attainable through these techniques under a va- 

iS ° n °? S ’ ^ t0 C ° dify that ex P (, rience into cal- 

200 tvn ? d f t0 f° S; SpectraI res Ponsivity of the QED- 
200 type detector ,s known, for example, to vary systemat 

ica fly w lth temperature (Kohler et al. 1990), and the spec 

" fanctions of the fibers in a self-calibrating 

radiometer must be re-characterized at a frequency that 

wfl guarantee the uncertainty of the calculated radiance. 
This frequency must be established through experience 
ut a reasonable first approximation is that filter trans- 
mission functions should be remeasured every few months. 

ne goal of the SeaWiFS mission is to base radiometric 
validation on shipboard, moored, and airborne radiometry 

achieved T'T*' ** * t0 b<? subst -tially 

achieved, then the work described above to establish new 

calibration standards and protocols must be pursued vig- 
orously over the next two years. S 

4.2.4 Traceability and Comparisons 

qilrp The r iety ° f instrum ents available for validation mea- 

brltkn? ma r k imperative that «^e common cali- 
bration traceability exists. Recognizing that it would be 

f —graphic 

tion facilities should be established,’ b 

and protocols used at these facilities should all be traced 
directly to those at the GSFC calibration facility Tht 
organizational structure is shown schematically in Fig 2 

int ; rCOmparison ma y include us'e of 
T calibrated filter radiometers to track and document 

this°int r i 6aCh faf:ll ‘ ty (^diometer wavelengths for 

Mnd cZ°rT nS ° n WiU ^ det ™ ed )- Round-robin 
blind calibration comparisons of a standard instrument 

sistenev of to b-Klmuu* ‘be internal co„- 

involvS * performed at ‘be various facilities 

4.2.5 Portable Standards 

Between radiometric calibration activities, stable lamp 
sources in rugged, fixed geometric configurations shoulZ 
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structure lor optical instrumentation characterization and calibration. 

not a substitute for complete calibrations. The .temporal 
record they provide will, however, be invaluable m cases 
where the pre- and post-cruise calibrations disagree or 
the instrument is disturbed, e.g„ opened between calibra- 
tions or if the data quality are otherwise suspect. These 
portable standards are an important part of the recom- 
mended instrument package. 

Although several manufacturers offer somewhat por 
able irradiance and radiance sources, there has be <^ 
little previous work to validate and use portable radiomet- 
ric standards to test oceanographic radiometers in thefiekb 
Therefore, detailed hardware specifications and procedura 
protocols must be developed through a series of laboratory 
and field tests using candidate equipment and standards. 

4.3 ABOVE- WATER RADIOMETRY 

In general, the protocols specified in Section 4.1 for 
in- water oceanographic radiance sensors are app ica e or 


used to track instrument performance. Irradiance chan 
nels can be monitored with irradiance sources at fixed di 
tances from the collectors, while radiance sources can be 
monitored by filling the FOV with 

in front of the irradiance sources, or by using integrati g 
cavity sources. In each case, careful attention mus 
given to fixing specific geometries of source and debtor 
in each use. The stability of the lamp output and the 
repeatability of measurement must be sufficient to e ec 
2% variations in an instrument's performance. An instra- 
ment should be connected to the portable standard and 
response recorded daily, keeping a record of instrument r^ 
sponsivity throughout an experiment. * 

sources would provide an essential warning of problems 

^The^ortable field reference source must be available 
when the complete radiometric calibrations are Performed 
so that a baseline may be established and maintained f 
each sensor channel (Section 4.1.1). T ese sources ar 
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characterizing and calibrating shipboard and airborne ra- 
diometers used to measure water-leaving radiance. Obvi- 
ous exceptions are that immersion and underwater FOV 
characterizations are not appropriate for above-water sen- 
sors. 

Polarization sensitivity is more critical in above-water 
radiometry than underwater radiometry. If a radiometer 
measures polarization components of radiance, then its 
responsivity and rejection of cross-polarization radiance 
must be characterized for each component channel. For 
above-water scalar radiance instruments, as with the Sea- 
WiFS radiometer, sensitivity to linear polarization must 
be less than 2%, and the actual degree of polarization sen- 
sitivity must be characterized for each channel. 

A generalized protocol for characterizing the polariza- 
tion sensitivity of a radiometer is given here. The in- 
strument should view a source of linearly polarized ra- 
diance, and its apparent radiance response Li(A) should 
be recorded. The instrument should then be rotated 90° 
about its FOV axis, still viewing the linearly polarized ra- 
diance source, and the apparent radiance response L 2 (A) 
should be recorded. The polarization sensitivity of the in- 
strument will be calculated as 


extremes of instrument operating temperatures, then po- 
larization sensitivity measurements should be made at sev- 
eral additional temperatures in that range. 

4.4 CALIBRATION OF IOP METERS 

Calibration of beam transmissometers has traditionally 
been carried out by means of in-air measurements, with a 
subsequent adjustment for changes in the Fresnel reflec- 
tions by the windows upon submergence into water. The 
660 nm transmissometers produced by Sea Tech, Inc. are 
independently calibrated at the factory after which cali- 
bration is maintained via frequent air calibrations (Sec- 
tions 5.4.1 and 6.7). This approach is adequate for this 
generation of 660 nm transmissometers. 

The newer generation of combined spectral absorption 
and beam transmission meters require a more accurate cal- 
ibration procedure. This class of instruments must be cal- 
ibrated with optically pure water, which has been freshly 
prepared using multiple-pass reverse osmosis filtration and 
analyzed with special optical techniques to verify optical 
purity. Ordinary laboratory grade pure water, prepared 
using distillation, deionization, or ordinary reverse osmo- 
sis filtration, will not ordinarily be optically pure enough 


to serve as a calibration standard. In general, it will not 
(15) be practical to perform this calibration procedure either 
at sea or in the average oceanographic laboratory. An- 
nual recalibration, either by the manufacturer, or another 
As required for the SeaWiFS radiometer, airborne and qua pfi ec ] laboratory, together with air calibrations every 
shipboard radiometers must satisfy P(A) < 0.02. three days during deployments, are recommended (Sec- 

A very simple, semi- quantitative test of a radiome- tion 5 2 . 4 ). The air calibration procedure for this class 
ter’s polarization sensitivity can be performed outdoors of i Ilstrumen ts is also more elaborate than for the red 
on a cloud- and haze-free day. The instrument should (gfiOnm) transmissiometers (Section 5.2.4). 
be pointed at the sky in the zenith-sun plane at an an- Procedures for calibrating a scattering meter depend on 
gle of approximately 90° from the sun, and its response t ^ e bailed design of the instrument. For example, hack- 
ly A) recorded. Since singly- scattered Rayleigh radiance scat t er ing meters are calibrated with Spectralon plaques 
is 100% polarized at a scattering angle of 90°, if aerosol of k nown reflectance. Detailed protocols for calibrating 
scattering is small, the sky radiance viewed at this angle sca ttering meters remain to be developed in a future revi- 
will be strongly polarized. If the instrument is then rotated s j on j- 0 document. 

90° about its FOV axis to measure L 2 (A), an approximate 


_ 2(L 1 (A)-L 2 (A)) 

1 j ” Li(A) + L 2 (A) 


estimate of P( A) may be computed, as above. 

Specification of detailed protocols for laboratory char- 
acterization of a radiometer’s polarization sensitivity will 
require more attention than is available here. In particular, 
protocols should be developed which describe in detail: 

1) laboratory setups for producing a stable, uni- 
form, extended source of linearly polarized radi- 
ance; and 

2) laboratory procedures for measuring the actual 
degree of polarization of the polarized radiance 
source and for determining the uncertainty of 
the polarization sensitivity estimate achieved us- 
ing a particular experimental setup. 

Temperature dependence of an airborne radiometer’s 
polarization sensitivity should initially be characterized at 
5° and 30° C. If significant differences in P(A) exist at these 


4.5 SUN PHOTOMETERS 

Sun photometer calibrations should be performed at 
least annually, when used consistently, through a Langley 
calibration procedure. In this procedure, the solar sig- 
nal transmitted through the atmosphere is measured over 
different air masses, i.e., at different solar zenith angles, 
throughout the course of the day. The validity of using 
these measurements as calibration of a sun photometer 
hinges strongly on the assumption that aerosols are uni- 
formly distributed and do not vary throughout the day. 
These Langley calibrations, therefore, should be performed 
in areas of atmospheric stability with low aerosol loading. 
Suitable locations include the astronomical observatories 
at Mauna Loa, Hawaii, and Kitt Peak, Arizona. In be- 
tween these calibrations, radiance calibrations with stan- 
dard lamps may be used as a stability check (Shaw 1976). 
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Temperature stability should be characterized for each 
instrument. Linearity and spectral calibrations should be 
performed with the same frequency as the absolute cali- 
bration; these characterizations must be performed in the 
laboratory. 

4.6 SKY RADIANCE CAMERAS 

Absolute and spectral calibrations should be performed 
on the radiance distribution camera before and after each 
cruise. A full characterization of the instrument should 
be performed initially, including camera lens roll-off char- 
acteristics for each camera (Voss and Zibordi 1989). If 
attenuation devices are used to prevent solar saturation, 
these should be calibrated frequently to track drift. Lin- 
earity calibrations should also be performed with the same 
frequency as the absolute and spectral calibration. Pro- 
cedures for characterizing this class of instruments are es- 
sentially the same as for other radiance detector systems. 
Each individual detector element in the detector array is 
essentially regarded as an independent radiometer. 

4.7 PIGMENT CALIBRATIONS 

HPLC equipment used to measure phytoplankton pig- 
ment concentrations is to be calibrated using standards 
distributed under the auspices of the SeaWiFS Program.! 
Concentrations of the pigment standards should be de- 
termined using a monochromator- based spectrophotome- 
ter and the extinction coefficients summarized in Bidigare 
(1991). The use of a diode array spectrophotometer can 


result in the underestimation of chlorophyll a concentra- 
tion because of errors associated with fluorescence con- 
tamination (Latasa et al. 1994). Bench fluorometers used 
to measure concentrations of extracted chlorophyll a and 
phaeopigments should be calibrated using authenticated 
standard chlorophyll a adopted for HPLC. In situ fluorom- 
eters should be calibrated against extracted chlorophyll a 
from concurrent bottle samples. 

4.8 CTD CALIBRATIONS 

The conductivity probe, temperature probe, and pres- 
sure transducer of the CTD should be recalibrated be- 
fore and after each major cruise by a properly equipped 
physical oceanographic laboratory, including those main- 
tained by many CTD manufacturers. In addition, the con- 
ductivity probe should be independently calibrated dur- 
ing the course of each cruise by obtaining salinity water 
samples (Section 5.2.3) simultaneous with CTD readings. 
These salinity samples are to be analyzed, either at sea 
or ashore, with a laboratory salinometer calibrated with 
International Association for the Physical Sciences of the 
Ocean (IAPSO) Standard Seawater. 

If simultaneous deployment of the CTD with optical 
instruments having independent pressure transducers is 
practical, the two depths measured by the different instru- 
ments should be compared over the range of the cast. If 
depth measurements disagree significantly, these compar- 
isons may be used to correct whichever transducer is found 
to be in error through analysis of pre- and post-cruise pres- 
sure transducer calibrations. 


t The HPLC pigment standards are available from Robert R. 
Bidigare, Department of Oceanography, University of Hawaii 
Honolulu, HI 96822. 
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Chapter 5 


Measurement Protocols 

Introduction 

This chapter describes the methods and procedures by which measurements are to be made during oceanographic 
deployments for SeaWiFS algorithm development and validation. Section 5.1 gives protocols for obtaining ra- 
diance and irradiance measurements from ships and aircraft. Section 5.2 gives protocols for measuring IOP, 
chlorophyll a fluorescence, and CTD profiles. Section 5.3 gives protocols for sun photometry and sky radiance 
distribution measurements. Section 5.4 covers protocols for obtaining and processing water samples for phyto- 
plankton pigment concentration measurements and spectrophotometric determination of absorption by particles 
and dissolved matter. Sections 5.5-5. 7 describe ancillary observations to be recorded at each station, radiometric 
measurements from moorings, and radiometry on drifting buoys. Finally, Section 5.8 gives protocols addressing 
sampling strategy, station location, and coordination of ship and aircraft measurements in Case-1 and Case-2 
water masses. 


5.1 IRRADIANCE AND RADIANCE 

Determinations of in-water spectral Ed, E u , and L u , 
both near the surface and as vertical profiles, are required 
for calibration and validation of the water-leaving radi- 
ance as retrieved from the SeaWiFS satellite sensor. Near- 
surface measurements should profile through at least the 
top three optical depths to reliably extrapolate to z = 0; 
it is essential to obtain a profile through at least the top 
optical depth. To better characterize the water column 
for remote sensing applications, e.g., primary productiv- 
ity estimation, deeper vertical profiles should be made to 
200 m, or seven optical depths whenever possible. Sea bed 
reflection influences on L u and E u should be avoided for 
SeaWiFS validation and algorithm development by collect- 
ing data only from water deeper than six optical depths 
for Ed (490); remote sensing applications for optically shal- 
low situations where bottom reflectance is present are not 
within the scope of these protocols. 

There are two primary sources of error in the deter- 
mination of these optical parameters: the perturbation of 
the in-water radiant energy field by the ship (Gordon 1985, 
Smith and Baker 1986, Voss et al. 1986, and Helliwell et al. 
1990), and atmospherically induced variability in radiant 
energy incident on the sea surface during in-water mea- 
surements (Smith and Baker 1984). The influence of ship 
shadow on the vertical profiles of Ed, E u , and L u is depen- 
dent upon the following variables: solar zenith angle, the 
spectral attenuation properties of the water column, cloud 
cover, ship size (length, beam, draft, and freeboard) and 
color, and the geometry of instrument deployment. At- 
mospheric variability is primarily dependent upon sun ele- 
vation and variations in cloud cover. The near surface in- 
water data also show variability caused by wave focusing, 


which can be minimized at a fixed depth by averaging over 
several wave periods, but which can pose severe problems 
in vertical profiles during which the instrument descends at 
speeds of 0.5-1 ms -1 . Raman scattering and fluorescence 
result in second-order errors near 490 nm (CDOM fluores- 
cence), and at longer wavelengths, contributions from phy- 
coerythrin and chlorophyll a fluorescence and water Raman 
scattering are significant. 

5.1.1 Ship Shadow Avoidance 

The complete avoidance of ship shadow, or reflectance, 
perturbations is a mandatory requirement for all radio- 
metric measurements to be incorporated into the SeaWiFS 
validation and algorithm database. The influence of ship 
shadow is best characterized in terms of attenuation length 
l/Ad(A) (Gordon 1985). Because L\y is required with an 
uncertainty of 5% or better, the protocol requires that ver- 
tical profiles be measured outside the effects of ship per- 
turbation to the radiant energy field. To accomplish this, 
the instrument must be deployed from the stern, with the 
sun’s relative bearing aft of the beam. 

Estimates of the minimum distance away from the ship, 
under conditions of clear sunny skies, are given below. The 
distances are expressed in attenuation lengths to minimize 
error. For Ed{ A) measurements, the general equation for 
distance away, f , in meters is given as 

sin (48 ^2 

* K d { A) ' k 

The distance from the ship is required to be 3/X u (A) m for 
E u {\) and 1.5/E L (A)m for L u { A) measurements. These 
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distances should be increased if the instrument is deployed 
off the beam of a large vessel. 

A variety of methods have been used to deploy opti- 
cal instruments beyond the influence of the ship. During 
CZCS algorithm development, floating plastic frames were 
equipped with small winches and instruments to obtain 
near surface optical profiles at some distance away from 
the ship. An umbilical cable provided power and data 
transfer. These platforms, while being somewhat difficult 
to deploy, worked well at avoiding ship shadow. Alterna- 
tively, extended booms can be used to deploy the instru- 
ment away from the ship and have the advantages of allow- 
ing relatively rapid deployment and simultaneous rosette 
bottle sampling. As a point of caution, however, very long 
booms may accentuate unwanted vertical motions due to 
ship pitch and roll. 

Waters et al. (1990) used an optical free-fall instrument 
(OFFI) that allows optical data to be obtained outside 
the influence of ship perturbation. In addition, the OFFI 
approach allows optical data to be obtained independently 
from violent ship motion, which may be transmitted to the 
instrument via the hydrowire, especially on a long boom. 
Yet another method for the deployment of optical sensors 
is via an ROV. Some groups, e.g., Smith (pers. comm.), 
have deployed a spectrometer on an ROV and obtained 
data completely free of ship influences. 

The above criteria for ship shadow avoidance are ad- 
mittedly very conservative. Unfortunately, the above cited 
models and observations provide only approximate guid- 
ance on minimum distances at which ship reflectance and 
shadow effects become insignificant under all circumstan- 
ces. Therefore, the SPSWG has adopted relatively extreme 
distance criteria, recognizing that in many specific combi- 
nations of lighting conditions, ships and optical properties, 
ship shadow, and reflection effects may become unimpor- 
tant much closer to the ship. 

The essential requirement is that each investigator es- 
tablish that any measurements of E d , E u , and L u , submit- 
ted for SeaWiFS validation and algorithm development, 
are free from ship-induced errors. The simplest way to do 
this is to adhere to the above distance criterion, which is 
not difficult when using either a tethered free-fall system 
or instruments mounted on an ROV. In other cases, it is in- 
cumbent on the investigator to otherwise demonstrate the 
absence of ship effects, e.g., through analysis of a series of 
profiles at increasing distance. 

At wavelengths where attenuation lengths are the same 
order of magnitude as, or less than, the size of the instru- 
ment package, e.g., in the ultraviolet-B (UVB) or red and 
near infrared spectral regions, care must be taken to con- 
sider possible perturbation of the radiant energy field by 
the instrument package itself. Methods of accounting for 
self-shadowing by the instrument are not well established 
and new measurement approaches must be developed (Sec- 
tion 5.1.6). 


5.1.2 Depth Resolution in Profiles 

The instrument sampling rate and the speed at which 
the instrument is lowered or raised through the water col- 
umn should yield at least two, and preferably six to eight, 
samples per meter. 

5.1.3 Instrument Dark Readings 

The dark current of optical sensors is frequently tem- 
perature dependent. As a consequence, accurate radiomet- 
ric measurements require that careful attention be given to 
dark current variability. It is recommended that each opti- 
cal measurement be accompanied by a measurement of the 
instrument dark current. When there is a large tempera- 
ture difference between the instrument on the deck and 
the water temperature, the instrument should be allowed 
to equilibrate with ambient water temperature at the be- 
ginning of each cast. 

Deep casts, e.g., 500 m, may permit the determination 
of the dark current in each optical channel at the bottom of 
each cast. Many instruments, however, are not designed to 
be safely lowered to 500 m, and this approach is usually not 
feasible. Furthermore, there is some intrinsic uncertainty 
over possible contamination by bioluminescence when dark 
readings are obtained in this way. If the instrument is 
equipped with a shutter, dark currents can be measured at 
any depth in the cast. If the dark current is not determined 
during the cast, it should be determined as soon as possible 
after the instrument is returned to the deck. 

Temperature effects on sensor responsivity can be sig- 
nificant and should not be ignored. Therefore, sensors 
should be equipped with thermistors on detector mounting 
surfaces to monitor temperatures for data correction. Oth- 
erwise, deck storage should be under thermally protected 
conditions prior to deployment and on-deck determination 
of dark voltages. 

5.1.4 Surface Incident Irradiance 

Atmospheric variability, especially under cloud cover, 
leads directly to variability of the in-water light field and 
must be corrected to obtain accurate estimations of optical 
properties from irradiance or radiance profiles. First order 
corrections for this variability can be made using above 
water (on deck) measurements of downwelling spectral ir- 
radiance, E s ( A) = E d ( 0 + , A). Smith and Baker (1984) and 
Baker and Smith (1990) theoretically computed the irra- 
diance just below the air- water interface, E d { 0“, A), from 
deck measurements to correct in-water profile data. 

The deck sensor must be properly gimballed to avoid 
large errors in E S (X) due to ship motion in a seaway. Im- 
proper gimballing can actually accentuate sensor motion 
under some circumstances, however, and this aspect of 
a shipboard radiometer system must be engineered with 
some care. 
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Recently, Waters et al. (1990) have demonstrated a 
method to more directly determine Ed( 0~,A) by deploy- 
ing an optical surface floating instrument (OSFI) to obtain 
continuous optical data just below the air-water interface. 
These Ed( 0“ , A) are used as a normalization factor to cor- 
rect for variations in irradiance during a vertical profile, 
or over the period of a day for a sequence of profiles. Re- 
search and development is needed to determine whether 
this should be the preferred approach for SeaWiFS valida- 
tion measurements. 

5.1.5 Instrument Attitude 

An instrument’s attitude with respect to the vertical is 
a critical factor in measurements of Ed{z , A) and E u (z , A), 
and is only slightly less critical for L u (z, A). Roll and pitch 
sensors must, therefore, be installed in the underwater ra- 
diometers used for the SeaWiFS program. The data from 
these attitude sensors are to be recorded concurrently with 
the data from the radiometric channels and are to be used 
as a data quality indicator. It is not deemed necessary to 
determine or control attitude determination errors result- 
ing from surface wave-induced accelerations at very shal- 
low depths. 

5.1.6 Instrument Self-Shading 

Gordon and Ding (1992) recently modeled the errors 
introduced by an instrument’s own shadow in direct mea- 
surements of L u ( A) and E u ( A). For this source of error to 
be less than 5%, without modeled corrections, the instru- 
ment radius r must satisfy r < [40a(A)] for E u ( A) and 

r < [l00a(A)] 1 for L u (\). They calculate for A = 865 nm 
in pure water, as an example, that the instrument radius 
must be approximately 0.3cm to measure E u (865) with a 
maximum of 5% error; the instrument radius must be sig- 
nificantly smaller for direct measurement error in L u (865) 
to be 5% or less. 

Gordon and Ding (1992) also propose a simple model 
for correcting L u ( A) and E U (X) for the self-shadowing ef- 
fect. They write 


MA) 


L U (X) 

l-e(A) 


(17) 


with instruments having maximum diameters of 24 cm for 
A < 650 nm, and with instruments of maximum diame- 
ter 10 cm for 650 < A < 700 nm at solar zenith angles 
(0 O ) greater than or equal to 20°, and maximum chloro- 
phyll concentrations of lOrngm -3 . To measure L u ( A) cor- 
rectable to less than 5% error at 6$ — 10° (with chlorophyll 
concentrations less than or equal to lOmgm -3 ), maximum 
instrument diameters are 12 cm for A < 650 nm and 5 cm 
for 650 < A < 700 nm. Even with these corrections, how- 
ever, instrument diameters of 1 cm or less must be used to 
assure self-shading L U (A) errors are 5% or less at 780 and 
875 nm. 

The Gordon and Ding (1992) model predictions have 
recently been compared to experimental measurements of 
L U (A) just beneath the sea surface, using a fiber-optic ra- 
diometric probe (Zibordi and Ferrari 1994). The exper- 
iment was performed in a lake, with solar zenith angles 
ranging between 25° and 50°, on several days with cloud- 
free skies. Spectrophotometric methods (similar to those 
in Sections 5.4.2 and 6.8.2) were used to measure absorp- 
tion by particles and Gelbstoff. At wavelengths of 500, 
600, and 640 nm, a series of discs was employed to vary 
instrument self-shading geometry in several steps over the 
range 0.001 < a(A)r <0.1. The Gordon and Ding (1992) 
model predicted self-shading radiance and irradiance ef- 
fects that may be applied as corrections, and which agreed 
with measured values within 5% and 3% respectively. The 
model corrections were all biased high relative to the mea- 
sured values. Zibordi and Ferrari (1994) chose to compare 
their measurements to the Gordon and Ding (1992) point 
sensor model, and use of their finite-sensor model results 
may have improved the comparisons. 

This initial confirmat ion of the Gordon and Ding (1992) 
instrument self-shading model is confined to clear-sky con- 
ditions, solar zenith angles greater than 25°, near-surface 
L u (0~,A) and £ U (0“,A), and a(A)r less than or equal to 
0.1. Additional theoretical and experimental research will 
be necessary to generalize this correction for cloudy sky 
conditions and for variations with depth in L u (z,A) and 
£ u (z,A) profiles. The above restrictions notwithstanding, 
the excellent agreement shown so far covers a very im- 
portant range of conditions for SeaWiFS algorithm devel- 
opment and validation. Therefore, a provisional proto- 
col for applying self-shading corrections to L u (0',A) and 
£u(b - ?^) is given in Section 6.1.7. 


and 

e(A) = 1 - e~ k ' a{x)r , (18) 

where ~ is the true value, ~ is the measured value, k' — 
y / tan 6q w * @ow is the refracted solar zenith angle and y is 
an empirical factor for which they give values determined 
by fitting their model results (y « 2). A similar correction, 
with a different table of values for y applies to E u ( A). 

When the above geometric corrections are applied, Gor- 
don and Ding (1992) estimate that errors less than or equal 
to 5% in L U (X) could be determined from measurements 


5.1.7 Airborne Ocean Color Radiometry 

Ocean color may be measured from aircraft using ei- 
ther imaging radiometers (usually flown at high altitude), 
or single FOV spectral radiometers (usually flown at low 
altitude to measure profiles of ocean color beneath an air- 
craft’s trackline). For SeaWiFS validation, either of these 
types of data can be an extremely valuable addition to 
shipboard measurements. 

Aircraft measurements observe the horizontal variabil- 
ity in ocean color radiance spectra on spatial scales which 


33 



Ocean Optics Protocols for SeaWiFS Validation, Revision 1 


are much smaller than SeaWiFS pixels; therefore, these 
data are more comparable to shipboard measurements. At 
a qualitative level, this information can indicate how well 
shipboard radiometric and bio-optical measurements can 
be compared to SeaWiFS data at greater than 1 km pixel 
resolution. For more quantitative work, if an airborne ra- 
diometer’s characteristics conform to specifications of Sec- 
tion 3.3, and if accurate corrections are applied for at- 
mospheric and surface reflection effects, the high resolu- 
tion airborne data may be spatially integrated to quan- 
titatively extrapolate shipboard bio-optical measurements 
to SeaWiFS pixel scales. (Note that the original data must 
be corrected for sun and sky glitter, as well as atmospheric 
transmission and path radiance, before spatial integration 
to SeaWiFS resolution.) Finally, in cases when airborne 
normalized water- leaving radiances estimated from aircraft 
measurements are independently validated as being in con- 
formity with SeaWiFS uncertainty goals, they may be used 
for direct radiometric validation of SeaWiFS radiometric 
performance. 

Methods for atmospheric correction and estimation of 
normalized water-leaving radiances from high altitude air- 
borne ocean color imagery are nearly identical to, and as 
challenging as, those methods which must be applied to 
SeaWiFS data itself (Carder et al. 1992 and Hamilton et 
al. 1992). These problems and their solutions lie beyond 
the scope of the ocean optical protocols per se, at least in 
this revision. The elements of ocean color radiometry from 
low flying aircraft are more closely associated with ship- 
board ocean optical measurements and will be addressed 
here. 

Under clear skies, or partial cloud cover, the radiometer 
should be pointed approximately 15-20° from nadir in the 
azimuth direction opposite from the sun, to avoid specular 
sun glint from wave facets. If seas are rough and the solar 
zenith angle is small, an even larger nadir angle may be 
necessary. Under completely uniform overcast conditions, 
the instrument should point to the nadir. If an imaging 
radiometer is fixed in the nadir plane, sampling should be 
scheduled to provide a solar zenith angle between 35° and 
45°, and all flight lines should be flown toward or away 
from the sun. 

For quantitative low altitude ocean color radiometry, 
solar downweiling spectral irradiance must be measured 
using a sensor mounted on top of the aircraft. Incident so- 
lar irradiance spectra are modified significantly by trans- 
mission through the atmosphere above the aircraft, and di- 
rect spectral irradiance measurements are absolutely essen- 
tial for computing normalized water-leaving radiances from 
flight level radiances. The measured flight- level down- 
welled irradiances should be used to model atmospheric 
transmission and path radiance corrections, which must 
be applied at altitudes of 100 m and higher. Reflected 
skylight must also be removed from the measured radi- 
ances. Finally, for use in algorithms related to SeaWiFS, 


the water- leaving radiances must be normalized by the ra- 
tio Fq(X)/E s (X) j where Fq(X) is the mean extraterrestrial 
solar flux, and E s ( A) is the measured incident spectral irra- 
diance propagated to the sea surface. This normalization 
is also critical for ratios of water-leaving radiance, because 
the spectral quality of incident irradiance is modified sig- 
nificantly, in a generally unpredictable way, by atmospheric 
transmittance. 

Under either clear skies or a completely uniform cloud 
cover, i.e., overcast skies, the downweiling irradiance mea- 
surements may be used directly to model corrections for 
low altitude atmospheric transmission and path radiance, 
and for reflected sky radiance. On the other hand, ocean 
color radiance data measured under variable cloud con- 
ditions are very difficult to correct and interpret. Under 
these conditions, the area on the sea surface viewed by the 
radiometer is often illuminated differently than is the in- 
cident irradiance sensor at the aircraft, and sky radiance 
incident on the surface also varies unpredictably. 

The requirement to directly measure incident spectral 
irradiance may sometimes be relaxed when the experimen- 
tal objective is simply to characterize spatial variability of 
bio-optical properties qualitatively. Under cloud- free con- 
ditions, moreover, it may often be possible to model surface 
incident irradiance accurately enough for semi-quantitative 
analyses of chlorophyll concentrations from low altitude 
airborne spectral upwelled radiance measurements alone. 

For atmospheric correction and analyses of aircraft ra- 
diance data, the following ancillary parameters must be 
recorded: date and time [Greenwich Mean Time (GMT)], 
azimuth and nadir angle of the radiometer, heading of the 
aircraft (and all other aircraft attitude angles, if available), 
the flight altitude, and air pressure at that altitude. The 
height of the inversion should also be recorded, as most 
of the aerosols will normally be concentrated under this 
height. Furthermore, a precise navigation recording is re- 
quired for identification of the aircraft track within the 
satellite image and for comparisons to shipboard measure- 
ments. In addition to the radiometer, a video camera is 
very helpful to check for the effects of cloud shadows, foam, 
and ocean fronts — effects which are difficult to identify in 
the horizontal radiance profile alone. 

5,1.8 Shipboard Radiance Measurements 

Water- leaving radiances can, of course, be measured 
from the deck of a ship using methods directly analogous to 
those of low-altitude airborne radiometry (Section 5.1.7). 
Given a very close proximity to the sea surface, no correc- 
tions are needed for atmospheric transmission and path ra- 
diance. Shipboard water-leaving spectral radiances must, 
however, be corrected for light reflected at the sea surface, 
and must also be normalized using direct measurements of 
E s (A). Performance characteristics for radiometers used to 
measure E a ( X) and L u { A) are specified in Sections 3.2-3. 3. 
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The best position for measuring the water-leaving ra- 
diance is usually near the bow of the ship. At this po- 
sition, the ship will have the least influence with respect 
to shading or reflection. It is also easy to point in a di- 
rection away from the sun to reduce specular reflection of 
sunlight. While underway, ocean color radiance measure- 
ments should always be made from the bow, because this is 
the only position where the water is normally undisturbed 
by the ship’s wake and without foam. 

To measure L w ( A), the radiometer should point to the 
sea surface with an angle of about 20° from nadir and away 
from the sun’s azimuth by at least 90°. (A viewing angle 
that is 180° away from the sun’s azimuth may be contam- 
inated by the glory phenomenon.) This is the point with 
the minimum in specularly reflected sun and sky radiance. 

nder overcast conditions, the instrument should point 
closer to nadir, but not so close that ship shadow or hull 
reflections affect the water-leaving radiance. (Indeed, some 
arger ships have a long bulbous bow nose, which must be 
avoided ) Under a high sun elevation (such as in tropical 
regions) and a rough sea, it may be very difficult, or even 
impossible, to avoid large amounts of specularly reflected 
sunlight, since large nadir angles cause severe problems due 
to variations in reflected skylight. In addition, foam and 
floating material must be avoided during measurements. 

Because of temporal variability due to waves, it is im- 
portant to record a number of spectra within a period of 
a few seconds, e.g., 30 spectra within 15s. Before calcu- 
* tln f mean and standard deviation spectra, outliers 
should be removed by computing initial estimates of these 
statistics and rejecting radiance spectra containing values 
ailing more than 1.5 standard deviations from the esti- 
mated mean. 

Incident spectral irradiance E a ( A) should be measured 
at a location free from shadows of, and reflections by 
the ship’s superstructure. As in low-altitude aircraft ra- 
diometry, E s ( A) must be used both to normalize mea- 
sured water-leaving radiance and to calculate spectral re- 
flectance. Some investigators (Carder and Steward 1985 
and Doerffer pers. comm.) advocate calculating spectral 
reflectance using relative E 3 { X) estimated by viewing a ref- 
erence diffuse reflectance plate. Spectralon plates are used 
most often, but Kodak gray cards and plates painted with 
NextelTM coating have good Lambertian reflectance char- 
acteristics (Doerffer pers. comm.). 

The contribution by specularly reflected sky light to 
a water-leaving radiance measurement can be calculated 
from measurements of the sky radiance, made either by 
looking at a horizontal first surface mirror (a mirror with 
no layers other than the reflective surface) at the same 
nadir and azimuth angles, or by pointing the radiometer 
into the sky with a zenith angle equal to the nadir angle 
and with the same azimuth angle. The specular reflectance 
of the water surface can be calculated using the Fresnel 
equation (reflectance is approximately 0.02 within 20° of 
nadir). 


The sequence of water-leaving radiance, diffuser plate, 
and sky radiance [L sky (A)j measurements can usually only 
be completed under clear sky or uniform overcast condi- 
tions. When coupled with concurrent direct measurements 
of E s ( A), however, L w { A) and L sky ( A) can probably be 
measured under variable cloud conditions, provided a se- 
quence of several [L w (\),E a (\)\ and (L sky (A), £ S (A)], sets 
of spectra can be measured in rapid succession without 
significant variations. 

5.2 ANCILLARY PROFILES 

Beam transmittance, CTD profiles, and chlorophyll a 
fluorescence should be measured at the same stations as the 
irradiance and radiance measurements. Preferably, these 
auxiliary profiles should be measured simultaneously with 
or otherwise immediately before or after, the radiometer 
profiles. If possible, these profiles should be made in con- 
junction with bottle samples. For the verification of the 
satellite sensor, these data will be used as a guide to the 
uniformity of the first optical depth and to determine water 
bottle sampling depths. 

The IOP, fluorescence, and CTD profiles will also be 
used as a guide for, and constraints on, the smoothing 
of K(z, A) from the radiometric profiles. The location of 
maxima and other features in the structure of these profiles 
identify inflection points for segmenting the optical profiles 
into finite depth elements (layers) for the analytical meth- 
ods described by Mueller (1991) or Petzold (1988). Both of 
these techniques use multiple segments for the statistical 
ht of analytic functions to the measured profiles. These 
data will also be used to develop and validate pigment and 
primary productivity algorithms. 

5.2.1 Beam Transmittance 

The windows on the beam transmissometer must be 
cleaned with lens cleaner or a mild detergent solution and a 
soft cloth or tissue, rinsed with distilled water, then rinsed 
with isopropyl alcohol and wiped dry. An approximate air 
calibration reading should be made before every cast to 
verify that the windows are clean. A transmissometer dark 
voltage should also be measured at this time. These on- 
deck air calibrations are not, however, very reliable mea- 
sures of temporal drift or degradation in the instrument’s 
source or detector. In the humid, or even wet, environ- 
ment on the deck of a ship, the windows are often quickly 
obscured by condensation, and the glass also tends to ab- 
sorb enough water to affect transmission slightly (Zaneveld 
pers. comm.). A very careful air calibration should be per- 
formed before and after each cruise under dry laboratory 
conditions. During an extended cruise, it is also recom- 
mended to remove the instrument to a dry location in a 
shipboard laboratory, and after allowing several hours for 
the windows to dehydrate, a careful air calibration should 
be performed. Only the laboratory air calibrations should 
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be used in the final processing of beam transmissometer 

dat Both the laboratory condition air calibration and dark 
voltages, and the factory calibration voltages, assume t e 
data acquisition system measures instrument response as 
true volts. It is imperative, therefore, to calibrate the 
end-to-end analog-to-digital (A/D) data acquisition sys- 
tem and characterize its response to known input voltages. 
Corrections are in the form of a linear function. 

V = g(T)V + f(T), ( 19 ) 

where T is temperature, must usually be applied to ex- 
ternal voltage inputs recorded with the A/D circuits o 
CTDs or profiling radiometer systems. The range depen- 
dent A/D bias coefficients should be determined at ap- 
proximately 5° C intervals, over the range from 0 25 V, 
to characterize the temperature sensitivity of the data ac- 

quisdion^sy development ^ bi(>optical algor ithms describ- 
ing the inherent and apparent optical properties of the wa- 
ter and for algorithms estimating primary productivity, 
more stringent requirements are recommended for trans- 
missometer calibration and characteristics. s Pe ctral mea- 
surements of beam transmittance should be made with ab- 
solute uncertainties of 0.1% transmittance per meter, or 
0.001 m -1 beam attenuation coefficient c(A). 

Different and more stringent calibration protocols are 
specified in Section 5.2.4 for the new combined spectral 
absorption and beam attenuation meters. Laboratory ca - 
ibrations of these instruments are based on measurements 
using optically pure water, produced using extraordinarily 
careful reverse osmosis procedures. It is ordinarily imprac- 
tical to produce and verify optically pure water either at 
sea or in the average oceanographic laboratory. As with 
the older red-wavelength beam transmissometers, air cal- 
ibrations are prescribed as a basis for correcting instru- 
mental drift between calibrations at a qualified laboratory. 
The increased sensitivity required for the new instruments 
makes it essential to carry out more frequent air calibra- 
tions at sea. As discussed above, accurate air calibrations 
must be carried out in a dry laboratory environment with 
controlled relative humidity— conditions which may be dif- 
ficult to implement on some ships. 

5.2.2 Chlorophyll a Fluorescence 


inputs characterized. During processing, a correction o 
the form given in (19) must usually be applied to values 
recorded with the A/D circuits of CTDs and profiling ra- 
diometer systems. As in beam transmittance, the range 
dependent A/D bias coefficients should be determined at 
approximately 5° C intervals over the range from 0 25 C 
to characterize the temperature sensitivity o t e a a ac 

quisition system. , 

Zero fluorescence offsets should be measured on deck 
before and after each cast; the optical windows should be 
shaded to avoid contamination of the zero offset value y 
ambient light. Before each cast, the fluorometer windows 
should be cleaned following the manufacturer s instruc- 

tions. 

For chlorophyll a determinations, the fluorescence mea- 
surements should be compared to HPLC and extracted pig- 
ment measurements from discrete water samples, to permit 
comparison with JGOFS standard measurements and his- 
torical databases. The in situ fluorescence measurements 
will be used to provide continuous vertical profiles of in- 
terpolated pigment concentration using the bottle samp es 
as tie points. 


5.2.3 CTD Profiles 

Vertical profiles of CTD should be measured to at least 
the depth of the deepest bio-optical profile. If the station 
schedule will permit it, sections of CTD casts extending 
to 500 m, or deeper, will be useful for computing relative 
quasi-geostrophic currents and shear, which may affect the 
advection and mixing of bio-optical properties during a 
cruise. A real-time analysis and display of the CTD pro- 
file together with displays of c(660) and m situ fluores- 
cence profiles, should be available as a guide in choosing 
the depths at which water sampling bottles will be closed 
If possible, a few deep (1,500 m depth or greater) CT D 
and bottle sample profiles should be made during each 
cruise to obtain data for calibrating the CTD s conductiv- 
ity probe. During these CTD calibration casts , water sam- 
ples should be taken at depths where the vertical gradient 
of salinity is very small. This practice will minimize errors 
in the conductivity calibration resulting from the spatial 
separation of the water bottle and CTD profile. The bot- 
tled salinity samples may be stored for post-cruise analyses 
ashore at a laboratory equipped with an accurate salmome- 

ter and IAPSO Standard Seawater, if suitable equipment 

j A fiVio cViin 


An in situ fluorometer should be employed to mea- 
sure a continuous profile of chlorophyll a fluorescence. The 
fluorometer should be mounted on the same underwater 
package as the transmissometer, CTD, and water sampler, 
if one is employed. If possible, the radiometer should also 

be on this package. . . 

The A/D channel used to acquire and record signa 

voltages from the in situ fluorometer must be calibrated, 
and its temperature-dependent response to known voltage 


5.2.4 Spectral Absorption Profiles 

Currently, there are two methods for measuring spec- 
tral absorption coefficients which use commercially avail- 
able equipment. The first is conventional benchtop spec- 
trophotometer measurements of spectral absorption by par- 
tides extracted from water samples and concentrated on 
filters, followed by measurements of the spectral absorption 
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by the dissolved component in the filtrate (Section 5.4.2). 
More recently, in situ profile measurements of spectral ab- 
sorption have been made using instruments which capture 
most scattered light using a reflecting tube. A combination 
of these two methods allows one to measure profiles of to- 
tal absorption and absorption by the dissolved component 
with the reflecting tube instrument, and to independently 
partition the absorption spectrum at discrete depths into 
its components due to particulate and dissolved substances 
using filters and spectrophotometry. 

In situ spectral absorption coefficient profiles can also 
be measured with spectral radiometers conforming to the 
performance specifications listed in Section 3 above, if the 
radiometric package is extended to measure E d ( A) and 
E u (\), as well as scalar irradiances £ 0 d(A) and £ 0 u(A). 
This combination may be approached either using hemi- 
spherical collectors to measure upwelling and downwell- 
ing hemispherical irradiances, or by using cosine collec- 
tors on one radiometer in tandem with spherical collec- 
tors on another radiometer. Given these irradiance com- 
ponents, spectral absorption is then computed using Ger- 
shun’s equation as 

<2o) 

where E{ A) = E d {X)-E u {X) is vector irradiance, and K{ X) 
is the vertical attenuation coefficient for vector irradiance. 

In a recent experiment, comparisons between absorp- 
tion profiles measured using Gershun’s equation with E( A) 
and £ 0 (A) (scalar irradiance) data, and absorption profiles 
measured with a reflecting tube instrument, agreed within 
8% (Pegau et al. 1994). This level of agreement is well 
within the calibration uncertainties of the particular pro- 
totype instruments used for that experiment, which were 
approximately 10% uncertainties in both the scalar irra- 
diance radiometer and in the reflecting tube instrument. 
Less than 5% uncertainty in absorption is expected in fu- 
ture experiments. In very clear oligotrophic water, how- 
ever, uncertainty in water absorption values may make it 
impossible to realize this level of relative agreement. Un- 
fortunately, an adequate commercially available scalar, or 
hemispherical, irradiance radiometer does not exist today, 
and development and characterization of suitable spheri- 
cal or hemispherical collectors is a prerequisite to applica- 
tion of this method. The present version of the protocols 
will, therefore, emphasize the combination of filter spec- 
trophotometry (Section 5.4.2) and in situ reflecting tube 
total absorption measurements. The present section cov- 
ers field calibration and sampling procedures for using the 
reflective tube absorption and beam transmission meter. 
Recommended methods of data analysis to derive vertical 
profiles of a(z , A) and c(z, A) are presented in Section 6.8.1. 

It is always best to determine optical properties in situ , 
if possible. Sampling variability, changes of light inten- 
sities, filtration procedures, and sample degradation over 


time all affect the particulate matter and distort its true 
optical properties as they existed in the ocean, and as 
they determine the remote sensing reflectance viewed by 
SeaWiFS. The reflecting tube method has been used to 
measure spectral absorption in the laboratory for many 
decades (James and Birge 1938). In recent years, this 
method has been adapted for use in the ocean (Zaneveld 
1992). Suitable instruments are now commercially avail- 
able and are coming into general use within the oceano- 
graphic community. 

The reflecting tube does not perfectly gather all scat- 
tered light and transmit it to the detector, and as a result, 
there is a scattering error on the order of 13% of the scat- 
tering coefficient. This error can be largely corrected if 
the beam attenuation coefficient is measured simultane- 
ously. In that case, the scattering coefficient is obtained as 
6(A) = c(A)-a(A). By assuming that the measured absorp- 
tion is due to water and scattering error at a wavelength in 
the infrared, and by subsequent correction at other wave- 
lengths using a provisional 6(A), it is possible to correct the 
spectral absorption to within a few percent of the scatter- 
ing coefficient. Only in waters with very high scattering 
and very low absorption would this error pose a serious 
absorption error. This correction method is described in 
more detail in Section 5.4.5 below. 

The in situ reflecting tube is normally used to mea- 
sure profiles of total spectral absorption, but it has been 
recently demonstrated that if the unit’s intake is fitted 
with a large area 0.2 //m filter, the spectral absorption 
of the dissolved component can be measured (Zaneveld 
pers. comm.). A pair of reflecting tube absorption meters 
can thus be used to determine the separate constituents of 
absorption due to particulate and dissolved substances — 
a distinction of fundamental importance in relating ab- 
sorption to remote sensing reflectance. More traditionally, 
the filtration and spectrophotometry techniques developed 
over the last decade also lend themselves well to this task. 
Using the methods described in Section 5.4.2, the spec- 
tral absorption coefficient is partitioned into components 
associated with Gelbstoff, pigments, and non-pigmented 
particles (the latter sometimes referred to misleadingly as 
detritus). 

There are two primary weaknesses of the filtration spec- 
troscopic method: 

1. Uncertainties associated with correcting the fil- 
ter absorption values for scatter within the fil- 
ter are not well characterized. The correction 
factor for this phenomenon, called the (3 factor, 
does not currently have a unanimous consensus 
value. 

2. The instrumental uncertainty in absorption due 
to the dissolved component, which attenuates 
weakly over the 10 cm path of a typical spec- 
trophotometer, is also not well characterized. 
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Note that some longer path spectrophotometers are in use 
experimentally, and may reduce this source of uncertainty 
following a suitable validation procedure. In any event, 
when the filter spectroscopy absorption values are com- 
bined with total absorption values measured with the re- 
flecting tube, the absolute values may be largely corrected 
to remove these uncertainties, if two reflecting tube meters 
are not available. Conversely, these comparisons serve to 
validate and control the quality of the a(A) profiles mea- 
sured in situ . 

5,2,5 Backscattering Profiles 

While the spectral absorption can now be measured 
with sufficient accuracy for the purpose of SeaWiFS, back- 
scattering is not routinely measured, yet it has a large 
influence on the reflectance. The now generally accepted 
dependence of the remote sensing reflectance on the ratio 
of the backscattering coefficient and the absorption coef- 
ficient has never led to a systematic investigation of the 
dependence of 65(A) on the various particulate parameters. 
Gordon (1989a) showed that the irradiance reflectance de- 
pends on the shape of the scattering function in the back- 
ward direction. Zaneveld (1994) derived the functional de- 
pendence of the remote sensing reflectance on backscatter- 
ing. 

There is little historic data on the variation of the shape 
of the volume scattering function, 0(0, A), in the backward 
direction. Petzold (1972) described /3(A) measured with 
the GASM. This reference is the one most widely used to 
describe shapes of 0(6, A). Since that time, only Balch et 
al. (1994) have published new data regarding the shape of 
0(6, A). Several single wavelength integrated backscatter- 
ing 65(0, A) meters have been constructed (Maffione et al. 
1991 and Smart 1992), all of which measure a weighted in- 
tegral of the scattering function. Maffione’s meter has the 
centroid angle near 150°, whereas Smart’s meter has the 
centroid angle near 170°. The 150° centroid angle is more 
appropriate for SeaWiFS. This observation also points out 
a problem in the reporting of results from scattering me- 
ters, in that results are commonly reported as simple 65(A). 
In the future, more detail should be provided to allow in- 
tercomparisons between results reported by different inves- 
tigators. 

All scattering sensors measure a parameter that is a 
weighted integral of the volume scattering function. The 
weighting function W(0) depends on the geometry of the 
device. Knowledge of the weighting function is critical in 
comparing scattering meters. For each scattering meter, 
documentation should be available describing this weight- 
ing function and how it was derived. The weighting func- 
tion can generally best be derived by moving a Spectralon 
target through the scattering volume space. 

The scattering meter measurement may be defined as 

7T 

0 b (6,A6,\,A\) = 2tt Jp{e,X)W {6) sin Ode, (21) 
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where 0 is the centroid angle computed as 
f 6W (9) sin 0 d6 

0 = , (22) 

f W (6) sin 6 dO 
o 

and A 6 is the FWHM bandwidth of the weighting function 
multiplied by sin(0), A is the central wavelength of the color 
filter, and AA is the FWHM bandpass of the filter. While 
the parameters do not give a complete description of the 
weighting function, they allow an informed comparison of 
scattering meter results. Investigators making backscatter- 
ing measurements are encouraged to report backscattering 
meter results in the format 0^(6, AO, A, AA). 

Backscattering measurement should be made in close 
proximity to absorption meters, if both are used, to avoid 
space-time aliasing. 

5.3 ATMOSPHERIC RADIOMETRY 

This section is concerned with two types of atmos- 
pheric radiometric measurements that are distinct from 
those described earlier in Section 5.1. The first of these 
is the photometric measurement of the direct solar beam 
to determine the optical thickness of the atmosphere (Sec- 
tion 5.3.1). The second is a measurement of the sky ra- 
diance distribution using a radiance distribution camera 
(Section 5.3.2). 

5.3.1 Sun Photometry 

Measurements of the direct solar beam, using the sun 
photometer, should be performed during the optical sta- 
tions. If sky radiance distribution measurements are per- 
formed, it is important that these measurements are per- 
formed concurrently. While the preferred method of deter- 
mining the optical thickness of the atmosphere is by mea- 
suring the solar transmission as a function of solar zenith 
angle, the atmospheric conditions are rarely stable enough 
at sea for this method to work. Thus, a stable, well cal- 
ibrated photometer can be used with measurements at a 
single zenith angle to obtain the solar transmission, and 
thus, the aerosol optical depth. 

Atmospheric measurements should be performed only 
when clouds, including high cirrus, do not obstruct the 
solar disc. Careful documentation of sky conditions is im- 
portant, as are accurate recordings of the time of day and 
the geographic location of the station. The latter data 
are important in determining the true solar zenith angle 
and, hence, the air mass in the solar path. It should also 
be obvious that care should be taken to avoid ship per- 
turbations (e.g., stack gas) from interfering with the mea- 
surements. Ancillary measurements, such as barometric 
pressure, are important in separating Rayleigh scattering 
from the aerosol scattering. 
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5.3.2 Sky Radiance Distribution 5.4 WATER SAMPLES 


Complete sky radiance distributions should be mea- Duplicate samples should be taken at each of 12 depths, 

sured with a radiance distribution camera during the Sea- including at least three depths within the first attenuation 
WiFS radiometric initialization and validation optical sta- length, 1/AT(490); however, in coastal areas with short 
tions. For this purpose, it is critically important that these attenuation lengths, this may not be possible. Samples 
measurements be obtained whenever totally clear sky con- should also be taken in the in vivo fluorescence and beam 
ditions persist. Coincident with these measurements, sun attenuation maxima. The remaining samples should be 
photometer measurements should be obtained. When lo- spaced throughout the water column using beam attenua- 
cating the camera system for these measurements, it is tion, m situ fluorescence, and CTD profiles as a guide, 
important that the FOV be as unobstructed as possible. 

While it would be optimum to have a completely unob- 544 Pigment Analysis 
structed FOV, this is often not practical. During measure- 

ments, therefore, at least one hemisphere (defined by the Water samples should be taken at the site of, and simu - 
sun-zenith plane) should be unobstructed; through sym- taneously with, the surface in-water upwelled radiance and 
metry, this should yield a complete radiance distribution, reflectance measurements, and at depth increments sufn- 
Ship perturbations, e.g., stack gas, must be avoided dent to resolve variability within the top optical depth, 
completely. It is important to document where the instru- The K(z, A) profiles over this layer will be used to com- 
ment is located and what possible perturbations might ex- pu te optically weighted, near-surface pigment concentra- 
ist even though these effects may be obvious in the data, tion for bio-optical algorithm development (Gordon and 


It would also be highly desirable to add sky radiance 
measurements to every SeaWiFS algorithm development 
and validation cruise. Gordon (1989b) rigorously demon- 
strated the importance of determining JL d , the mean co- 
sine for downwelling radiance (Morel and Smith 1982), for 
the bio-optical interpretation of K d (z, A). Gordon (1989b) 
also showed, for cloud-free skies, how to obtain a reason- 
able estimate of /J d (0+) from spectral irradiance deck cell 
measurements with, and without, the sun blocked from the 
view of the irradiance collector. This procedure should be 
done routinely whenever it is practical to do so. Unfortu- 
nately, the collective scientific experience is that cloud-free 
skies rarely occur at ocean optical stations. 

If a calibrated spectral radiance distribution camera 
is available, then L sky (A,6>,0) should be measured several 
times during each spectral radiometer cast and used to 
compute Jj^(0 + , A). Radiance distribution cameras are ex- 
pensive to build, however, and one is not likely to be avail- 
able aboard every SeaWiFS validation vessel. A recom- 
mended alternative approach is to acquire: 

1) all-sky photographs taken either with a conven- 
tional camera or, preferably, a digitally recorded 
camera system (some research must be done to 
develop procedures for attenuating, or blocking, 
the sun’s image and for using filters); 

2) measurements, with a narrow FOV spectral ra- 
diometer of L s ky(A, Qi, <t>i) at several discrete an- 
gles Oi and <t>i\ and 

3) E a { A) measurements of ^sun (A) 4- £sky(A) with 
the sun’s image blocked by shadowing the deck 
cell irradiance collector (Gordon 1989b). 

During the prelaunch experiments, Md(0 + ,A) values esti- 
mated from measurements of these types should be com- 
pared with Ji d {0 + ,X) values determined from direct mea- 
surements of sky radiance using a spectral radiance distri- 
bution camera system. 


Clark 1980). 

When possible, samples should be acquired at several 
depths distributed throughout the upper 200 m of the wa- 
ter column [or in turbid water, up to seven optical depths, 
ln(( E( 0) /E{z)) = 7], to provide a basis for relating chloro- 
phyll a fluorescence signals to pigment mass concentration. 

For low uncertainty determinations of chlorophylls a, 
b } and c, as well as carotenoid pigments, HPLC tech- 
niques are recommended. It should be noted, however, 
that the reverse- phase C18 HPLC method recommended 
by the Scientific Committee on Oceanographic Research 
(SCOR) (Wright et al. 1991) is not capable of separat- 
ing monovinyl chlorophyll a from divinyl chlorophyll a nor 
monovinyl chlorophyll b from divinyl chlorophyll b. 1 his 
method, therefore, only provides estimates of total chloro- 
phyll a and total chlorophyll b concentrations, respectively. 

Divinyl chlorophyll a, the major photosynthetic pig- 
ment found in prochlorophytes, accounts for 10-60% of 
the total chlorophyll a in subtropical and tropical oceanic 
waters (Goericke and Repeta 1993, Letelier et al. 1993, 
Bidigare et al. 1994, and Bidigare and Ondrusek 1994). 
Divinyl chlorophyll a is spectrally different from normal 
(monovinyl) chlorophyll a, and its presence results in a sig- 
nificant overestimation of total chlorophyll a concentration 
as determined by the conventional HPLC methods (Goer- 
icke and Repeta 1993, Letelier et al. 1993, and Latasa et 
al. 1994). To avoid these errors, it is recommended that 
monovinyl and divinyl chlorophyll 0 be spectrally resolved 
or chromatographically separated in order to obtain an un- 
biased determination of total chlorophyll a (that is, total 
chlorophyll a equals divinyl chlorophyll a plus monovinyl 
chlorophyll a) for the purpose of ground-truthing SeaWiFS 
imagery. These co-eluting chlorophyll species can be re- 
solved spectrally following C18 HPLC chromatography 
(Wright et al. 1991) and quantified using dichromatic equa- 
tions at 436 and 450 nm (Goericke and Repeta 1993 and 
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Latasa et al. 1994). Alternatively, these two chlorophyll 
species can be separated chromatographically and individ- 
ually quantified using the C8 HPLC technique described 
by Goericke and Repeta (1993). (Cl 8 and C8 designate 
column packing materials used in HPLC.) 

These protocols, to be employed in the SeaWiFS valida- 
tion program for HPLC pigment analyses, are prescribed in 
the JGOFS Core Measurement Protocols (JGOFS 1991). 
These protocols include: 

a) use of Whatman GF/F glass fiber filters, ap- 
proximately 0.7 /an pore size, 

b) extraction in 90% acetone, and 

c) calibration with authenticated standards. 

It is recommended that seawater samples not be pre- 
filtered to remove large zooplankton and particles as this 
might result in the exclusion of pigment-containing colo- 
nial and chain- forming phytoplankton, e.g., diatoms and 
Tnchodesmium sp. Large zooplankton can be removed 
following filtration using forceps. 

When it is appropriate to measure phycoerythrin con- 
centration, it is important to discriminate between Tn- 
chodesmium and Synechococcus phycoerythrins, as they 
have different spectral absorption characteristics. The con- 
centration of Trichodesmium phycoerythrin may be deter- 
mined using the glycerol uncoupling method described by 
Wyman (1992). Separating Tnchodesmium from Syne- 
chococcus may require pre-filtering. 

In addition to HPLC analyses, it is recommended that 
the standard fluorometric methodology used for measur- 
ing chlorophylls and phaeopigments (Yentsch and Menzel 
1963, Holm-Hansen et al. 1965, and Strickland and Par- 
son 1972) also be applied to the same extracted pigment 
samples used for HPLC analysis (Section 6.6.1). This ad- 
ditional analysis will enable a direct link to the historical 
bio-optical algorithms and database development during 
the CZCS validation experiments. 

5*4.2 Spectrophotometric Absorption 

This section describes laboratory protocols for separat- 
ing the total spectral absorption coefficient, a(z, A), into 
its components due to absorption by pure water, in vivo 
phytoplankton pigments [a0(A)], dissolved material [a ff (A), 
GelbstofF] , and other particulate material [a t (A), tripton] 
by spectrophotometric measurements of filtered discrete 
water samples. Associated data analysis protocols are pre- 
sented in Section 6.8.2. 

All equipment utilized in water sample processing for 
GelbstofF samples should minimize contamination by or- 
ganic or colored material and should protect samples from 
photodegradation. Glass provides the best material for fil- 
tration units arid sample storage containers. Glass filtra- 
tion units should have stainless steel filter supports and not 
ground glass frits, since these tend to clog over time and 
change particle retention efficiencies of the units. Amber- 
colored borosilicate glass bottles are preferred for samples 


and filtrate collection because they protect the dissolved 
samples from ambient light. Caps lined with PTFE are 
preferred to minimize contamination. Prior to each exper- 
iment, all filtration and storage bottles should be thor- 
oughly cleaned with detergent and rinsed several times 
with distilled water (DIW). 

For measurement of absorption by dissolved organic 
matter (DOM), samples should be filtered through 0.2 pm 
membrane filters, e.g., Nuclepore™ polyester 0.2/rni fil- 
ters. Membrane filters minimize the contamination of the 
filtrate by filter fibers. It is recommended that separate 
water samples be filtered for DOM, particulate matter, and 
pigment analyses. Combining several analyses in one fil- 
tration can introduce errors, because each analysis has a 
specific filtering protocol. 

The procedure for preparing a sample for DOM mea- 
surements is as follows: 

1. Before filtering a sample for the measurement, 
filter at least 100 ml of DIW through the filter to 
remove wetting agents or other organic material. 

2. Swirl and rinse the filter trap. 

3. Filter 30-50 ml of the available volume, and cap- 
ture in a trap. 

4. Rinse the filter and flask a second time with a 
second aliquot of sample, or, if sample is scarce, 
refilter the first-pass filtrate. 

5. Swirl and rinse the trap after the second filtra- 
tion, and discard the filtrate. 

6. Filter the remaining volume, capture in the trap, 
and refilter the sample. 

7. Rinse a clean, amber glass, storage bottle three 
times with filtrate, then fill, cap and store up- 
right. 

If it is necessary to delay the measurement of the Gelb- 
stofF sample for more than a few hours, freeze each filtrate 
sample in an amber bottle, and store upright and in the 
dark. 

An alternate simple and rapid filtration method for 
measuring absorption spectra of dissolved organic matter 
in water samples has been proposed by Yentsch and Phin- 
ney (pers. comm.). A high-volume 0.22 pm Sterivex™ fil- 
ter cartridge! is used to perform a one-step filtration with 
the sample introduced directly into the cuvette from the 
filter cartridge. The sample is pushed through the car- 
tridge using a 50 cc sterile plastic syringe attached by a 
Luer lock connection. Sample handling is minimized and 
preparation time, including rinses, is only a few minutes 
per sample. The cartridges cost approximately $2.00 each, 
which makes them expensive to use for only one sample.’ 
Yentsch and Phinney found, however, that rinsing between 
samples with distilled water permits the use of one car- 
tridge per station with no measurable contamination ef- 
fects. They verified non-contamination by measuring the 

t A membrane filter made by Millipore Corporation. 
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absorption spectrum of a final distilled water rinse after 
all samples for a station were filtered, and then compar- 
ing the absorption spectrum of this sample to the baseline 
spectrum measured prior to the sample analysis. 

A potential disadvantage of the Yentsch and Phinney 
method results from the need to detach and reattach the 
cartridge to the syringe between sample and rinses. Back 
pressure applied while removing the syringe plunger with 
the cartridge attached can rupture the membrane. This 
problem can be avoided by using a peristaltic pump to in- 
troduce the sample and rinse water directly from separate 
bottles. 

The samples should be allowed to equilabrate to room 
temperature before spectrophotometer measurements are 
made, and they should be protected from ambient light. 
If the Yentsch and Phinney suggestion to inject the sam- 
ple directly into a cuvette is followed, this implies use of 
several cuvettes per station and some type of dark storage 
arrangement. Otherwise, the samples should be injected 
into amber bottles which have been prerinsed in the sam- 
ple filtrate. It may be possible to wait approximately 20 
minutes before filtering, to allow the samples to warm to 
room temperature, provided an investigator demonstrates 
that this delay does not allow modification of the sample 
filtrate. 

To summarize, the following method can be employed 
on samples collected and stored using any currently estab- 
lished technique: 

1. Prerinse a Sterivex cartridge using 100 ml of 
Nanopure or similarly polished DIW. 

2. For spectrophotometric baseline determination, 
rinse the reference and sample cuvettes 2 or 3 
times with filtered DIW and fill the cuvette, 
measure the absorption spectrum, and store it 
as the baseline. 

3. For each successive sample, evacuate the car- 
tridge, and rinse it with 50 ml DIW. Then, evac- 
uate the cartridge and rinse it with a 50 ml sam- 
ple. Refill the syringe with a sample, rinse the 
sample cuvette, and inject the sample from the 
syringe through the filter cartridge into the cu- 
vette. Measure the absorption spectrum in a 
spectrophotometer (see below). 

4. After all samples have been analyzed, rinse the 
cartridge with 50 ml DIW, refill the syringe with 
DIW, rinse the sample cuvette and fill it with 
filtered DIW, and measure the absorption spec- 
trum. Compare this absorption spectrum to the 
baseline. 

To measure the absorption of the dissolved fraction, 
minimize the distance between the detector and the cu- 
vettes in the spectrophotometer, and use 10 cm quartz cu- 
vettes to maximize the absorption signal. Warm the fil- 
trate samples to room temperature while protecting them 
from direct light. Filter reverse osmosis or Milli Q purified 


water (ROW) through a rinsed GF/F filter for reference 
as a blank, discarding the first 100 ml. Begin the analyses 
when the sample and reference material are both at room 
temperature. Perform a new baseline measurement range 
using a ROW blank each time the instrument has been 
powered up or the configuration changes. 

Using a dual beam spectrophotometer, scan the instru- 
ment’s baseline optical density spectrum OD fc (A) from ap- 
proximately 300-700 nm, beginning at the lowest ultravio- 
let (UV) wavelength the instrument will measure without 
cuvettes in the light path. Place the reference blank ROW, 
filtered as described above, in both sample and reference 
cuvettes and scan the reference optical density spectrum 
OD r (A). Rinse the sample cuvette with approximately 
5 ml of the sample three times. Rinse and fill the sample 
cuvette with the DOM sample filtrate and scan OD g (A). 
Replace reference water blanks every few samples to pre- 
vent bubble formation in the reference cuvette. Between 
samples, block the reference light and open the sample 
compartment door to minimize temperature changes in the 
ROW reference. 

The GF/F filter (which is binder-free and combustible, 
with a nominal pore size of 0.7 fim) is the best choice for 
particle absorption sampling. This filter has a large scat- 
tering coefficient at visible wavelengths, which increases 
the optical path length of the photons and improves the ac- 
curacy of particulate absorption measurements. This type 
of filter is also recommended by (JGOFS 1991) for carbon, 
hydrogen, and nitrogen (CHN); pigment; and primary pro- 
ductivity analyses. The optical transparency of the GF/F 
filter decreases significantly below 380 nm. 

Some authors have reported that particulate material 
less than 0.7/zm in size will not be retained by the GF/F 
filter, and that this fraction may contain up to 10-15% 
of the phytoplankton biomass as measured by chlorophyll 
concentration. Chavez et al. (1995), however, found no 
statistical difference between GF/F and 0.2 /rm filters for 
chlorophyll and productivity measurements. 

The goals for filtration of particulate samples are to 
minimize contamination and particle degradation, maxi- 
mize retention, and concentrate an adequate amount of 
particles on the filters to permit accurate spectrophoto 
metric measurements. Vacuum pressures below 5 inches of 
mercury, i.e., 117 mm Hg, are recommended to reduce the 
chances of particle breakage. It is currently recommended 
to set the filtration volume for each sample to approxi- 
mate an optical density value of 0.25 at 400 nm (Roesler 
and Perry 1992). Small area (13 mm 2 ) filtration units may 
be used to reduce water sample volumes. 

Optical density spectra of the filters should be mea- 
sured as soon as possible, because pigment decomposition 
may occur (Stramski 1990). If necessary to store filters, 
place the unfolded filters into polypropylene tissue capsule 
containers and store in liquid nitrogen. Freezing in conven- 
tional freezers causes ice crystals to form on the particles 
and leads to cell breakage and pigment degradation. 
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A new instrument baseline scan should be measured 
each time the spectrophotometer is powered up or any- 
time its configuration has been changed. To measure the 
particulate absorption spectrum, first prepare a blank fil- 
ter by soaking it in 50 ml DIW for 20 minutes (remember- 
ing to select it from the same box as the sample filters). 
An alternative approach is to use prefiltered seawater (pre- 
filtered using a 0.2 pm Millipore GS filter) to wet the GF/F 
filter blank (Cleveland pers. comm.), for the reason that 
the sample filters are also wetted with seawater. If using 
a single beam spectrophotometer, scan the instrumental 
baseline OD b (A) without the filter. Then insert the wet- 
ted filter blank to obtain the reference spectrum OD r (A). 
With a dual beam spectrophotometer, two wetted refer- 
ence filter blanks must be used to measure the reference 
spectrum, and one is left in the reference beam during sam- 
ple measurements. Some investigators recommend aligning 
the filter fibers in a fixed orientation each time a filter is 
inserted. 

If the samples were frozen, remove the sample filter 
from the liquid nitrogen. Place a few drops of DIW or 
filtered seawater on a glass (or Plexiglas) slide and then 
place the sample filter, bottom side down, on the water. 
This moistens the filter and hastens thawing. Place the 
filter in the sample compartment with the filtered mate- 
rial facing the source. Scan the OD p (A). Place the sample 
filter back in the filter cup. Gently pour 30-50 ml of hot 
spectrophotometric- or HPLC-grade methanol over the fil- 
ter, taking care not to lift any particulates off the filter 
(Kishino et al. 1985). Let the filter sit 10-15 minutes un- 
der the methanol, filter through, and repeat the process. 
Rinse with 10-20 ml of 0.2 ^m DIW. Seawater is not used 
here, because methanol may cause precipitates to form 
on the filter. Position the filter in the spectrophotome- 
ter and scan OD^(A). This spectrum is associated with 
raw, non-pigmented particulates, often referred to as de- 
tritus or tripton . The absorption spectrum should decrease 
exponentially with wavelength. If there is a residual chlo- 
rophyll q absorption peak at 676 nm, repeat the extraction 
process until the peak disappears. 

Recommended methods for analysis to convert the mea- 
sured optical density spectra to absorption spectra are dis- 
cussed in Section 6.8.2. 

5.4.3 Total Suspended Matter 

All suspended particulate material (SPM) dry weight 
(mgl -1 ) will be determined gravimetrically as outlined in 
Strickland and Parsons (1972) and as specified in JGOFS 
(1991). In general, samples are filtered through 0.4 ^m 
preweighed polycarbonate filters. The filters are washed 
with three 2.5-5.0ml aliquots of DIW and immediately 
dried, either in an oven at 75° C, or in a dessicator. The 
filters are then reweighed in a laboratory back on shore 
using an electrobalance with at least seven digits of preci- 
sion. 


5.5 ANCILLARY OBSERVATIONS 

Ancillary observations are often of key importance in 
flagging and interpreting apparently aberrant data. The 
minimal set of ancillary supporting observations must in- 
clude: 

1) date and time (both GMT and local); 

2) geographic location, using the Global Position- 
ing System (GPS) if possible, before and after 
each cast and at times of satellite and aircraft 
overpasses; 

3) solar azimuth and zenith angles, as calculated 
from position, date, and time; 

4) position of the optical cast in relation to the ship 
orientation and position of the ship relative to 
the sun (a sketch in the field notes is recom- 
mended); 

5) sea state (photographed if possible) with ap- 
proximate swell height, direction, and notes on 
presence and density of whitecaps; 

6) quantitative measurements of surface wire an- 
gles during deployments of the instrument pack- 
age; 

7) time of cast (begin and end), as well as time and 
depth of water samples collected; 

8) percent cloud cover and cloud type, and solar 
occlusion conditions; and 

9) wind direction and velocity. 

Desirable additional ancillary measurements include: 

a) an all-sky photograph plus a photographed time 
history of sea surface is advised for radiometric 
stations; and 

b) Secchi depth. 

5.6 PROTOTYPE OPTICAL BUOY 

A prototype optical measurement system designed for 
long-term buoy deployment with a satellite data telemetry 
capability is presently under development and is focused 
on satisfying the SeaWiFS optical data requirements. The 
concept is constrained by the requirement that the in- 
strument be capable of maintaining measurement integrity 
while being unattended for long periods of time. This 
constraint has lead to a design that minimizes the num- 
ber of moving parts to one, and has resulted in the spec- 
trographic application of concave holographic diffraction 
gratings. These holographic gratings provide an approxi- 
mate flat focal field to the degree that planar silicon pho- 
todiode arrays may be used as detectors. Inherent within 
this technology are the features of simplicity, compactness, 
durability, and stability. 

The optical system utilizes two spectrographs, in tan- 
dem with a dichroic mirror specifically designed to avoid 
measurements in the water absorption region, in order to 
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measure radiometric properties with high spectral resolu- 
tion and stray light rejection. The dichroic mirror is de- 
signed to transmit the red (630-900 nm) and reflect the 
blue portions (380-600 nm) of the spectrum, making the 

— IT fr ° m ri refleCtanCe t0 tran smittance between 600 
and 640 nm The potential for stray light is greatly re- 
duced by splitting the visible spectrum at the beginning 
of the water absorption region, since most of the short 
wavelength energy is diverted from the entrance slit of the 
long wavelength spectrograph. The splitting also allows 
the spectrographs re., gratings and sampling periods, to 
e optimized for the two distinctive spectral domains. A 
further reduction of stray light for the long wave spectra 
graph will be achieved by utilizing a minus blue filter 
The optical system will be deployed on a slack-line 
moored wave rider buoy that has a 10-20 m optical bench 
attached. Apparent optical properties will be measured 
by a series of remote collectors that are coupled to the 

storId men H f W u° PtiCS - ° ata wU1 be impressed, 

stored, and forwarded through a NOAA Geostationary Op- 
erational Environmental Satellite (GOES) and an ARGOS 
telemetry link. 

This type of optical mooring represents a new and chal- 
lenging technology. Detailed protocols for deploying and 
maintaining this type of mooring, and for evaluating its 
data quality, must be developed in light of the experience 
to be gamed over the next 2-3 years. 

The current practice is to apply marine antifouling com- 
pounds (such as OMP-8) to prevent the growth of marine 
organisms on the windows of moored radiometer systems 
This approach is less satisfactory for IOP instruments, be- 
cause for collimated light, transmission characteristics of 
the optical wmdows can be adversely affected by the layer 
of the anti-foulmg material. * 


5.7 DRIFTING OPTICAL BUOYS 

anH^h fting T iCal instruments are a recent development 

ratio h f e 18 t, m ° St n ° hlSt ° ry of their quantitative appli- 
cation to problems in ocean color algorithm development 

and remote sensing radiometric validation. It is probable 
that significant experience in the uses and limitations of 
such instruments will be gained in SeaWiFS related exper- 

Zmt1 Unng 1 Ke Prelaunch P eriod - The critical questions 
about these instruments, which should be answered during 
prelaunch work, include: g 

1. How accurately can L w { A) be estimated from 
L 'A z ' X ) at a single near surface depth, using 
on y an estimate of K L ( A) obtained from ocean 
color ratios? Is 5% uncertainty feasible? 

2. How accurately can a normalized water-leaving 
radiance ratio L WN {Ai)/ LwN ( \ 2 ) be estimated 
with these instruments using only a single chan- 
nel b d {\) measurement for normalization, which 
is contemplated for instruments currently being 
developed (Section 2.6)? Are the normalized 


water-leaving radiance ratios from clear sky and 
overcast conditions comparable enough that the 
drifter data can provide a basis for interpolating 
SeaW iFS data through cloudy periods? 

When answers to these fundamental questions are in hand 
it will be possible to draft and implement more detailed 
p otocols for the use of optical drifters in SeaWiFS radio- 
metric validation and algorithm development. 

Many potential applications of optical drifters in ocean- 
ographic research using SeaWiFS data are more obvious 
but from a radiometric standpoint are less stringently de- 
manding^ Protocols for those applications are, however 
beyond the scope of this report. 

5.8 SAMPLING AND VALIDATION 

The following discussion of bio-optical sampling pro- 
ocols is organized into three subtopics: sampling for the 
initial and ongoing validation of the SeaWiFS radiomet- 
ric system performance (Section 5.8.1), algorithm develop- 
ment and validation in Case-1 waters (Section 5.8.2), and 
algorithm development and validation in Case-2 waters 
.Section 5.8.3). The distinction between the first subtopic 
and the second two is clear-cut, but what precisely is meant 
by Case-1 and Case-2 water masses? 

In its literature and reports, the ocean color research 
community has formally adopted definitions originally due 
to Morel and Prieur (1977), who stated: Y 

Case- i is that of a concentration of phytoplankton 
(which is] high compared to that of other particles 
The pigments (chlorophyll, [and] carotenoids) play 
a major role in actual absorption. In contrast, the 
inorganic particles are dominant in Case-2, and pig- 
ment absorption is of comparatively minor impor- 
tance. In both cases, [the] dissolved yellow sub- 
stance is present in variable amounts and also con- 
tributes to total absorption. 51 

the pvrq ? iCe ’ bowever - only those water masses where 
the CZCS-type blue-green ratio algorithms for phytoplank- 
ton pigment concentration (chlorophyll a + phaeopigment 
a) work reasonably well have been treated as Case-1 All 
other water masses have often been loosely lumped into the 
Case-2 definition, albeit with considerable confusion over 
how to categorize coccolithophorid blooms, and similar 
phenomena normally classified as Case-1 waters, in which 
strong concentrations of Gelbstoff vary independently from 
chlorophyll a concentration. * 

In the present discussion of sampling protocols, Case-I 
will be considered to refer to what might be called ordinary 
open ocean Case-1 waters, wherein scattering and absorp- 
t on are dominated by phytoplankton, pigments, and Gelb- 
stoff concentrat'on 8 and where ghb&l b , ue _ green co]or r& _ 

a gorithms for chlorophyll a concentration and #(490) 

cTe W I' area f in the deep ° Cean belon S to this 

case. Water masses which do not satisfy this criteria will 
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be grouped under the heading Case-2. Within Case-2 by 
this definition, water masses with a wide diversity of 10 - 
optical characteristics will be found. Prominent subcate- 
gories include: 

1) coccolithophorid blooms, wherein the detached 
coccoliths dominate light scattering and remote 
sensing reflectance independently from pigment 
concentration; 

2) coastal areas, wherein DOM of terrestrial ori- 
gin contributes a strong absorption component 
which does not co-vary with pigment concentra- 
tion; 

3) phytoplankton blooms with unusual accessory 
pigment concentrations, e.g., red tides, which re- 
quire the use of special regional or local ocean 
color algorithms; and 

4) classical extreme Morel and Prieur (1977) Case- 
2 waters where optical properties are dominated 
by inorganic particles, with many possible varia- 
tions in chemical and geometric characteristics. 

It is important to recognize that some aspects of the 
water mass distinctions given above are dependent on the 
spectral regions in which measurements are to be made. 
Strong absorption at UV, red, and near-IR wavelengths 
requires the use of radiometric techniques similar to those 

required for Case-2 waters. * 

In addition to determining the bio-optical category a 
characteristics of a particular water mass, the validation 
sampling strategy must be concerned with spatial and tem- 
poral variability. Spatial and temporal variability in bio- 
optical properties will profoundly affect the validity of com- 
parisons between SeaWiFS and in-water optical measure- 
ments A single SeaWiFS instantaneous FOV measure- 
ment will integrate L W {\) over approximately a square 
kilometer, or a larger area at viewing angles away from 
nadir. Furthermore, the location uncertainty for a sing e 
pixel may be several kilometers, except in near-shore areas 
where image navigation can be improved by using land- 
navigated anchor points. 

Bio-optical profiles measured at a single station are rep- 
resentative of a spatial scale that is only a small fraction of 
a kilometer. Data from a grid of several station locations 
may be required to estimate the spatial averages of opti- 
cal properties represented by a SeaWiFS pixel, or a block 
of pixels. Because the ship measurements over the grid 
are not instantaneous, temporal variability in bio-optica 
properties can add additional uncertainty to the compar- 
isons. Aircraft radiometric observations can, conceptually, 
be used both to locate comparison sites away from areas of 
strong spatial variability and to document changes in t e 
pattern of spatial variability over the period required for a 
ship to occupy all stations in a comparison grid. 

Vertical stratification of water temperature, salmi y, 
and density often affect the vertical structure of variability 
in bio-optical properties. This variability, in turn, affects 


the remote sensing reflectance. Vertical stratification o 
the water column becomes especially important in many 
Case-2 waters, where the top attenuation depth may be as 
shallow as 1-2 m and the entire euphotic zone may be con- 
fined to less than 10 m depth. It is important therefore, 
to minimize ship-induced disruption of vertical st ratifica- 
tion in the water column. Whenever possible, the s p 
should be maneuvered as little as possible while on sta- 
tion with its propellers and bow thruster, and the prac- 
tice of backing down hard to stop quickly when on station 
should be strongly discouraged. If wind and sea conditions 
permit, the preferred method of approaching a station is 
to take enough speed off the ship to coast to a stop over 
approximately the last 0.5 km of approach to the station. 
The approach should be planned to allow the s ip ° e 
turned, preferably using only the rudder, to place the sun 
abaft the beam or off the stern, depending on where the 
radiometers will be deployed. It must be realized, how- 
ever that depending on wind and sea conditions, and a 
particular ship’s hull and superstructure configuration, l 
may not be possible to maintain an acceptable orientation, 
with respect to the sun, while the ship is adrift. In these 
situations, some use of the engines to maintain an accep - 
able ship’s heading may be unavoidable. 

The chief scientist should also consult with the ship s 
captain and chief engineer to avoid, or at least minimize, 
overboard discharges while the ship is on station. Material 
from a ship’s bilge or sewage treatment system can signi - 
cantly change near-surface chemical and optical properties 
if discharged near the immediate site of a bio-optical pro e 

or water sample. , . t , . 

In some coastal areas, where a relatively transparent 

water mass overlies a highly reflective bottom, LwW in- 
cludes light reflected from the sea floor. These cases require 
special treatment of bottom reflectance effects whether the 
local water mass regime is Case-1, Case-2, or a combina- 
tion of both. Methods of measurement, experiment design, 
and sampling strategies to study bottom reflectance effects 
are beyond the scope of this revision to the ocean optics 
protocols. There is a significant current research effort fo- 
cused in this area (Carder et al. 1992, Hamilton et al. 1992 
and Lee et al. 1994), and it is expected that experience will 
evolve the basis for new protocols in this topic area to be 
included in the next revision of this document. 

The bottom reflection of areas with a water depth ex- 
ceeding 30 m normally does not contribute to the water 
leaving radiance, L w { A)- Areas with a depth shallower 
than 30 m will be flagged in the SeaWiFS level -2 data 
product. Pixels covering very turbid waters may, however, 
even be usable even in shallower areas. As a genera ru e, 
the water depth should be deeper than 2.5 attenuation 
lengths, 1/AT(490), at all SeaWiFS algorithm development 
and validation stations. The prime exception to this rule 
is in developing local SeaWiFS algorithms where bottom 
reflectance contributions must be taken into account (Lee 
et al. 1994). 
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Variables to be measured at each validation station are 
summarized in Table 1, and methods for making each mea- 
surement are discussed in the previous subsections of Sec- 
tion 5. Methods of data analysis and reporting are dis- 
cussed in Section 6. 

5.8.1 Initialization and Validation 

Data intended for direct comparisons between observed 
L U (A) and SeaWiFS L w { A) estimates should usually be ac- 
quired in areas where bio-optical variability is known to be 
very small. This will ordinarily dictate that such data be 
acquired from optically clear and persistently oligotrophic 
Case-1 water masses. Potentially suitable sites include the 
northeastern Pacific central gyre off Baja, California (to 
the southwest), and the central Sargasso Sea. When plan- 
ning validation cruise locations and timing, seasonal and 
regional cloud cover statistics should also be considered in 
order to maximize the likelihood of simultaneous SeaWiFS 
and shipboard observations. A semi-oligotrophic site in 
the northeast Pacific, near Hawaii, is the prime candidate 
for placing a moored radiometer for continuous time-series 
radiometric comparisons with SeaWiFS Lyv(A) estimates. 

A series of radiometric comparison stations should be 
made over a wide range of latitude in both the Northern 
and Southern Hemispheres, to look for evidence of cyclic 
thermal sensitivity affecting SeaWiFS. The spacecraft and 
instrument will be heated by sunlight throughout the de- 
scending (daylight) data acquisition segment of each orbit 
and will be cooled by thermal radiation while in the Earth’s 
shadow throughout the remainder of the orbit. This cy- 
cling is likely to induce transient thermal gradients in the 
instrument, as well as a time varying cycle in the tempera- 
tures of its detectors and other components; these thermal 
variations could affect the spectral bandpass or responsiv- 
ity of one or more SeaWiFS channels. Unfortunately, a 
set of stations covering the full range of latitudes cannot 
all be sited in regions where mesoscale variability in ocean 
optical properties can be neglected. As when acquiring 
data for developing and validating Case-1 bio-optical algo- 
rithms (Section 5.8.2), a significant effort must be exerted 
to quantify spatial variability in normalized water-leaving 
radiance. When possible, airborne radiometer data, in 
combination with careful characterization of atmospheric 
aerosol and cloud conditions, should be employed to aug- 
ment shipboard radiometry at the stations selected for this 
aspect of the validation. If aircraft support is not available, 
semi-synoptic shipboard transects covering a 20x20 km 2 
grid should be used to characterize spatial bio-optical vari- 
ability near a sampling station. 

The minimum set of variables to be measured for radio- 
metric validation are those from the Product Verification 
column of Table 1. Measurements used to calculate nor- 
malized water-leaving radiance for direct comparison to 
SeaWiFS radiances must be made under cloud-free condi- 
tions and within five minutes of the satellite overpass. 


5.8.2 Case-1 Water Protocols 

In open-ocean oligotrophic water, it is usually practi- 
cal to assume that a station is in a Case-1 water mass, 
although some caution must be taken to detect coccol- 
ithophorid blooms and suspended coccoliths. In more tur- 
bid coastal transition regimes, however, the classification 
of the local water mass as Case-1 or Case-2 may be less ob- 
vious. In this environment, moreover, Case-1 and Case-2 
water masses may both be present in the domain sampled 
by a ship. One example of this situation would be Case- 
1 water within an eddy-like intrusion from offshore into 
coastal areas normally occupied by Case-2 water masses. 
Another would be Case-2 waters in a major river plume 
intruding into an ambient Case-1 water mass regime. 

In general, a water mass may be categorized as Case-1 
if: 

1) Gelbstoff absorption at 380 nm, a 9 (380), is less 
than 0.1m” 1 (Sections 5.2.4, 5.4.2, 6.8.1, and 
6 . 8 . 2 ); 

2) total SPM concentration is less than 0.5 mg l” 1 
(dry weight) (Section 5.4.4); 

3) measured Lw{ A) values, used in the SeaWiFS 
Case-1 algorithm, predict measured fluorometric 
chlorophyll a concentration within 35%; and 

4) measured Lw{ A), used in the SeaWiFS algo- 
rithm, predicts measured remote sensing K( 490) 
within 20%. 

The determination of criteria 1 and 2 above (Doerffer 
pers. comm.) will ordinarily require retrospective analy- 
sis. On the other hand, radiometric profiles and pigment 
samples can ordinarily be analyzed on board to allow de- 
termination of criteria 3 and 4 shortly after the samples 
are acquired. 

SeaWiFS Case-1 algorithm development and validation 
requires measurements from Case-1 water masses span- 
ning a wide range of optical properties and phytoplank- 
ton pigment concentrations. In optically transparent low- 
chlorophyll oligotrophic water masses, spatial variability is 
usually small and a station location and sampling strategy 
like that discussed in Section 5.8.1 is appropriate. 

In high-chlorophyll mesotrophic Case-1 water masses 
with increased turbidity, mesoscale and smaller scale vari- 
ability is often significant . In very productive Case-1 water 
masses, station placement and many other aspects of sam- 
pling schemes are similar to those appropriate for Case-2 
water masses (Section 5.8.3). At algorithm development 
stations, where measurements need neither be coincident 
with, nor matched to, SeaWiFS observations, it will be nec- 
essary to characterize spatial and temporal variability only 
over the relatively short scales distinguishing the separate 
in- water radiometric, optical, and pigment measurements. 
Airborne ocean color, or lidar characterizations, of spatial 
variability in the vicinity of these stations will not usually 
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be essential, although such additional information may be 
very helpful. 

At stations where data are acquired for algorithm val- 
idation, and where a match to concurrent SeaWiFS mea- 
surements is required, it will be necessary to determine 
the patterns of spatial variability over a domain extending 
approximately 20x20 km 2 centered at the station, and to 
place the ship in a 2x2 km 2 domain over which A(490) 
and chlorophyll concentrations vary less than 35% about 
the mean. Within a few hours before and after a Sea- 
WiFS overpass, in-water measurements should be made at 
several random locations to characterize variability within 
the 2x2 km 2 validation comparison site. In some cases, 
it may be possible to determine spatial variability ade- 
quately from ship station data and along-track measure- 
ments alone. In regions of strong mesoscale variability, 
however, concurrent aircraft ocean color or lidar measure- 
ments should be used both as a guide for selecting the 
ship’s location, and as a basis for spatially extrapolat- 
ing the in- water measurements to match the much coarser 
resolution of the SeaWiFS measurements (Section 5.1.7). 
These topics are developed in more detail, with respect to 
Case- 2 waters, in Section 5.8.3. 

5.8.3 Case-2 Water Protocols 

Although coastal and continental shelf areas comprise 
only 10% of the total ocean area, they provide roughly 
half of the oceanic new production and most of the se- 
questerable DOC (Walsh et al. 1981). These areas are 
typically higher in phytoplankton pigment concentration, 
and may include colored terrigenous constituents such as 
DOM and suspended sediments. In these Case-2 waters, 
the global color ratio algorithms break down because two 
or more substances with different optical properties are 
present which do not co-vary with chlorophyll a concen- 
tration (Section 5.8 above). These might be waters with 
exceptional plankton blooms (such as red tides), areas dis- 
colored by dust transported by the wind from deserts into 
the sea, or coastal areas influenced by river discharge of 
mineral and organic suspended materials, and DOM, i.e., 
Gelbstoff, such as humic acids. 

It is not always easy to decide to which case a wa- 
ter mass belongs. As a starting point, the water belongs 
to Case-2 if any of the four Case-1 criteria, set forth in 
Section 5.8.2, are not satisfied. For Case-2 waters defined 
by any one of these criteria, it remains a further problem 
to determine the specific bio-optical characteristics which 
distinguish it from Case-1. Case-2 sampling must usually 
include a more complete subset of the variables in Table 1, 
in addition to a ff (A) (Sections 5.2.4 and 5.4.2), SPM con- 
centration (Section 5.4.3), chlorophyll a concentration, and 
radiometry. For example, it may be necessary to determine 
complete pigment composition and other optically impor- 
tant characteristics of exceptional phytoplankton blooms 
for such planktonic groups as coccolithophorids, diatoms, 
cyanobacteria, or microflagellates. 


To achieve valid comparisons between the ship and sat- 
ellite data, sharp horizontal gradients and sub-pixel patchi- 
ness must be avoided, and image navigation must have land 
anchor points near the study site. Suitable landmarks are 
usually available in near-shore coastal waters. The other 
conditions are difficult to meet in Case-2 water masses, 
where mesoscale and sub-mesoscale variability is typically 
very strong. Sub-pixel variations of no more than ± 35% of 
the mean pixel chlorophyll will be tolerated, but variabil- 
ity must be measured and taken into account statistically 
in the analysis (see below). 

From the above generalities, it is clear that signifi- 
cant problems are encountered in near-shore coastal wa- 
ters characterized by small scale patches and dynamic vari- 
ability due to tidal currents. A particular problem occurs 
in the shallow areas which are influenced by strong tidal 
currents — areas that are normally well mixed during part 
of the tidal cycle. In the slack water tidal phase, how- 
ever, a vertical gradient of the suspended matter concen- 
tration may form, which may cause problems in relating 
water-leaving radiance to the concentration of suspended 
matter. During calm periods with strong insolation, even 
water that is normally well mixed can become stratified. 
In these cases, the formation of very dense phytoplankton 
blooms, such as red tides, can be observed. Such blooms 
will occur in coastal seas when nutrient concentrations are 
elevated by the influx of river water. In these circumstan- 
ces, it is especially critical to avoid disturbing the vertical 
stratification of the water column with the ship’s propellers 
(Section 5.8 above). 

5*8.4 Case-2 Sampling Strategy 

One approach to sampling in this environment has been 
suggested by R. Doerffer (pers. comm.). In order to get a 
good statistical base, water samples are first taken in a ran- 
dom order within the area under research. The concentra- 
tions derived from the SeaWiFS data are then compared 
with the ground truth data by statistical parameters, such 
as the mean, median, standard deviation, and the shapes 
of histograms (frequency distribution). For this type of 
statistical comparison, only sections of SeaWiFS images 
which match the area covered by the ship should be ana- 
lyzed. Water samples and satellite data should also be tem- 
porally concurrent within the same tidal phase in order to 
avoid biases due to temporal variability. In these regimes, 
analyses to validate algorithms cannot be based on Sea- 
WiFS data directly, but must instead be based on water- 
leaving radiance spectra measured in situ (Sections 5. LI- 
ST. 6) or from a ship (Section 5.1.8). This approach has 
the advantage that water samples and radiance spectra are 
taken nearly simultaneously. 

Using either flow-through pumping systems or systems 
towed outside the ship’s wake, fluorometry can be used 
to assess chlorophyll patchiness if frequent, i.e., every 10- 
15 minutes, chlorophyll fluorescence-yield calibration mea- 
surements are performed. Towed absorption, scattering, 
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reflectance, and beam transmission meters can also be used 
to characterize spatial variability. Within a few hours of 
the overpass, the ship should occupy several stations at 
random locations within a 2x2 km 2 area central to the 
area selected for comparison with SeaWiFS data. Sam- 
pling stations placed across a tidal front during a SeaWiFS 
overpass may help to identify two different water masses 
even when the front has moved. Comparisons between in 
situ and SeaWiFS data in patchy coastal areas may be 
enhanced by using horizontal radiance profiles measured 
from an aircraft flying at low altitude (Section 5.1.7). Sub- 
sets of such airborne profiles allow direct comparisons with 
shipboard data. A corresponding profile may then be ex- 
tracted from the SeaWiFS data for a direct comparison to 
the aircraft track line profiles. In Case-2 situations, such 
direct radiometric comparisons are valuable for validating 
and tuning local algorithms, but are not appropriate for 
SeaWiFS system validation per se. 

To validate the SeaWiFS atmospheric correction, water- 
leaving radiances measured in situ from the ship should 
be compared with those derived from the SeaWiFS data. 
Sample matching problems aside, Case-2 waters are of- 
ten characterized by strongly varying patchiness in optical 
properties, pigment concentrations, and remote sensing re- 
flectance at spatial scales smaller than a SeaWiFS pixel. 
Because of the nonlinear relationship between 65(A) /a( A) 
and remote sensing reflectance, the pigment concentration 
derived from spatially averaged SeaWiFS Lwn( A) values 
will systematically underestimate the true spatial average 
concentration by as much as a factor of 2 when sub-pixel 
variability is significant. It is, therefore, essential to de- 
scribe sub-pixel scale variability in Case-2 waters both 
statistically and in terms of organized structure. Such a 
description may be accomplished through rapid sampling 
at closely spaced ship stations in combination with air- 
borne ocean color measurements — for this purpose, track 
line data from low altitudes and high-resolution imagery 
from high altitudes are both acceptable. 

Absorption coefficients are large enough in all Case- 
2 waters to require instrument self-shadow corrections to 
L u (0",A), even though the correction model (Gordon and 
Ding 1992) has been experimentally verified only for the 
case where a(A)r is less than 0.1 (Sections 5.1.6 and 6.1.7). 
In extreme Case-2 waters, large values of spectral absorp- 
tion may confine the first optical attenuation depth to the 
top 1-2 m, where it is difficult to measure remote sens- 
ing reflectance in situ. Such short absorption scale lengths 


lead to instrument self-shading effects in L u ( A) which are 
correctable only for instruments with diameters no larger 
than approximately 1cm (Gordon and Ding 1992). Ra- 
diometers with such a small shadow cross section are con- 
ceptually feasible, and a few prototype instruments ex- 
ist which may be suitable, but they are not commercially 
available, and self-shadow sensitivities have not yet been 
experimentally verified for these extreme conditions. In 
these extreme cases, direct m situ measurements of a( A) 
and c(A) (Sections 5.2.4 and 6.8.1) and #,(#, A0, A, AA) 
(Section 5.2.5), together with above-water measurements 
of Lw(A) or may provide the only practical means 

of validating semi-analytic Case-2 algorithms. This topic 
is an important area for near-term research and develop- 
ment. 

5.9 VICARIOUS CALIBRATIONS 

An important obligation of any flight project is the 
production of a high quality, calibrated, Earth located 
(level-1) data set. Consequently, the production of a cali- 
brated set of SeaWiFS radiances that have been verified 
through direct, or vicarious, calibration techniques was 
recommended by the original workshop participants. 

One potentially useful technique follows the approach 
currently in use to verify the responsivity of the AVHRR 
Television and Infrared Observation Satellite (TIROS) sat- 
ellite instruments. Twice a year, an aircraft instrument, 
which has been recently calibrated directly to laboratory- 
based NIST traceable standards, should be used to obtain 
simultaneous views of a particular ocean scene. The air- 
craft scene must be obtained from a high altitude aircraft, 
such as the NASA Earth Resources-2 (ER-2), and flown at 
an altitude above most of the terrestrial atmosphere. The 
existing data sets of the AVHRR-NASA Aircraft/Satellite 
Instrument Calibration (NASIC) Project demonstrate a 
capability to limit the uncorrected trend in the AVHRR 
data sets to under 2% over two years. This concept may 
allow an independent verification of the atmospheric ra- 
diative transfer models used to compute ocean biological 
quantities, when the aircraft data are used in conjunc- 
tion with the surface truth campaign measurements of 
those same ocean biological quantities. The SeaWiFS ab- 
solute uncertainty requirements, however, are more strin- 
gent than those associated with the AVHRR, and a cor- 
respondingly more accurate airborne radiometer system 
(Sections 3.3 and 4.3) must be used (Section 2.3). 
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Chapter 6 


Analytical Methods 

Introduction 

This chapter describes methods for processing, analyzing, and reporting SeaWiFS validation data sets acquired 
following the protocols of Chapter 5. Where a consensus from the community exists for a preferred method, it 
has been adopted and protocols are given accordingly, e.g., for phytoplankton pigments (Section 6.5). In other 
instances, different methods are in current use, and as yet, there is insufficient evidence to give protocols specify- 
ing a preferred approach, e.g., the different methods of A'-analysis described in Section 6.1.4. Intercomparisons 
axe in progress to evaluate this particular topic. 


6.1 IN- WATER RADIOMETRY 

Spectral surface irradiance, E s ( A); upwelling irradiance, 
E u (z, A); downwelling irradiance, Ed(z, A); and upwelling 
radiance, L u (z, A), measurements should all be recorded 
and archived at four levels. 

a) level-0: raw instrument digital output; 

b) level-1: instrument output in volts, or frequency 
(if appropriate), and depth; 

c) level -2: calibrated irradiance and radiance, an- 
cillary measurements in appropriate geophysical 
or biological units, and depth, with ail variables 
corrected for dark or zero offsets; 

d) level-3: Smoothed profiles of K(z, A) and as- 
sociated Ed(z, A) or L u (z, A) with irradiance or 
radiance normalized by measured surface irradi- 
ance; and 

e) level -4: level -3 data normalized to clear-sky, a 
zenith sun at the mean Earth-sun distance, and 
spectrally adjusted to match the actual reference 
wavelengths and FWHM bandwidths. 

The formats of these data sets will vary somewhat be- 
tween individual instruments. It is anticipated that the 
SPO will promulgate suitable standard format specifica- 
tions, or guidelines, to facilitate database management and 
interchanges of level- 1, level -2, level -3, and level -4 data. 
These data files should each contain a header record iden- 
tifying as a minimum: 

1) date and time, i.e., GMT, of the station; 

2) geographic location (latitude and longitude in 
decimal degrees to the nearest 0.001); 

3) cloud cover and sky conditions; 

4) identification of each variable, including units 
and wavelengths, for radiometric channels; 

5) source of dark (zero-offset) data; 


6) calibration date and file identification; 

7) instrument identification; 

8) method for determining K (level -3); 

9) normalization algorithm (level -4); 

10) Secchi depth; 

11) depths of associated water samples, if any; and 

12) depth offsets (to nearest cm) between the pres- 
sure transducer and all sensor probes, including 
L u window, Ed and E u collectors, and all ancil- 
lary probes on a package. 

In addition to profile files, each data set should contain: 

i) calibration files used to compute level-2 data; 

ii) level -0 and level -1 dark files, and an average 
dark voltage file used for computing the corre- 
sponding level -2 files (in some cases a dark value 
may be extracted from the deep portion of a pro- 
file); 

iii) files with data from comparisons with a portable 
irradiance and radiance reference standard made 
in the field and used to track the instrument’s 
stability during a deployment; and 

iv) anecdotal and environmental information about 
each profile, either in the header, or in an ac- 
companying American Standard Code for Infor- 
mation Interchange (ASCII) text file. 

In general, the data should be retained at full reso- 
lution as recorded, but with noise-contaminated records 
removed through level -2. If the data are binned prior 
to A'-determination, as is sometimes done in the deriva- 
tive method (Section 6.1.4), the binned representations 
should be recorded as a level -2a file, in addition to the full 
resolution level -2 file. Depth offsets between the sensor 
port (where depth is determined) and the sensor apertures 
[Ed (A) and £ U (A) collectors and L u ( A) window] must be 
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recorded as header information in level -2 data files, and 
applied to the data prior to K-analysis (see Section 6.1.4 
below). 

6.1.1 Instrument Calibration Analysis 

Instrument data from pre- and post-deployment cali- 
brations should be compared with: 

1 ) each other; 

2 ) the long-term history of an instrument’s calibra- 
tions; and 

3) the record of comparisons with a portable field 
irradiance and radiance standard, to be made 
frequently during a cruise. 

Based on this analysis of the instrument’s history, a cali- 
bration file will be generated and applied to transform the 
data from level -1 to level- 2 . This analysis, and the ratio- 
nale for adopting a particular set of calibration coefficients, 
both for responsivity and wavelength, should be fully de- 
scribed in the documentation accompanying the data set, 
preferably in an ASCII file to be retained on line with each 
data set. 


6.1.3 Normalization by Surface Irradiance 

The dominant errors in measured K(z, A) profiles re- 
sult from changes in cloud cover. Cloud cover variabil- 
ity causes strong variations in incident surface irradiance, 
E 3 (X,t), measured at time t, during the time required to 
complete a radiometric cast. In present usage, E s (X,t) 
refers to incident spectral irradiance measured with a deck 
cell aboard a ship. Smith and Baker (1984 and 1986) dis- 
cuss a method for propagating E s ( A) through the sea sur- 
face to estimate £*(()-, A), and they also present a model 

for adjusting .Ed(0~,A) to compensate for solar zenith an- 
gle. 

An alternative, and conceptually better, scheme for es- 
timating E d ( 0 , A) is to measure E d (z r ,X) using a radiome- 
ter which has been floated away from the ship and held at 
a shallow depth, z r , during a cast (Waters et al. 1990). In 
either case, the record of E s (X,t) or E d {z r ,t) is recorded 
together with profiles of E d (z,X,t), E u (z, A), and L u (z, A). 
Assuming that transmission of E s ( A, t) through the surface 
does not vary with time, then a simple and effective nor- 
malization of the profiles is obtained as 


6.1.2 Raman Corrections 


E' (z X) = 

’ E 3 (X,t) 


(25) 


Marshall and Smith (1990), and the references cited 
therein, show transpectral Raman scattering contributes 
significantly to measured irradiance between 500-700 nm. 
At a particular wavelength, the Raman contribution is 
excited by ambient irradiance at a wavenumber shift of 
3,400 cm 1 . For example, Raman scattering at a wave- 
length of 500 nm (20,000 cm' !), is excited by light at wave- 
length 427 nm (23,400 cm' 1 ), and at 700 nm (14,286 cm' 1 ) 
by light at 565 nm (17,686 cm -1 ). Marshall and Smith 
(1990) give a transverse Raman scattering cross section fat 
90°) of 8.2 x 10 ~ 30 cm -2 molecule -1 sr' 1 , a value within 
the range of other published observations. By integration, 
they derive a total Raman scattering coefficient of: 

6 r (488) = 2.6 xlO^m' 1 . (23) 

The wavelength dependence of the Raman scattering cross 
section is theoretically about the same as that for Rayleigh 
scattering 

MA) = M488) AM , (24) 


where E s (X,t) is the deck cell irradiance measured at the 
time t when the radiometer was at depth z and E s { 0", A) is 
the measurement when the radiometer was at the surface. 

Some previous investigators have used E s {X,t) at a sin- 
gle reference wavelength, e.g., 550 nm, to normalize pro- 
files, and have thus ignored the usually small spectral vari- 
ations in incident irradiance. For SeaWiFS validation and 
algorithm development, however, the recommended proto- 
col is to use multispectral E s { A, t), or possibly near-surface 
Ed(z r ,X,t), to determine E d (z,X,t) at each wavelength. 

Because of spatial separation between the surface and 
underwater radiometers, cloud shadow variations are not 
measured, either identically or in phase, by the two instru- 
ments. The E s (X,t) or E d (z r ,X,t) profiles should, there- 
fore, be smoothed to remove high frequency fluctuations 
while retaining variations with periods of 15 seconds or 
greater. The smoothed E,(0~ , X)/E,{X,t) profiles should 
then be applied as a normalizing function to the irradiance 
and radiance profiles. 


although this has not yet been experimentally confirmed. 

A method for applying Raman corrections to measured 
profiles of irradiance and radiance is suggested and applied 
to homogeneous clear-water profiles by Marshall and Smith 
(1990). Additional work is needed to develop a robust 
Raman scattering correction model for general application 
in more turbid and vertically stratified water masses. The 
relative magnitude, and thus importance, of the Raman 
signal at each wavelength in the upper three attenuation 
lengths should also be investigated more thoroughly than 
has been done to date. 


6.1.4 Kr Analysis 

Normalized and Raman corrected profiles of E d {z y A), 
Eu{z, A), and L u (z , A) (with 2 corrected for pressure trans- 
ducer depth offset relative to each sensor) should be fit to 
the equations 

- f K d (z' ,\)dz' 

Ed(z, A) = E d (0~ , A) e 0 , (26) 

— J KxA{z\\)dz' 

E u {z,\) = £’ u (0“,A)e 0 } ( 27 ) 
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and 

- J K L (z\\)dz' 

L u (z, A) - L u ( 0-,A)e ° , (28) 

respectively. The vertical profiles of attenuation coeffi- 
cients K d {z, A), K u (z, A), and K L (z, A), together with the 
respective values of £d(0~,A), E u { 0 , A), and L u (0 } A) 
at the surface, provide the needed specifications for the 
smoothed irradiance and radiance profiles. 

If the natural logarithm of (26), (27), and (28) is taken, 
equations of the following form are obtained. 

z 

- J K{z)dz = ln(£(z)) - In (£(<T)), (29) 

0 


so that 


K(z) 


d ln(P(z)) 
dz 


z 


(30) 


The traditional method of A'-analysis, e.g., Smith and 
Baker (1984 and 1986), is to estimate K(z) as the local 
slope of measured ln[A(z)] in an interval of a few meters 
centered at depth z m , he., at depths near depth z m , 


at least 5-8 spectral bands. By using a multispectral ra- 
diometer (such as the MER class of instruments) during a 
profile, the data in all channels will be sampled contem- 
poraneously and recorded digitally 2-10 times per meter. 
These are level -1 data and are stored in files for subsequent 
processing and analysis. 

The level -2 through level -4 data give increasingly re- 
fined information in each successive processing level, which 
requires various amounts of intervention from the analyst. 
After appropriate editing to remove artifacts, such as the 
effects of ship shadow, vertical profiles of K are computed 
from the logarithmic decrement with depth of the radio- 
metric profiles. Direct derivative method calculations of 
K profiles using computer techniques (see above) require 
the use of a depth interval so large, frequently 20 m, that 
information about the slope, and hence, about K near the 
top and bottom of the profile, is lost. Averaging over such 
a large interval causes the slopes in sharply defined layers, 
e.g., regions of high gradients, to be poorly represented. 
Attempts to reduce these effects by using a significantly 
smaller depth interval result in unacceptably noisy K pro- 
files. 

An improved approach, suggested by Petzold (1988), is 
to fit a series of analytic functions to the radiometric data 


In [E{z)} “ In [E(z m )] “ ( z “ z m )K(z m )- (31) 

The unknowns ln[A(z m )] and K(z m ) are determined as 
the intercept and (negative) slope of a least-squares re- 
gression fit to measured ln[.E(z)] data within the depth 
interval z m - Az < z < + Az. The half-interval Az 

is somewhat arbitrary. Smith and Baker (1984 and i486) 
suggest a Az of approximately 4 m, but for noisy profiles, 
a Az as large as 10 m may be needed to smooth over in- 
cident irradiance fluctuations left as residuals by the deck 
cell normalization. 

When this method js used, the shallowest possible val- 
ues in the smoothed E(z) and K(z) profiles are at depth 
Az m, and the deepest values are Az m above the deepest 
measurements in the profile. If obvious ship shadow effects 
are present in the data, the shallowest valid smoothed data 
point will be at depth z a 4- Az, where z 5 is the depth to 
which the data are regarded as contaminated and are ex- 
cluded from the analysis. 

It is often convenient, although not necessary, to pre- 
average radiometric data into, e.g., 1 m, bins prior to per- 
forming the least-squares analysis. If this is done, the data 
should be pre-filtered to remove any noise spikes and then 
averaged before it is log-transformed. 

A corollary to having a large database is the need to 
facilitate its manipulation and analysis in order to make 
its application to various tasks feasible. In the proposed 
SeaWiFS effort, for example, radiometric measurements 
from many oceanographic stations will be examined. Each 
station will require one or more vertical profiles, from the 
surface to depths of up to 200 m, of downwelling irradiance, 
upwelling radiance, and upwelling irradiance, measured in 


using non-linear least-squares regression fitting techniques. 
The profiles are broken up into as many layers as required, 
and functions are fit to each layer, with the constraint that 
the functions and the derivatives of the functions be every- 
where continuous and finite. It is found that the logarithm 
of the radiometric data versus depth can be fit by a series 
of hyperbolic tangents superimposed on straight lines us- 
ing this technique. The data for a profile consisting of two 
layers can be matched by using the analytic expression 
together with the values of five parameters derived from 
the regression fitting procedure. In analyses of 2,100 pro- 
files, the most complicated profiles encountered required 
20 parameters — most required 5 or 10. 

With the analytic form of the curve fitting the data, 
it is a simple matter to differentiate the function to deter- 
mine the slope of the radiometric profile and obtain noise- 
free profiles of K. Using this technique, Petzold is able to 
store a very large database in a very compact form, stor- 
ing only the parameters and the program for reconstruct- 
ing the data. Additional analyses can be easily performed 
using the analytic representation of the data in lieu of the 
original large discrete data files. 

The basic functional form of the expression used to fit 
the data is 

y = p x + p 2 b + P 3 ^~ f • (32) 

A + A 

where B = P 4 - X, A = e B ! p \ Pi through P 5 are coeffi- 
cients to be determined in the analysis, X is the depth in 
meters, Y is the base 10 logarithm of the radiometric mea- 
surement, i.e., the downwelling irradiance (Ed), upwelling 
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irradiance (E u ) or upwelling radiance (L u ). The form of 
(32) is a hyperbolic tangent superimposed upon a straight 
line. It has a point of inflection at X = P±,Y — Pi and ap- 
proaches the asymptote Y — Pi + P^B ± P3 as X becomes 
larger or smaller than P4. 

The first derivative is 


dY 

dX 



(33) 


When X = P 4 ,B = 0, and Y = P u then 


dY 

dX 



(34) 


It is also possible to define K as the slope of the plot of 
the natural logarithm of the measured radiometric variable 
against depth, or 


K 


2.3026 


dY 

dX‘ 


(35) 


At depth K(z, A) is expressed as 

K(z,X) = tfo(Aboi(0 + d z K 0 (\hn(0 

+ tfi(Abo 2 (0 + a,^i(Abi2(0, 

where K o and K\ are values of K , and d z K o and d z K\ 
are its vertical derivatives at the two nodes of the depth 
element containing z. With this representation of K (2, A), 
it is possible to write (37) for each measured depth z m as 
the weighted sum 

N 

A) — ^ ^ -f* X n (39) 

n=0 

for the n = 0, 1, ... ,7V nodes dividing the water column 
into TV depth elements. The coefficients h x j are obtained as 
analytic integrals over the Hermitian polynomials 7 \j for 
the finite elements above and including depth 2; h t j = 0 
for elements below the one containing 2. Since such an 
equation may be written for every measured optical depth, 
the profile may be represented in matrix form as 


At the point of inflection ( X — P4 and Y = Pi), 


f = h k, 


(40) 


K 


2.3026 P 2 + 



(36) 


In the limit K — ► 2.3026 P2, and K will not exceed some 
finite value. 

The method described above was applied to radiomet- 
ric profile data from six cruises covering a wide variety of 
ocean regimes and latitudes from 24.0-77.4° N. Approxi- 
mately 2,100 profiles were fitted, and typically, the stan- 
dard deviation of the ratio between the function derived 
from the regression method and the original radiance and 
irradiance data was 6% or less. For a large fraction of the 
fitted profiles, the standard deviations were between 1-3%. 

An alternative method of determining A"-profiles was 
recently developed by Mueller (1991). Radiometric profiles 
are represented in terms of optical depth, r, which from 
(26), (27), and (28) is 


r(z,A) 


J K(z',X)dz' 


0 


In 


' £(0~, A) ~ 
. E(z, A) ^ 


(37) 


The Ff-profile is represented analytically by Hermitian cu- 
bic polynomials, 7fj(£), over finite depth elements. The 
argument £ is a local coordinate such that £ = 0 at the 
center of a finite depth element, £ = — 1 at the shallow end 
point (node) of the element and £ = 1 at the deep node. 
[Hermitian cubic polynomials are defined in any text on 
finite element modeling, e.g., Pinder and Gray (1977).] 


where f is the vector of measured optical depths, H is the 
matrix of coefficients h X j, and K is the vector of K n and 
d z K n at the TV nodes. The least-squares solution for the 

unknown vector K is obtained as 

= [H r E] _1 H T f, (41) 

which with (38) yields the complete profile K(z). 

The surface boundary condition assumed by Mueller 
(1991) is that K(z) is constant between the sea surface 
(node 0) and the first subsurface node (node 1). If obvious 
or suspected ship shadow effects are present in the upper 
profile, the depth of node 1 is set immediately below the 
affected area and the data in that top element are excluded 
from the fit; the solution to (41) at nodes 0 and 1 is, in this 
case, determined entirely by the data from depths below 
node 1. 

The solution at the deepest node is not constrained 
and depends only on the observations in the depth ele- 
ment immediately above it. The one-sided solution to (41) 
is often unstable at this node. If two nodes are placed close 
together at the bottom of the cast, then the unstable so- 
lution is confined to only the bottom node, which may be 
discarded after (41) is solved. 

In order to solve (41), the surface values of 1^(0", A), 
is u (0~,A), and L u (0 _ ,A) must be independently deter- 
mined or specified. At present, this is done iteratively 
by requiring the solution to closely approximate the mean 
value of measurements in the top 1-2 m of the profile. Data 
from this near-surface layer are usually not significantly af- 
fected by ship shadow, but they may be severely affected 
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by irradiance fluctuations associated with light focusing by 
surface waves. Additional research is needed to develop a 
more objective method of determining these surface values. 

In the present implementation (Mueller 1991), place- 
ment of nodes is largely subjective, even when guided by 
structure in accompanying c(660) and chlorophyll fluores- 
cence profiles. Qualitatively, the integral solutions mimic 
the structure in the c(660) and fluorescence profiles more 
faithfully than do the derivative solutions; they also do 
a better job of filtering irregularities that are apparently 
associated with large fluctuations in deck cell irradiance. 
Quantitative evaluation of sensitivity to exact node place- 
ment is in progress. Development of objective criteria for 
node placement will require further research. 

6.1.5 Finite Bandwidth Correction 

Siegel et al. (1986) and Marshall and Smith (1990) dis- 
cuss the effects of finite spectral FWHM bandwidth, and 
the normalized response function, on determination of the 
attenuation coefficient, K( A), for a vertically homogeneous 
water column. Given a channel's nominal wavelength, A', 
and normalized response function, h( A), the apparent at- 
tenuation coefficient measured in a homogeneous water col- 
umn is approximately 

CO 

f K{X)h{X)e~ K ^ z dX 

K s (z, A') = - 2 -35 • ( 42 ) 

/ h(X)e- K (» z dX 
0 

Marshall and Smith (1990) applied a correction for this 
effect to clear-water profiles of Ed(z,5 89). In general, cor- 
rection of K s (z, A') for finite bandwidth effects associated 
with K for pure water is straightforward. Additional re- 
search will be needed to model, from the spectral irradi- 
ance data itself, additional bandwidth effects associated 
with attenuation by phytoplankton and other particles, 
and to correct K s (z, A) accordingly. 

6.1.6 Extrapolation to the Sea Surface 

Because of surface waves, it is rarely possible to mea- 
sure Ed, E u , or L u at depths that closely approximate 
z = 0". The shallowest reliable readings typically oc- 
cur at depths ranging from 0.5-2 m. The data from this 
zone usually exhibit strong fluctuations associated with 
surface waves, and thus require some form of smoothing 
or averaging. It is almost always necessary to apply some 
means of extrapolating the data upward to the sea surface. 
Whatever method is used should reconcile extrapolated 
E d { 0~,A) with deck measurements of E s ( A). 

If K(z) profiles are determined using the derivative 
method, the shallowest smoothed estimates will occur at 
depth zo = Az, if there are no ship shadow effects. The 
usual procedure is to extrapolate values to z = 0 _ as 

E d ( 0-,A) = E d (z 0 ,X)e Kd{io ' X)zo , (43) 


Eu( 0-,A) 

= E u (z 0 , \)e K ^ Zo ' X)zo , 

(44) 

M o-,A) 

= L u {z 0 ,X)e I<L( - Zo ' X)zo - 

(45) 


If ship shadow is present, zq may be 20 m or more, and the 
extrapolation becomes somewhat tenuous. 

If K(z) profiles are determined by means of the inte- 
gral method, then £d(0~,A), E u (0~ ,X), and L u (0~, A) are 
automatically determined as part of the fitting procedure. 
The surface values thus obtained are not necessarily su- 
perior to those obtained by extrapolating the derivative 
method solutions, but they do have the advantage of rep- 
resenting an internally consistent fit to the entire profile 
beneath the surface boundary layer. 

By either method, extrapolation of measured Ed(z , A), 
E u (z, A), and L u (z, A) to z = 0“ becomes very difficult at 
A > 670 nm. At these wavelengths, the rapid decrease in 
daylight over an extremely shallow first attenuation length 
may compete with an increase in flux with depth due to 
chlorophyll fluorescence. Additional research is needed 
to address measurement and estimation of Ed(0~,X) and 
L u (0“,A) at these wavelengths, especially in chlorophyll- 
rich Case-2 waters. 

6.1.7 Instrument Self-Shading Corrections 

A provisional protocol is given here for radiometer self- 
shading corrections to L u (0~,X) and E u ( 0“,A) measure- 
ments (see also Section 5.1.6). The protocol is based on the 
model of Gordon and Ding (1992) and limited experimen- 
tal confirmation by Zibordi and Ferrari (1994). Although 
additional research is necessary to extend and verify these 
correction algorithms, the results published to date show 
clearly that even a provisional correction will significantly 
improve L u (0“, A) and F u (0“ , A) estimated from underwa- 
ter measurements. 

It is first necessary to estimate the spectral absorption 
coefficient a(A), preferably using measurements following 
the protocols of Sections 5.2.4, 5.4.2, 6.8.1, and 6.8.2. It is 
also possible to estimate a(A) using other approximations 
suggested by Gordon and Ding (1992), based either on 
measurements of phytoplankton pigment concentrations or 
of irradiance attenuation coefficients. 

It will also be necessary to measure, or estimate, the 
direct solar, E sun (X), and skylight, E s k y (A), components of 
incident spectral irradiance, E s ( A), where ^(A) = E sun (X) 
+£ s ky(A). Zibordi and Ferrari (1994) describe one method 
of estimating the ratio E s k y (X) / E sun (X) , and Gordon and 
Ding (1992) suggest other alternatives. 

Following Zibordi and Ferrari (1994), the coefficients, 
K f , given in Table 2 of Gordon and Ding (1992), are fit 
to linear regression models as functions of the solar zenith 
angle 9q in the range 30° < 9o < 70°. The results given 
for L u (0~,A), with sun only, for a point sensor may be 
computed as 

K ;un,o tan ^Ot« = 2.07 + 5.6 x 10 -3 #Oi (46) 
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and for a finite sensor occupying the full diameter of the 
instrument, 

Kiun^tanflou, = 159 + 6.3 x 10" 3 # 0 > (47) 

where 6o and 0ow are the solar zenith angles in air and 
water, respectively, measured in degrees. In practice, the 
diameter of the radiance sensor aperture is usually a small 
fraction of the instrument diameter. In the results reported 
by Zibordi and Ferrari (1994), the point sensor model al- 
ways overestimated e, and use of the finite sensor model 
(47) will always yield a lower estimate of e. Pending new 
insights from future theoretical and experimental work, it 
is suggested to estimate 


so that 

K sun = (1~/) K sun,0 + / K sun,J> ( 57 ) 

where / is the ratio of the diameter of the irradiance col- 
lector to that of the instrument. 

For the sky component, /t' ky is defined as 

*' ky = 2.70 - 0.48/. (58) 

Values of /c' un and < ky from (57) and (58) are then sub- 
stituted in equations (50) and (51) to obtain e S un(A) an ^ 
fsky(A), which are then used in (52) to solve for e(X). Fi- 
nally, corrected upwelled spectral irradiance E u (0“,A) is 
estimated as 


k' s un tan e 0 w = (l-/)«sun,O tan0 O + /<un,/ tan 6 l 0 , ( 48 ) 

where / is the ratio of sensor-to-instrument diameters. The 
coefficient, «' ky , for the self-shading effect on L u ( 0~, A) 
caused by incident diffuse skylight is similarly estimated 
as 

< ky = 4.61 - 0.87/, (49) 

where the coefficients are derived from values given in 
Table 3 of Gordon and Ding (1992). Self-shading errors 
e sun (A) and e sky (A) for F S un(A) and £ sky (A) components, 
respectively, are then computed as 

e su „(A) = 1 - (50) 


E u ( 0-,A) 


Ki o-,A) 

1 - *(*)’ 


(59) 


where E^(0~ , A) is determined from the upwelled spectral 
irradiance profile (Section 6.1.6). 

It is recommended that this correction algorithm be 
applied to all L u (0 ,A) and F u (0“, A) measurements used 
for SeaWiFS validation and algorithm development. Rec- 
ognizing the provisional nature of the correction, how- 
ever, the uncorrected measured values must also be re- 
ported. Moreover, the method and data used to estimate 
a(A),E S un(A), and E sky {X) must be documented and re- 
ported with all data sets corrected using this protocol. 


and 


Csky(A) 1 


e K' ky a(A)r^ 


(51) 


where r is the instrument radius in meters, and the ab- 
sorption coefficient a( A) is in units of m -1 . 

The self-shading error in L u (0~,A) is then calculated 


as 


e(A) 


^sun (A) e s ky(A)/i 

1 + h 


(52) 


where 

h ^skyW (53) 

F'sun(A) 

Finally, the corrected radiance L U (0~ ,X) is estimated 

< 51 > 

where L^(0~, A) is determined by analysis of the measured 
upwelled radiance profiles (Section 6.1.6). 

Similarly, for E u { CT,A), the values given in Tables 2 
and 3 of Gordon and Ding (1992) determine that for a 
point irradiance sensor, 


K 


sun,0 


3.41 - 1.55 x 10 -2 # 0 * 


(55) 


6.1.8 Spectral Adjustments 

New methods must be developed to reconcile in-water 
measurements L u ( 0“,Ao + AAi) integrated over a sensor 
response function hi(\) with SeaWiFS measurements of 
L t (Ao + AA 2 ) integrated over a wider sensor response func- 
tion h -2 ( A ) . The challenge is to account for differing ra- 
diometric sensitivities to fine-scale Fraunhofer structure 
in extraterrestrial solar spectral flux, Fo(A), as modified 
by atmospheric spectral transmittance, t( A), and oceanic 
spectral reflectance, R L { A). By assuming F 0 ( A) is exactly 
known, and that over the wavelength range defined by 
/ 12 (A) and / 11 (A), both t(X) and R L ( A) vary slowly with 
wavelength, it should be possible to adjust the Lw( A 0 + 
AAi) derived directly from the in- water instrument, to es- 
timate the water-leaving radiance L\y (A 0 + AA 2 ). This will 
be transmitted through the atmosphere and contribute to 
Lt { Ao + AA 2 ) measured by SeaWiFS. At the very least, 
this type of correction should be practical for a given at- 
mosphere. Prelaunch radiative transfer model sensitivity 
studies and experimental verifications should be done to 
determine the magnitudes and uncertainties of such cor- 
rections for the various SeaWiFS bands. 


For an irradiance collector with a diameter equal to that 
of the instrument, 

«sun,i = 2 76 - 1.21 X 10 _2 6>o, (56) 


6.1.9 Normalized Water- Leaving Radiance 

To standardize the in-water SeaWiFS algorithms, it is 
necessary to normalize measured Lw{ A) to those values 


53 


Ocean Optics Protocols for SeaWiFS Validation, Revision 1 


that would be measured were the sun at the zenith, at the 
mean Earth-sun distance, and with the effects of the at- 
mosphere removed. Following Gordon (1988), normalized 
water-leaving radiance as measured by the satellite sensor 
may be defined as 


Lwn(X) 


Lw / r_ \ 2 

£(A,0 o )(l - p(#o)) cos 0 O '-R' 


(60) 


where 9 0 is the solar zenith angle, p(9 0 ) is the air- water 
Fresnel reflectance for incident angle Q 0y t(\,0 o ) is the at- 
mospheric transmittance and R is the mean Earth-sun dis- 
tance. The Earth-sun distance on the day of the measure- 
ment, r, is given by 


1 + 0.0167 cos(^)’ 

where D is the sequential day of the year. 

The ( r/R ) 2 adjustment was not employed by Gordon 
and Clark (1981) or Gordon (1988) because it cancels in 
ratio algorithms, and the measurements they used were all 
taken within the span of a few months, so this source of 
variation w_as very small in their data. The range of varia- 
tion in (r/ R) 2 is approximately 6% over a full annual cycle. 
This adjustment should be made, nevertheless, for it be- 
comes important in algorithms predicting absolute values 
of Lw( A), as in the clear- water radiance model of Gordon 
and Clark (1981), and in algorithms for either estimating 
or detecting anomalously high water reflectances in, for 
example, a coccolithophore bloom. 

To obtain L\y{ A) it is necessary to propagate L u ( A) 
upward through the sea surface as 


L w (0\X) = L u (0-,X) — 2 P j*: g) , (62) 

n w (A J 

where p(X,9 and n u .(A) are the Fresnel reflectance and re- 
fractive index of sea water, respectively. Normalized water 
leaving radiance is then computed as 

L\vn{ A) = Lw{ A) (63) 

&S \ A) 

w r here Fq(X) denotes the mean extraterrestrial solar irra- 
diance (Neckel and Labs 1984). 


where p(A, 0) is Fresnel reflectance of the water surface. 
If L' ky (A,0,0 4- 7 r) is measured by viewing a horizontally 
oriented mirror of reflectance p rn (9 ), then (64) must be 
modified to 


L W (X) = L(X,e,4>) + AML L ' siiy (X,8,4> + 7r). (65) 

The remote sensing reflectance, Rl( 0 + , A), may be cal- 
culated from L\v(X) as 


Ri(o *’ A) = <66) 

where E s ( A) is incident spectral irradiance, preferably mea- 
sured concurrently with a calibrated irradiance meter, or 
estimated with greater uncertainty by viewing a horizon- 
tally oriented gray reflectance plaque (Section 5.1.8). Nor- 
malized water-leaving radiance is similarly calculated using 
(63). 

An alternative approach may be used to calculate re- 
flectance from an uncalibrated radiometer as 


R'l( 0 \X) 


S W (X) - p(X,9)S sky (X) 

7rS G (A)p 5 (A,0) 


(67) 


where 5 U ,(A), S sky (A), and Sg( A) are the radiometer sig- 
nals measured viewing the water, sky, and a reflectance 
plaque of approximately 10% reflectance (e.g., Spectralon), 
respectively; p(0) is the Fresnel reflectance of the surface; 
and p y (A) is the reflectance of the reference plaque. To 
remove residual surface reflectance due to wave facets, it 
is assumed that the residual signal at 750 nm is entirely 
due to surface reflection, so that the corrected reflectance 
is calculated as 


Rl( A) = R' l ( A) - i?' L (0+,750). (68) 

When using this technique, it is critical to measure all 
radiances well clear of the ship’s superstructure and to hold 
the reference plaque horizontally. 

It is more difficult to make verifiable corrections to 
L{ A, 9 , 4>) when measured from low flying aircraft (Sec- 
tion 5.1.7). In this case, atmospheric transmittance and 
path radiance are not entirely negligible, and water-leaving 
radiance must be calculated as 


6.2 ABOVE- WATER RADIANCE 

When upwelled radiance above the water surface is 
measured from a ship (Section 5.1.8), the measured ra- 
diance, L(A,0, 4>), must be corrected for the reflected sky 
radiance, L sky (A, 9, <f>). If T sky (A, #, </> 4- 7 r) is measured by 
directly viewing the sky, water-leaving radiance should be 
calculated as 

L\v{ A) = L(A, #,0) + p(A, #)L sky (A, 9, 4> 4- 7r), (64) 


L w {\,9,4>) = t- l [L{\,0,<t>)-L*{\6,<t>)} 

- p(A,0)L sky (A,0, (p 4- 7r), 

where t(X,9) is the atmospheric transmittance from the 
surface to the aircraft altitude, L*(A, 0,<f>) is the atmos- 
pheric path radiance at flight altitude, and L sky (A, 9, 0+tt) 
is the sky radiance incident at the surface. At altitudes 
up to 300 m, £(A,0) is close to 1, and may usually be es- 
timated with relatively low uncertainty. Unfortunately, 
there is no straightforward way to measure T aky (A,0, <fi + 
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7 r) or L*{\,0,<t>) from an aircraft. These variables must, 
therefore, normally be estimated using an atmospheric ra- 
diative transfer model. Verification of model estimates 
of L s ky {\,0,4> + 7 r) and L*{\,6,<p) is difficult under ideal 
cloud-free sky conditions, and is wholly impractical under 
variable cloud cover. 

Incident irradiance E d (z,\) measured at flight altitude 
must also be propagated to the surface to compute Rl(\) 
and LwnW fr° m airborne measurements, i.e., 



L w { A) 

RlW ~ E d (z,\)t d (z, A)’ 

(70) 

and 



, \ Fo(A) 

LwnW - LwnW Ed ( z ,\)t d (z,\y 

(71) 


where t d {z, A) is the downward spectral irradiance trans- 
mittance from flight altitude 2 to the surface. 

The authors defer any recommendation for specific pro- 
tocols for applying atmospheric radiative transfer models 
to estimate Lyv(A), Rl(A), and Lwtv(A) with (69)— (Tl) , re- 
spectively. It is incumbent on each investigator to report 
in detail the methods used to analyze (through measure- 
ments, models, or both) the necessary correction terms, 
and to include a detailed error analysis of the results which 
are obtained. Aircraft, radiometric data sets that do not 
include this information should not be used for SeaWiFS 
algorithm development or validation. 

6.3 MOORED RADIOMETRY 

Methods are not highly developed for analyzing data 
from moored radiometers to calculate LwnW- The prin- 
ciples of this analysis are well understood, but the commu- 
nity has had little experience with moored measurements 
of L u (z, A), determination of K L (z, A), and extrapolation 
to L u ( O', A). The moored optical system being developed 
by D. Clark of the NO A A National Environmental Satel- 
lite Data Information Service (NESDIS) for SeaWiFS and 
MODIS is the first system to be specifically engineered to 
address this problem. 

Smith et al. (1991) successfully acquired a nine-month 
time series of spectral Ed(z,\,t) and L u (z,A,t) at three 
depths (32, 52, and 72 m). They placed an additional 
above-water radiometer on a surface float, but this unit 
failed and provided no data. Smith et al. (1991) analyzed 
the Kd( 441) time series using a broad-band irradiance mea- 
surement to estimate E d ( 0~,441) for the 0-32 m depth in- 
terval, and using measured £d(z,441) over the 32-52 and 
52-72 m depth intervals. They did not, however, estimate 
K l ( A) for L u { A). They also developed and evaluated al- 
gorithms for estimating phytoplankton pigment concentra- 
tions from spectral reflectance and from naturally stim- 
ulated (by incident daylight) chlorophyll fluorescence at 
L u (683). Smith et al. (1991) demonstrated that continu- 
ous time series of K d ( A, t) and pigment concentration may 
be measured using this type of moored system. 


Smith et al. (1991) and Dickey et al. (1991) together 
illustrate methods that can be used to specify protocols 
for oceanographic analyses of bio-optical time series mea- 
sured using moored optical systems. Such protocols would 
be very valuable for planning and executing oceanographic 
studies using data from moored systems together with Sea- 
WiFS time series data; they are not, however, directly rel- 
evant to SeaWiFS validation. It is anticipated that optical 
protocols for US and International JGOFS will be pub- 
lished by working groups convened by these programs. 

6.4 AEROSOL OPTICAL DEPTH 

If multiple measurements of the solar beam are ob- 
tained during stable atmospheric conditions, then the Lan- 
gley method can be used to obtain the atmospheric trans- 
mittance. This method consists of plotting the natural 
logarithm of the voltage from the sun photometer versus 
the inverse of the cosine of the solar zenith angle. Ihe 
slope of this straight line is the total optical depth of one 
atmosphere. If only a single measurement is obtained, the 
instrument calibration is applied to determine radiance, 
which can be combined with the extraterrestrial solar ir- 
radiance to calculate the atmospheric optical depth. 

To obtain the aerosol optical depth, total optical depth 
must be used with computed optical depths due to molec- 
ular scattering (Rayleigh optical depth), and absorption 
by ozone and other important gases (NO 2 for some spec- 
tral bands). By subtracting the optical depths of these 
well-mixed gases from the total measurements, the aerosol 
optical depth can be determined. 

6.5 SKY RADIANCE 

Sky radiance distributions, L sky , measured with a cal- 
ibrated radiance distribution camera, perhaps augmented 
by sun photometry or narrow FOV L s k y discrete measure- 
ments in the zenith-sun plane, will be used to estimate 
the aerosol phase function (Voss and Zibordi 1989). De- 
velopment of detailed protocols and methods of analysis, 
including new inverse modeling techniques, for estimating 
aerosol optical depths and phase functions will require new 
research. The spectral mean cosine p^(0 + , A) for downwell- 
ing radiance at the sea surface will be calculated directly 
from radiance distribution camera data, when available. 
Under cloud-free conditions, ~p d (0 + ,\) can also be esti- 
mated by measuring £ s k y (A) 4* £ S un(A) with an irradiance 
deck cell. The algorithm for these computations is given 
by Gordon (1989b). 

When a spectral radiance distribution camera system 
is not available and skies are not cloud free, it may be 
possible to estimate Md(0 + ,A) f rorn some combination of 
deck cell unshaded £ s k y ( A) + E &nn (X) and shaded £ s ky(A) 
measurements, all-sky photographs, and measurements of 
L sky (A,0i,<M made at discrete angles with a hand-held 
radiometer. Additional research will be required to de- 
velop and test viable protocols for jxj(0 + ,A) estimation 
from these types of measurements. 
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6.6 PHYTOPLANKTON PIGMENTS 

This section specifies methods of analysis for determin- 
ing pigment concentrations by HPLC and fiuorometry (fol- 
lowing acetone extraction) from filtered water samples ob- 
tained at discrete depths, and for using in situ profiles of 
chlorophyll a fluorescence to interpolate pigment concen- 
tration over depth. 

6.6.1 HPLC Pigment Concentration 

The JGOFS protocols and standards for HPLC pig- 
ment concentration analysis (JGOFS 1991) will be the pri- 
mary method of determining pigment concentrations for all 
SeaWiFS algorithm development and validation activities. 

6.6.2 Fluorometric Determination 

Protocols for fluorometric determination of the concen- 
trations of chlorophyll and phaeopigments were developed 
initially by Yentsch and Menzel (1963) and Holm-Hansen 
et al. (1965), and are described in detail by Strickland and 
Parsons (1972). Although these measurements have been 
shown to contain errors as compared to HPLC determina- 
tions, e.g., Trees et ai. (1985), the CZCS phytoplankton 
pigment concentration algorithms were based on them en- 
tirely. The SeaWiFS protocols for this analysis will be 
those given in Strickland and Parsons (1972) as updated 
by Smith et al. (1981). 

6.6.3 In Situ Chlorophyll a Fluorescence 

In situ fluorometers produce nearly continuous profiles 
of artificially stimulated chlorophyll a fluorescence. Level -1 
fluorometer data (in volts) should be converted to level -2 
simply by subtracting an offset, determined by shading the 
instrument on deck. For qualitative guidance in K - profile 
analysis, level -2 (or even level- 1) fluorometer profiles are 
adequate. 

To produce vertical profiles of pigment concentration, 
HPLC-derived pigment concentrations from water samples 
taken at discrete depths should be interpolated, with the 
aid of in situ fluorescence profiles, for SeaWiFS bio-optical 
algorithm development. These fluorescence interpolated 
profiles should then be used with K d (z , A) profiles to com- 
pute optically weighted pigment concentration over the top 
attenuation length (Gordon and Clark 1980). 

6.7 BEAM ATTENUATION 

Raw beam transmissometer voltage profiles, V(z ), are 
first corrected for any range-dependent bias of the A/D 
data acquisition system (Section 5.2.1). The corrected 
voltages, V(z ), are then further adjusted for instrument 
drift (occurring subsequent to the factory calibration) with 
the equation 

V(z) = (v(z)-V datk )^, (72) 

' *air 


where Vd&rk is the instrument’s current dark response with 
the light path blocked, and V^ r and are, respectively, 
the current air calibration voltage (Section 5.2.1) and the 
air calibration voltage recorded when the instrument was 
calibrated at the factory. V(z) is then converted to trans- 
mittance, T(z, A), over the transmissometer’s path length, 
r, following the manufacturer’s instructions for the partic- 
ular instrument. 

The beam attenuation coefficient c(z, A) is then com- 
puted as 

c(z,A) = ~ 1 ln(T(z, A)), (73) 

which has units of m -1 . The apparent values of c(z,A) 
should be further corrected, again following the manufac- 
turer’s instructions, for the finite acceptance angle of the 
instrument’s receiver; this is usually a small, but signifi- 
cant, correction. Finally, the beam attenuation coefficient 
due to particles is computed as 

c P (z, A) = c(z, A) - c w (A), (74) 

where c w ( A) is the beam attenuation coefficient, i.e., the 
sum of absorption, a w ( A), and scattering, b w (A), for pure 
water. Smith and Baker (1981) tabulate a w (\) and b w ( A) 
over the spectral range of interest here. 

6.8 SPECTRAL ABSORPTION 

Protocols are given above in Section 5.2.4 for measuring 
absorption profiles in situ using reflective tube meters, and 
in Section 5.4.2 for laboratory spectrophotometric analyses 
of water samples captured at discrete depths. 

Data from a reflective tube absorption and beam atten- 
uation meter may be analyzed to obtain vertical profiles of 
a (z, A), a g (z, A), and c(z,A), and by difference 6(z, A) = 
c(z, A) — a(z, A) and a p (z, A) = a(z, A) — a^(A, t), using the 
protocols of Section 6.8.1. Optical density spectra for fil- 
trate and filtered water samples (Section 5.4.2) may be an- 
alyzed to obtain independent measures of a g (z , A), a p (z , A), 
and a t (\, t), and by difference a/ (z, A) = a p (z, A)-a*(z, A), 
using the protocols of Section 6.8.2. Methods for merging 
and comparing the two independent types of absorption 
measurements, and for interpreting the results in terms of 
remote sensing reflectance, are the subject of currently ac- 
tive research by several investigators. The next revision 
to this document may be expected to contain extensive 
modifications and extensions of these protocols. 

6.8.1 Reflective Tube Measurements 

At the time of this writing, there is only one commer- 
cially available reflective tube absorption and beam atten- 
uation meter, the AC-9, manufactured by Western Envi- 
ronmental Technology, Inc., of Philomath, Oregon. The 
following data analysis protocol is, therefore, based on the 
characteristics of that instrument (Zaneveld and Moore 
pers. comm.). 
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The AC-9 output in each a( A) or c( A) channel consists 
of the digital responses of a reference detector and a wa- 
ter transmission detector, C r ( A) and C t ( A). During the 
instrument’s factory calibration using optically pure water 
(Sections 4.4 and 5.2.4), the responses C"( A) and C'(A) 
are recorded and entered on the instrument’s calibration 
sheet, along with the instrument temperature V at which 
the calibration was performed. Calibration factors for each 
channel are then determined as 

A( x a ) = ( 75 ) 

for absorption at wavelength A a , and 

a(\ c ) = ^T) e ~ c "' (Xr)x ' (76) 


for beam attenuation at wavelength A c , where x is the 
optical pathlength in meters, and au;(A a ) and c w (A c ) are 
the absorption and beam attenuation coefficients of opti- 
cally pure reference water, respectively (Smith and Baker 
1981). The instrument’s factory calibration also includes 
determination of a temperature coefficient £*(A n ) for each 
channel. 

In the subsequent discussion in this section, the sub- 
script i designates initial values of the absorption beam at- 
tenuation coefficients derived directly from measurements. 
The subscript e denotes the corrected values of the same 
coefficients due to (extra) substances other than water. 

When data are recorded in the field, apparent absorp- 
tion coefficients a z (A a ) due to substances other than pure 
water are calculated as 


a t {X a ,T) = - In 


C t (\ a ,T) 


TA( A a ) 


C r (A ai T) 

+ (T'-T)B t (X a ) 

+ [ a air(A a ) — G air(^o)] 


(77) 


where V and T denote the instrument temperature, and 
a^ r (A 0 ) and a air (Ao) are air calibration values at the times 
of the instrument factory calibration and field data record- 
ing, respectively (Section 5.2.4). A similar equation is 
applied to obtain apparent beam attenuation coefficients 
c 2 (A c ) for substances other than water at each wavelength. 

Apparent absorption coefficients a*(A,T) must be cor- 
rected for the fraction of scattering that is undetected by 
the reflective tube and detector assembly (Section 5.2.4). 
An initial estimate of particle scattering at each wave- 
length bi(X) is obtained as 


bi{ A) = Ci( A) - aj(A,T), (78) 


recalling that attenuation by pure water already has been 
removed from c z (A) and a, (A ,T) using the instrument’s 
calibration coefficients. It is assumed that: 


1) the uncertainty in c t (A) due to forward scattered 
light, which is viewed by the detector (Section 6.7), 
is a wavelength-independent fraction k c of total par- 
ticle scattering b p ( A); 

2) the fraction h of 6 p (A) included in a t (A) is also in- 
dependent of wavelength; and 

3) at a near-IR wavelength A r , absorption is entirely 
due to pure water, and the apparent absorption 
a l (\ r ,T) is completely due to undetected scatter- 
ing kb p (X r ) and temperature dependent variations 
in water absorption (Pegau and Zaneveld 1993) rel- 
ative to absorption at T f (the calibration tempera- 
ture). 

From assumptions 1) and 2), it follows that 

bp(X i ) _ bi(X\) (79) 

M a 2 ) “ b t (x 2 y 

From assumption 3), it follows that 

, _ a t (A r ,T) 4- S(X r )(T - T') 

b t ( X r ) ’ V ' 

where 5( A r ) represents the coefficient of water tempera- 
ture variation in a w (X,T). At A r = 712 nm, for example, 
5(712 = 0.0035 m -1 °C _1 , and 5(A) is negligible at A < 
700 nm (Pegau and Zaneveld 1993). It is then useful to 
combine (79) and (80) to determine the absorption coeffi- 
cient a e (A) due to substances other than water (Zaneveld 
pers. comm.) as 


a e (A) = a z (A,T) 

[^(Ar.n-s^xr-roj&^A) (si) 

MA r ) 

When equations (75) through (81) are applied to unfiltered 
total absorption measurements, it follows that 

a e (A) = n p (A) -1- a g ( A), (82) 

where a p (A) and a g (X) are respectively the particulate and 
Gelbstoff (dissolved substances) components of absorption. 
If a 0.2 pm intake filter canister is used to exclude the par- 
ticulate fraction from the water, then a e (A) is due to a g ( A) 
alone. If profiles are made with two AC-9 meters (one fil- 
tered and one unfiltered), the total absorption coefficient 
may, therefore, be partitioned as 

a(X) — d w (X) + &p(A) + (83) 

To further partition a p (A) into a^(A) ( phytoplankton pig- 
ments) and a t ( A) (tripton) absorption coefficient compo- 
nents, it is necessary to obtain discrete water samples at 
several depths distributed over the profiles, and to apply 
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the filtration and spectrophotometry laboratory analyses 
described in Sections 5.4.2 and 6.8.2. 

The fraction of total particle scattering, k c b p ( A), in- 
cluded in the beam transmission measurement may be esti- 
mated from the detector acceptance angle and path length 
geometry of a particular instrument. Using this coefficient, 
the corrected non- water beam attenuation coefficient c e (A) 
may be calculated as 


some investigators have used apparent optical densities at 
750 nm for this correction, but recent work has shown this 
wavelength to be inappropriate in some oceanic waters. 

The prevalent current practice is to estimate (3 empiri- 
cally through a function that may be expressed in the form 

/T 1 = Ci + C 2 [OD p (A)-OD r (A)], (90) 


Ce(A) = Cj(A) + fc c [c e (A) -a e (A)], 

(84) 

or 


„ t w g(A) fc c a e (A) 

Q “ 1 - k c 

(85) 

Therefore, the corrected particle scattering 
may be estimated as 

coefficient 

MA) ^e(A) n e (A). 

(86) 


6.8.2 Analysis of Absorption Spectra 

Spectral absorption coefficients a g ( A) (Gelbstoff), a p 
(total particulates), and a t (tripton) in units of m -1 are 
calculated from optical density spectra measured using the 
protocols of Section 5.4.2. For Gelbstoff, 

2 303 

M A ) = ~ — [OD s (A) - OD r (A)], (87) 

where l is the cuvette pathiength (0.1 rn), and the other 
terms are introduced in Section 5.4.2. 

An equation similar to (87) is used to calculate ab- 
sorption coefficients a p (A) and a t { A), but in these cases, 
a correction must be applied to account for the effective 
increase in pathiength due to scattering in the GF/F fil- 
ters (Mitchell and Kiefer 1984). The nominal absorption 
pathiength l s of the filtered material in suspension is given 

by 

l ° = A' (88) 

where V is the volume of water filtered and A is the clear- 
ance area of the filter. Scattering of light within the GF/F 
filter lengthens the absorption pathiength. This factor is 
termed the (3 factor by Mitchell and Kiefer (1988) after 
Butler (1962), which is not to be confused with the volume 
scattering coefficient (3{6) used elsewhere throughout this 
document. The absorption coefficient of filtered particles 
must, therefore, be corrected for pathiength amplification 
and the equivalent absorption in suspension is computed 
as 

2 3034 

M A ) = —pv [ODp(A) — OD r (A)] . (89) 


where C\ and C 2 are coefficients of a polynomial regression 
fit. For example, Mitchell (1990) gives C x = 0.392 and 
C ‘2 = 0.665, while Cleveland and Weidemann (1993) give 
Ci = 0.378 and C 2 = 0.523. 

Spectrophotometers that disperse white light into its 
spectrum using a monochromator prior to illuminating the 
sample are termed pre-sample monochromators, while in- 
struments that pass the light through a monochromator 
after it has passed through the sample are termed post- 
sample monochromators. The value of the pathiength 
amplification factor has been found to be independent of 
post-sample optical geometry for a variety of standard 
laboratory instruments with pre-sample monochromators 
(Mitchell 1990, Stramski pers. comm., and Cleveland pers. 
comm.). Clark (pers. comm.) and Cleveland (pers. comm.) 
have found, however, that diode array spectrophotometers 
with a post-sample grating require the use of an integrating 
sphere attachment which may affect the /? factor. Lee et al. 
(1994) use a diode array system, which collects and focuses 
the post-sample light and gives raw OD p values similar 
to standard instruments with pre-sample monochromators 
(Cleveland pers. comm, and Carder pers. comm.). 

Substituting (90) in (89) directly yields the spectral 
absorption coefficient for particles in suspension as 


a p (A) — 


2.303,4 


[ODp(A) - OD r (A)] 

+ C 2 [OD p (\) - OD r (\)Y 


(91) 


and for tripton 


a p( A) — 


2.303,4 


C\ [OD t (A) — OD r (A)] 

+ C 2 [OD t (A)-OD r (A)] : 


(92) 


Finally, the spectral absorption coefficient for in vivo phy- 
toplankton pigments is computed as 

n^(A) = a p (A) ut(A). (93) 


To correct for residual offsets in the filter spectrophoto- The best overall performance of this algorithm will be 
metric measurement, it is assumed that at 850 nm the true achieved when filtered sample density yields OD (675) be- 
optical densities are zero. The apparent measured values tween 0.05 and 0.1. The volume of water filtered^or parti- 
of OD p (850) and OD r (850) are, therefore, subtracted from cle absorption measurements should therefore be adjusted 
ODp(A) and OD r (A) prior to applying (89). Previously, accordingly (Section 5.4.2). 
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For purposes of data reporting and archiving, the ab- 
sorption coefficients a 9 (A),a p (A), and a*(A) will be reported 
in m _1 as computed using equations (87), (91), and (92), 
respectively. Uncorrected optical density spectra for the 
blank filters and distilled water blank must be recorded 
and provided with absorption coefficients to allow reversal 
of the correction algorithms. The values of OD p (850) and 
OD r (850), which have been subtracted to remove resid- 
ual instrument offset, must also be reported with each ab- 
sorption spectrum. The pathlength amplification factor, 
together with a description of the method and date of its 
determination, must also be reported with all particle ab- 
sorption coefficients based on filter spectrophotometry. 

Roesler and Perry (pers. comm.) have proposed a theo- 
retically derived expression for /3, with a formulation differ- 
ent from (90). At this writing, however, this approach has 
not gained widespread acceptance in the community. Fur- 
ther review of this protocol is warranted for possible future 
revision. Sosik and Mitchell (1994) have proposed a spec- 
trofluorescence excitation method to make direct estimates 
of absorption by the photosynthetically active pigments 
[°p«(A)] ■ They found a ps (A) < a*(A), where the difference 
is attributable to photo- protective pigments in the phyto- 
plankton. Carder and some of his associates (pers. comm.) 
recommend the determination of ajf (666) in the methanol 
extract using a spectrophotometric method at 666 nm with 
850 nm zero. This method sets an upper limit to a^>(675) 
for the unpackaged state, and allows an alternate estimate 
of chlorophyll a. 

6,9 CTD PROFILE ANALYSES 

Each CTD profile should be prefiltered to remove any 
depth reversal segments resulting from violent ship or hy- 
drowire motions. This will remove many instances of salin- 


ity spiking , an artifact which occurs when water tempera- 
ture changes at a rate faster than the conductivity probe 
can follow. The CTD data should then be processed to 
profiles of potential temperature (°C), salinity (Practical 
Salinity Units [PSU] based on the Practical Salinity Scale 
of 1978, PSS78), and density (kgm -3 ) using the algo- 
rithms which have been endorsed by the United Nations 
Educational, Scientific, and Cultural Organization (UN- 
ESCO)/SCOR/International Council of Exploration of the 
Seas (ICES)/IAPSO Joint Panel on Oceanographic Tables 
and Standards, and also by SCOR Working Group 51 (Fo- 
fonoff and Millard 1983). 

At this stage, each set of CTD profiles should be care- 
fully examined to detect any significant static instability 
artifacts resulting from salinity spiking. After any such 
major artifacts are removed by editing, the data should be 
further smoothed by averaging temperature and conduc- 
tivity data into 2 m depth bins, and the final profiles of 
salinity, density, and other derived parameters should be 
recomputed using the smoothed CTD profile. 

For any hydrographic station, descriptive hydrographic 
analyses should include T-S profile characterizations of wa- 
ter masses. Features in the density profile which appear 
to be related to physical mixing and stability should be 
compared with features in the corresponding bio-optical 
profiles. CTD profiles from horizontal transects (i.e., two- 
dimensional grids) should be used in the computation of 
tw r o-dimensional sections, or three-dimensional gridded ar- 
rays, for such variables as geostrophic currents, tempera- 
ture, salinity, and the density anomaly a t . These analy- 
sis products, together with corresponding two- or three- 
dimensional representations of bio-optical variability, can 
be used to estimate the relative importance of advection 
and isopycnal mixing in redistributing or modifying upper 
ocean optical properties during a cruise. 
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Glossary 

A/D Analog- to- Digital 

ALSCAT ALPHA and Scattering Meter (Note: the symbol 
a corresponds to c(A), the beam attenuation coeffi- 
cient, in present usage.) 

AOL Airborne Oceanographic Lidar 
ARGOS Not an acronym: the name given to the data col- 
lection and location system on NOAA Operational 
Satellites 

ASCII American Standard Code for Information Inter- 
change 

AVHRR Advanced Very High Resolution Radiometer 
AVIRIS Advanced Visible and Infrared Imaging Spectrome- 
ter 

BSI Biospherical Instruments, Inc. 

CDOM Colored Dissolved Organic Material 
CHN Carbon, Hydrogen, and Nitrogen 
CTD Conductivity, Temperature, and Depth 
CW Continuous Wave 
CZCS Coastal Zone Color Scanner 

DIW Distilled Water 

DOC Dissolved Organic Carbon 

DOM Dissolved Organic Matter 

ER-2 Earth Resources-2, a research aircraft 

FEL Not an acronym; a type of standard lamp for irra- 
diance and radiance calibration 
FOV Field-of-View 
FWHM Full-Width at Half-Maximum 

GAC Global Area Coverage 
GASM General Angle Scattering Meter 
GF/F Not an acronym; a specific type of glass fiber filter 
manufactured by Whatman 
GMT Greenwich Mean Time 

GOES Geostationary Operational Environmental Satellite 
GPS Global Positioning System 
GSFC Goddard Space Flight Center 

HPLC High Performance Liquid Chromatography 

IAPSO International Association for the Physical Sciences 
of the Ocean 

ICES International Council on Exploration of the Seas 
IFOV Instantaneous field-of-view 
IOP Inherent Optical Properties 
IR Infrared 

JGOFS Joint Global Ocean Flux Study 

MARS Multispectral Airborne Radiometer System 
MER Marine Environmental Radiometer 
MERIS Marine Environment Research Imaging Spectrora- 
diometer (French) 

MODIS Moderate Resolution Imaging Spectrometer 

NAS National Academy of Science 
NASA National Aeronautics and Space Administration 
NASIC NASA Aircraft/Satellite Instrument Calibration 
NESDIS National Environmental Satellite Data Information 
Service 

NIST National Institute of Standards and Technology 
NOAA National Oceanic and Atmospheric Administration 
NOARL Naval Oceanographic and Atmospheric Research 
Laboratory 
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OCTS Ocean Color and Temperature Sensor (Japanese) 
OFFI Optical Free-Fall Instrument 
OMP-8 Not an acronym; a type of marine anti-biofouling 
compound 

OSFI Optical Surface Floating Instrument 

PAR Photosynthetically Available Radiation 
POC Particulate Organic Carbon 
POLDER Polarization and Directionality of the Earth Re- 
flectance (French) 

PON Particulate Organic Nitrogen 
PSU Practical Salinity Units 
PTFE Polytetrafluoroethylene, commonly known by the 
trade name Teflon® 

QED Quantum Efficient Device 

ROSIS Remote Ocean Sensing Imaging Spectrometer, also 
known as the Reflecting Optics System Imaging 
Spectrometer (German) 

ROV Remotely Operated Vehicle 
ROW Reverse Osmosis Water 

SCOR Scientific Committee on Oceanographic Research 
SeaWiFS Sea-viewing Wide Field-of-view Sensor 
SIRREX SeaW it S Intercalibration Round-Robin Experiment. 
SNR Signal-to-Noise Ratio 
SPM Suspended Particulate Material 
SPO SeaWiFS Project Office 
SPSWG SeaWiFS Prelaunch Science Working Group 
SST Sea Surface T ’em per at ure 

' F- S Tern p era t u re- S al i n i 1 y 
TIROS Television Infrared Observation Satellite 
TSM Total Suspended Material 

UNESCO United Nations Educational, Scientific, and Cul- 
tural Organizations 
UV Ultraviolet 
UVB Ultraviolet- B 

WMO World Meteorological Organization 

Symbols 

A Fitting coefficient for P 4 - AC or clearance area of 
a filter, depending on usage. 

A(X a ) AC-9 instrument calibration factor for absorption. 
A(A C ) AC-9 instrument calibration factor for beam atten- 
uation. 

a (A) Total absorption coefficient. 
a(z, A) Spectral absorption coefficient. 

(A) Absorption coefficient due to substances other than 
water. 

a/(z,X) a p (A) - a f (z, A). 
a gW Gelbstoff spectral absorption coefficient. 

Ot(A a , 7 ) Initial estimate of the apparent absorption coeffi- 
cient; used for determining the apparent absorption 
coefficient for substances other than water. 
a 4 > (^) Phytoplankton pigment spectral absorption coeffi- 
cient determined in methanol extract. 
a p(A) Particulate spectral absorption coefficient. 
a pa (A) Photosynthetically active pigment spectral absorp- 
tion coefficient. 

a t(A) Tripton spectral absorption coefficient. 
a w(X) Absorption coefficient for pure water. 
a <t>W Phytoplankton pigment spectral absorption coeffi- 
cient. 


B Fitting coefficient for e B/Ph . 
b(z,X) Total scattering coefficient. 
b(z, A 0l 0) Volume scattering coefficient. 
bb(z, A) Spectral backseat tering coefficient. 

6 l (A) Initial estimate of the particle scattering coefficient; 
used for determining the apparent particle scatter- 
ing coefficient for substances other than water. 
b p ( A) Total particle scattering. 

M A) Total Raman scattering coefficient. 
bw (A) Scattering coefficient for pure water. 

c(^,A) Spectral beam attenuation coefficient. 
e(z,660) Red beam attenuation (at 660 nm). 

c e (A) Corrected non- water beam attenuation coefficient. 
t: i(A) Initial estimate of the beam attenuation coefficient; 
used for determining the apparent beam attenua- 
tion coefficient for substances other than water. 
c p( A) Beam attenuation coefficient due to particles. 

C «’(A) Beam attenuation coefficient for pure water equal 
to a u ,(A) + b w ( A). 

C\ Polynomial regression factor. 

Additional polynomial regression factor. 

C f/ (A) AC-9 factory calibration coefficient. 

C r (A) Additional AC-9 factory calibration coefficient. 
CV(A) Digital response of reference detector. 

UV(A) Digital response of water transmission detector. 

I) Sequential day of the year. 
d Distance of lamp source from collector surface. 

P(X) I r radiance. 

B(z, m Smoothed (estimate of irradiance obtained by least- 
squares regression fit in the center of a depth inter- 
val. 

B a (X) Irradiance in air. 

Pa ii C calibration source irradiance. 

Ed(0“.A) Incident spectral irradiance. 

P\i(z, A) Normalized downwelled spectral irradiance 
Ed ( 2 , A ) Downwelled spectral irradiance 
Ps(X) Surface irradiance. 

Esky(A) Spectral sky irradiance. 

Psun(X) Spectral sun irradiance. 

T u ( 2 .A) Upwelled spectral irradiance. 

Ew(zi A) Irradiance in water. 

E’o(A) Mean extraterrestrial solar flux. 

Pi (A) Immersion correction factor. 

F v ( A) Fiekl-o f-view coefficient. 

/ Ratio of sensor-t.o-instrument. diameters. 

/(T) Offset voltage correction from the linear function 
characterizing temperature response. 

G(z.X) Solid angle dependence with water depth. 

g(7 ) Coefficient of a linear function characterizing tem- 
perature response. 

HI Matrix of coefficients h tJ . 

h tJ Analytic integral coefficients over the Hermitian pol- 
ynomials 7 *j . 

h(X) Normalized response function. 

K Vector of K n . 

K(X) Generic irradiance attenuation coefficient. 

K(z, A) Diffuse attenuation coefficient . 

Pd(z, A) Vertical attenuation coefficient for downwelled irra- 
diance. 
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Kl(z, A) Attenuation coefficient upwelled radiance. 

K„ K at node depth z n determined, with its vertical 
derivative by least-squares fit to radiometric pro- 
files. (See equations 37-41 and surrounding text.) 

K s (z, A') Apparent attenuation coefficient measured in a ho- 
mogenous water column. 

K u {z, A) Vertical attenuation coefficient for upwelled irradi- 
ance. 

k Fractional factor of total particle scattering. 

k f y/ tan 0o w • 

k c Wavelength independent fraction. 

I Cuvette pathlength. 

l s Nominal absorption pathlength. 

L(z,O t <t>) Submerged upwelled radiance distribution. 

L*(A,0,0) Atmospheric path radiance at flight altitude. 

Leal Calibration source radiance. 

I/sky (A, 0,0) Spectral sky radiance distribution. 

Lt Radiance at top of atmosphere. 

L u (z, A) Upwelled spectral radiance. 

L U (A) True upwelled spectral radiance. 

L u ( A) Measured upwelled spectral radiance. 

Lw{ A) Water- leaving radiance. 

Lwn( A) Normalized water-leaving radiance. 

Li(A) Apparent radiance response to a linearly polarized 
source. 

L 2 (A) Orthogonal apparent radiance response to a linearly 
polarized source. 

n Starting index in measurement for angular measure- 
ments, or node index for the integral K analysis. 

n 9 (A) Index of refraction of Plexiglas. 

riu,(A) Index of refraction of water. 

N Ending index in measurement sequence for angular 
measurements. 

OD b (A) Baseline optical density spectrum. 

ODg(A) Optical density of soluble material (Gelbstoff). 

ODp(A) Optical density spectra of filtered particles. 

ODr(A) Optical density reference for filter or distilled water. 

OD((A) Optical density of non-pigmented particulates (trip- 
ton). 

P(A) Polarization sensitivity. 

Pi Fitting coefficient for i = 1 to 5. 

Q( A) L u (0~ , A) to £ U (CT,A) relation factor (equal to 7 r 
for a Lambertian surface). 

~R Mean Earth-sun distance. 

R! l Reflectance from an uncalibrated radiometer. 

R l (z, A) Spectral reflectance. 

r Radius, or Earth-sun distance depending on usage. 

S( A r ) A coefficient of water temperature variation in 
a w (A ,T). 

S G ( A) Radiometer signal (uncalibrated) measured viewing 
a reflectance plaque. 

S 9ky Radiometer signal (uncalibrated) measured viewing 
the sky. 

SwW Radiometer signal (uncalibrated) measured viewing 
the water. 

T' Instrument temperature during calibration. 

T g ( A) Transmittance through a glass window. 

T S (X) Transmittance through the surface. 

T w ( A) Transmittance through a water path. 

t d (z, A) Downward spectral irradiance transmittance from 
flight altitude z to the surface. 


V Volume of water filtered. 

V True voltage. 

V Measured voltage. 

V(z) Transmissometer voltage. 

V^ ir Current transmissometer air calibration voltage. 

V{9) Normalized measured value for a cosine collector. 
V(0t) Mean normalized measured value of instrument re- 
sponse. 

V a i r Factory transmissometer air calibration voltage. 

Vaark Transmissometer dark response. 

We Weighting function. 

X Depth in meters. 

y Empirical factor. 

V Base 10 logarithm of the radiometric measurement 
Ed , E uy or L u - 

z Vertical coordinate. 

z ' Corrected depth for pressure transducer depth offset 
relative to a sensor. 

Zm Centered depth. 
z r Shallow depth. 

z s Exclusion depth due to data contamination. 

0 Filter absorption correction factor for scattering 
within the filter. 

j3(z, 0, A) Spectral volume scattering function. 

~P b The measured integral of the volume scattering func- 
tion in the backward direction. The symbol in the 
text — /3f,(0, A0, A, A A) — is for the measurement in 
the backward direction over the angle 0 at a given 
wavelength A. 

y t j (£) Hermitian cubic polynomial. 

A z Half-interval depth increment. 

A 0 Angular increment. 

6 Cosine response asymmetry. 

e Cosine collector response error. 
e S un Self-shading error for E su n- 
tgky Self-shading error for E s k y . 
e(A) 1 - e~ k ' aWr . 

6 Centroid angle of the scattering measurement. 

8 a In-air measurement angle. 

0 t Any nominal angle. 

On Angular terminus. 

0 n Angular origin. 

8t Tilt angle. 

0 U , In- water measurement angle. 

Oow Refracted solar zenith angle. 

01 Lower integration limit. 

02 Upper integration limit. 

k Self-shading coefficients. 

A' Channel of nominal wavelength. 

Ao Center wavelength. 

A n Any nominal wavelength. 

A r Near-IR wavelength. 

V'di z \ A) Spectral mean cosine for downwelling radiance at 
depth z. 

£(A) Minimum ship-shadow avoidance distance. 
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p( X) Bidirectional reflectance. 
p^(A) Gray card or plaque reflectance. 

a t Density anomaly. 

f Vector of measured optical depths. 

t( 2,A) Spectral optical depth. 
t 3 (A) Spectral solar atmospheric transmission. 

<j> Azimuth angle. 
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